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The Low Frequency Fundamental Bands of Methyl Chloride, Methyl Bromide and 
Methyl Iodide 


E. F. BARKER AND E. K. PLYLER, University of Michigan and University of North Carolina 
(Received May 1, 1935) 


The parallel type bands v; for CH;Cl, CH:Br and CHI have been found at \13.7, 16.4 and 
18.84, respectively, the rotation lines of the chloride band being partially resolved. The iodide 
band has a particularly sharp zero branch suitable for use in calibration. Both the chloride and 
bromide bands are complicated by isotopic effects, the separations being 6 cm™ and 1 cm“, 


HE infrared absorption spectra of the 
methyl halides out to 144 have been ob- 
served and correlated by Bennett and Meyer.' 
They were unable to reach the lowest frequency 
bands »v; for the bromide and the iodide. The 
corresponding chloride band was observed but 
not resolved, while that of the fluoride was com- 
pletely analyzed. Since these observations were 
made the apparatus has been improved some- 
what in sensitivity, and the NaCl fore prism has 
been replaced by one of KBr, extending the 
possible range of measurement out to beyond 
20u. In order to complete the series of methyl 
halide bands, and also to’ provide convenient 
reference wavelengths for the calibration of 
prism spectra in this region, we have observed 
the ». bands of CH;I and CH:Br. We have also 
remapped the CH;Cl band obtaining somewhat 
bette resolution than before. The results are 
indic.\ted in Fig. 1, where the ordinate is per- 
cent. “e absorption and the zero branches extend 
bout 80 percent in each band. 
CH;Cl band represents absorption by a 
n of gas 26 cm long at a pressure of 3 cm, 
it width being equivalent to 0.5 cm”. 


nett and Meyer, Phys. Rev. 32, 888 (1928). 


The rotation lines are partially resolved, es- 
pecially on the low frequency side. The irregulari- 
ties observed are not due to lack of resolving 
power, but to superposition of isotopic bands. 
The line spacing, which cannot be determined 
with great precision, is approximately 1.1 cm™!. 
This indicates a moment of inertia A about an 
axis through the centroid and normal to the 
symmetry line, of about 5010-” g cm’. The 
moment of inertia C with respect to the sym- 
metry axis is, of course, very nearly equal to that 
of CH,. The vibration »; to which this band is 
due consists primarily of a motion of the halogen 
atom toward and away from the CHs group. 
It is therefore particularly sensitive to isotopic 
differences in the mass of the halogen atom. In 
first approximation the methyl group may be 
treated as a rigid unit, and the reduced mass 
computed as for a diatomic molecule. The differ- 
ence in frequency thus predicted for the two 
bands due to Cl** and Cl*’ is 6 cm~. This agrees 
exactly with the observed difference between the 
principal zero branch at 733 cm™ and the minor 
one at 727 cm“. 

For the observations on methyl bromide a 
pressure of 5 cm was used in the 26 cm cell, and 
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Fic. 1. The parallel type absorption band »;. 


the slit included a spectral range of 0.3 cm7. 
The two isotopic bands are of nearly the same 
intensity in this case, and are only slightly 
separated, the zero branches appearing at 610 
and 611 cm™. The expected displacement for 
Br* from Br7*, computed as before is 1.2 cm—. 
The effect of partial superposition is to make the 
maxima appear somewhat closer together than 
they actually are. The interval between rota- 
tion lines should be approximately 0.5 cm-, 
which is very close to the limit for this spectrom- 
eter. However, it is clear that the bands super- 
pose in such a way that no resolution is to be 
expected. 

The third curve represents the absorption of 
methyl iodide at 20 cm pressure in the 26 cm 
cell. The zero branch is very sharp and intense, 


and constitutes an excellent reference point. 
It appears at 532.9 cm~. Although there is no 
complication due to isotopes the rotation lines 
cannot be distinguished since their separation 
is below the limit of resolution of the spectrom- 
eter. There is no indication of a transition be- 
ginning at the first excited state, which must 
have an appreciable population. This suggests 
that the anharmonic constants are small. Strong 
absorption due to water vapor in this region may 
account for some of the small irregularities in 
the band envelope. 

A comparison of the infrared observations 
with data obtained by Raman scattering has 
already been presented.? 


? Adel and Barker, J. Chem. Phys. 2, 627 (1934). 
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The Zeeman Effect of the Absorption Lines of Crystalline KCr(SO,).-12H.O 


F. H. SPEDDING AND G. C. NuttinG, Department of Chemistry, University of California, Berkeley 
(Received April 22, 1935) 


The Zeeman effect of the absorption lines of chrome alum 
crystals cooled to 20°K has been photographed for field 
strengths of from 12,000 to 25,000 gauss. Each of the three 
strong lines splits into three elliptically polarized compo- 
nents, whose separation of very close to 2w suggests ‘‘spin 
only” orientation in the field. At 50°K, but not at 20°, the 
pattern depends upon the orientation of the crystal 
lattice with respect to the magnetic field. Interaction be- 
tween closely spaced levels seems to be small and is evi- 
denced only by slight intensity dissymmetries. No shift 
of the central components from the positions of the field- 
free lines is observed, and the magnitude of a possible 
second-order Zeeman effect is estimated as less than 10 


percent that of the first-order effect. The selection rule for 
the magnetic quantum number, Am= +1, 0, seems to be 
obeyed, making it impossible for the present to decide 
whether the upper levels belong to the quartet or doublet 
system. Various possible coupling schemes are considered 
in the light of the experimental evidence and it is concluded 
that the actual one existing in chrome alum is either the 
strong field case in which L and S do not interact appre- 
ciably and L is almost completely quenched, or the very 
strong field case in which the coupling of the individual 
l’s is broken down and the individual orbital moments are 
quenched. 





ECENTLY there has developed consider- 

able interest in the nature of the energy 
levels of solids. Already a number of experi- 
mental and theoretical papers in this field have 
appeared and many more will appear in the near 
future. This subject is of the utmost importance, 
for the majority of the physical and chemical 
properties of a substance are intimately con- 
nected with the energy states of the outermost 
electrons of its atoms. 

The application of group theory to the wave 
mechanics has made it possible to treat the 
problem from the theoretical standpoint. The 
theoretical papers that have appeared so far can 
be divided into two groups. The first group deals 
with solids in which the lowest state is an ionized 
or metallic one,! and the second considers those in 
which the basic state of at least one of the ele- 
ments is atomic. Because the treatment is sim- 
pler, the former class has received by far the 
more attention, but discussion of it will not be 
pursued here since chrome alum is not a solid of 
this type. Among the papers in the second 
group is one by Bethe? in which he shows that 
the splitting of an energy level in the electric 
fieli of a crystal is determined primarily by the 
symmetry of the atomic pattern. Van Vleck,’ 


_ | by an ionized state we mean one in which an electron 
is not quantized with respect to any specific atom, yet 
remains in the crystal lattice. 

*hethe, Ann. d. Physik 3, 133 (1929). 

®Van Vleck, Electric and Magnetic Susceptibilities, 
Oxford, 1932. 


utilizing the work of Bethe, and making certain 
simple and reasonable assumptions as to the 
behavior of atoms in solids, has been able to 
explain the experimentally determined magnetic 
susceptibilities of salts of the iron group and the 
manner in which the susceptibilities depend 
upon the crystal structure. Schlapp and Penney* 


have applied Van Vleck’s theory to the specific 
case of chrome alum and have been able to 
account for the observed magnetic susceptibility 
very satisfactorily. At the same time their 
theory necessarily predicts the nature and 
position of the energy levels in the solid. Thus 
far the theoretical physicists have been greatly 
handicapped in their treatment of this problem 
by the meager amount of good experimental data 
upon which to base the model from which they 
start their calculations, and, naturally, this 
same lack of data has also prevented them from 
adequately checking their conclusions. 

For some years we have been interested in the 
nature of energy states in solids from the experi- 
mental standpoint. In a recent paper® we pointed 
out that the experimentally determined energy 
levels of KCr(SO,4)2-12H2O do not seem to be in 
agreement with the theories and calculations of 
Van Vleck, Schlapp and Penney. This paper is a 
continuation of the previous research presenting 
further evidence as to the real nature of the 
energy states in chrome alum and indicating 


4Schlapp and Penney, Phys. Rev. 42, 666 (1932). 
5 Spedding and Nutting, J. Chem. Phys. 2, 421 (1934), 
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where changes in Schlapp and Penney’s model 
will have to be made before their theories can be 
brought into accord with experiment. 


EXPERIMENTAL METHODS AND OBSERVATIONS 


The experimental methods have been de- 
scribed previously.’ Our earlier difficulty in cool- 
ing the chrome alum crystals to 20° when they 
were contained in a narrow Dewar flask has been 
overcome, partly by using thinner crystals (1.8 
mm thick) and partly by mounting them in a 
holder of improved design. The spectrograph 
had this time seven instead of five prisms and the 
dispersion was 3.8A per mm at 6700A. A Nicol 
prism was used in two positions, the one allow- 
ing passage of light polarized in a direction 
perpendicular to the direction of the magnetic 
field, the other light polarized parallel to the 
magnetic field. Various field strengths from 
12,000 to 25,000 gauss were used. 

The three absorption lines of chrome alum at 
14,932, 14,925 and 14,858 cm~ have absorption 
coefficients many times larger than those of the 
surrounding lines of the 6700A multiplet. With 
thick crystals the intensities of the absorption 
lines differ little from one another since all or 
nearly all of the incident light of the correspond- 
ing frequencies is absorbed. With thin crystals, 
however, while v14,932, 14,925 and 14,858 
attain complete absorption, the weaker ones 
do not and appear much fainter. Thus, by taking 
crystals of the proper thickness, the three 
strongly absorbing lines can be made to appear 
intense and at the same time the weaker neigh- 
bors at 14,935, 14,928.5 and 14,854 cm™' be 
made so faint that they can scarcely be seen. 
The latter then completely disappear into the 
background when a magnetic field is applied and 
the original intensity is spread over three or 
four components. One therefore observes Zeeman 
patterns of only the three intense lines, but it is 
not to be overlooked that the levels giving rise 
to the faint lines are still present and are still 
close to the levels giving rise to the strong lines. 
If these levels exert any influence upon one 
another it must be reflected in the Zeeman 
patterns observed. When a magnetic field is 
applied to the crystal cooled to 20°K, the strong 
lines split into triplets of approximately equal 
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separation, the separation being twice the norma! 
Larmor splitting. 

Photographs were made with the light ray 
and the magnetic field mutually perpendicular 
and passing along fourfold axes of the chrome 
alum octahedron; and with the light and field 
perpendicular, the light passing along a four- 
fold axis and the field directed perpendicular 
to a natural crystal edge. Schnetzler® has in- 
vestigated at liquid air temperature the Zeeman 
effect of the corresponding lines of potassium 
chromium selenate and he reports that there is 
dependence of the Zeeman pattern on _ the 
position of the lattice in the field. We have not 
observed such an effect in the lines of potassium 
chromium sulfate at 20°K.’ However, in a group 
of photographs of thicker crystals, 
temperature we estimate to be about 50°, the 
patterns observed are more complicated and 
more varied. One expects this, for in the absence 
of a magnetic field the lines are still diffuse be- 
cause of thermal agitation. In addition, the 
“high temperature”’ line’ v14,928.5 falling be- 
tween the two lines under observation, while 
absent at 20°, is at 50° almost as intense as these 
other two, and the Zeeman patterns of the three 
are to some extent superposed. One accordingly 
observes a pattern which is the resultant of all 
three lines. This with the temperature broaden- 
ing almost certainly prevents true resolution of 
the components, so that the pattern observed is 
in reality the sum of the envelopes of the three 
Zeeman patterns. This summed intensity itself 
gives rise to components which should show 
peculiar (nonlinear) behavior in fields of differ- 
ent strength. They should also show different 
behavior with rotation of the crystal axes in the 
field, since the relative intensity of the ‘“‘high-” 
and “low temperature” lines even at zero ficld 
strength varies for different crystal directions 
due to different polarization caused by the 
electric field of the crystal. It has actually been 
observed that, with the crystal at 50° and with 


whose 


6 Schnetzler, Ann. d. Physik 10 (3), 373 (1931). 

7 The low frequency component of the line v1 
appears in slightly different positions in the photogr :phs 
taken in the two directions of the magnetic field. Since 
this is the only difference greater than our probable : ‘ror 
of either the intensities or positions of the compon: "ts, 
and since this difference is only about twice the pro! :ble 
error, we have not attached much weight to it. 
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the field and light ray directed along four- 
fold axes, the composite pattern of v14,932, 
14,928.5 and 14,925 consists of four lines, and on 
changing the orientation of the light and the 
field with respect to the crystal the four become 
six lines. 

At 20°K as the magnetic field is applied the 
lines broaden, and, with the instrument used, 
first resolve at about 18,000 gauss. With 
I1= 24,800 the components are well separated 
from one another. At this latter field strength, 
even though the outer components of the strong 
lines must have crossed the outer components 
of the weak lines lying close to them, the Paschen- 
Back effect is extremely small. Since these lines 
are absorption lines of a solid at a very low 
temperature, the levels from which they arise 
lie equally close, yet the interaction between 
them must be very small. That there is some 
interaction can be seen from the fact that there is 
unsymmetrical distribution of intensity about 
the central components. This dissymmetry in 
intensity is appreciable for the line v14,858 but 
for the other two lines is small. There appears also 
to be slight dissymmetry in the position of the 
components, amounting to about 10 percent of 
the distance between them. This dissymmetry 
too is greatest in the line 14,858, but since its 
components are very weak, the errors in measure- 
ment are also largest here. In any case the amount 
of the dissymmetry cannot be stated surely 
because it is of the same order of magnitude as 
the possible error of measurement. 

No certain differences in the separation of the 
components of lines photographed with their 
light polarized parallel to the field and those 
with light perpendicular to the field have been 
detected. The intensities, however, are different. 
With the light polarized perpendicular to the 
field the unshifted components, corresponding 
to \m=0, are more intense than the outer 
components, for which Am=-+1; and of these 
latter the low frequency ones of v14,858 and 
14,935 and the high frequency one of 14,932 are 
the weaker. The low frequency component of 
v14.-58 is so faint that it was measured only with 
difficulty. With the light polarized parallel to 
the ‘eld the components are of about the same 
intensity, with the high frequency components 
bein. perhaps the strongest. 


OF CRYSTALLINE 


KCr(SO,4)2°12H20O 


Fic. 1. Enlarged reproductions of the 6700A multiplet 
of chrome alum. All temperatures are 20°K. 1. H=0; 
thickness of crystal, 4.2 mm. 2. H=18,700 gauss; thick- 
ness, 1.8 mm; s components. 3. H=24,850; thickness, 
1.8 mm; ~ components. 4. H = 24,850; thickness 1.8 mm; 
Ss components. 


No shift of the central components with re- 
spect to the field-free lines has been noted® and 
therefore any second-order Zeeman effect must 
be less than 10 percent of the first-order splitting 
and any relative shift less than 5 percent. The 
fact that the components arising from transi- 
tions in which Am=0 are more intense when the 
light is polarized perpendicular to the field is 
surprising and must be attributed to the pre- 
dominating influence of the electric field of the 
crystal. Since all the components can be seen 
with both directions of the Nicol prism, the 


8 Our earlier report (see reference 5) of the positions of 
the absorption lines of chrome alum at 20° in the absence 
of a magnetic field was incorrect. The values given for the 
6700A multiplet were too high by 1 cm™. 
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lines are elliptically polarized, as one might 
expect when electric and magnetic fields are 
crossed. 

Table I summarizes the observations on the 


TABLE I. 








Plate 1 
1 


2 3 4 
II + II 





14855.2f 
58.0f 
60.9f 

14922.7 


14855.2f 
57.8 
60. lof 

14922.4f 
25.0 
27.4 

14929.7 
32.0 
34.7f 


14855.1f 
57.5f 
60.4f 

14922.5 


14855.4f 
57.8 
60.30f 

14922.8f 
25.2 
27.4 

14929.8 
32.1 
34.8f 


24.9 
27.5 


25.0 
27.5 
14929.7 
32.1 
34.8 


14929.8 
32.1 
34.9 








spectra. In making plates 1 and 2 the mag- 
netic field was directed down a fourfold axis, 
in plates 3 and 4 perpendicular to a crystal 
edge. The letters f and vf indicate components 
which are faint and very faint, while 1 and | 
denote that the plane of polarization of the light 
ray is perpendicular and parallel, respectively, 
to the direction of the magnetic field. The 
probable errors in v and Av for the components of 
v14,858 are estimated at +0.2 cm, and for 
14,925 and 14,932 at +0.1 cm™. In a field of 
24,800 gauss twice normal splitting is 2.3 cm™. 


DISCUSSION OF RESULTS 


We shall discuss theoretically the Zeeman 
patterns which would arise from the various 
coupling which one might imagine to occur in 
a solid. We shall discuss these in light of the 
experimental evidence and point out wherein 
such types of coupling are supported by that 
evidence, and wherein they are unlikely or are 
completely ruled out. 

In our previous paper we have examined the 
various theories as to the nature of the sharp 
absorption lines of chrome alum. We concluded 
that they very probably arise from transitions 
between states originating from a 3d* electron 
configuration. From magnetic considerations it 
was evident that the basic state must arise from 
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a level belonging to the quartet system, and that 
the other low lying levels which were shown at 
that time to be present in chrome alum must also 
belong to this system. The experimental evi- 
dence, however, was not sufficient to determine 
whether or not the electric fields are strong 
enough to break down the 2/=Z coupling, and 
thus whether the lines are transitions from the 
4F to the *P, *H, ?G, *F, 27D and ?P terms or 
whether they are from a strong field ‘I term to 
another ‘I or 7’ term. It was shown that the 
weak field case in which J is still a good quantum 
number is excluded. 
Various possible methods of coupling are: 


Case I 


The coupling is similar to that which takes 
place in a gaseous ion and the effect of the elec- 
tric fields on the levels is either negligible or very 
small. This suggestion was put forward by 
Deutschbein® and is of course ruled out by the 
magnetic considerations. Bose and Datta! have 
suggested, however, that while most of the ions 
have their orbits bound so that only the spins 
can orient, a very few may be situated in pockets 
or holes in the solid in which they are free to 
orient like gaseous ions, and that it is these that 
give rise to the sharp absorption lines of chrome 
alum. Apart from the difficulty in conceiving of 
such “‘holes’’ in solids at the temperature of 
liquid hydrogen, the possibility is set aside for 
the following reason. If we grant that inter- 
combination jumps may be made freely and 
that the AJ selection rule is invalid, the lines at 
20°K could arise from transitions from a */’3/2 
term to a term such as *P5/2, ?G7/2, 7D5/2 and 
others. The g factors of all these terms are difler- 
ent, being 2/5 for 4F 3/2, 8/5 for *P5/2, 8/9 for 
°G7/2, and so on. The Zeeman patterns of the 
absorption lines would then be very complicated 
and very different from the ones observed, 
even though the levels are considerably altered 
by the electric fields of the lattice. 


Case II 


The electric fields of the crystal are sufficient 
to split the levels considerably, but the amount 
of the splitting is still small compared to the 


9 Deutschbein, Zeits. f. Physik 77, 489 (1932). 
10 Bose and Datta, Zeits. f. Physik 80, 376 (1933). 
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multiplet separations, so that J remains a good 
quantum number. Each original state is split 
into J+3 doubly degenerate levels.!! As soon 
as a magnetic field is applied there is a fixed 
reference direction in space, and each of these 
levels splits into two new levels. For the new 
levels there is a variety of g values and conse- 
quently the separation between pairs of them is 
different.” 

This case is ruled out for the same reasons as 
Case I. 


Case III 


The electric fields of the crystal are strong 
enough to break down the coupling between 
L and S so that J no longer has meaning. We 
may distinguish three variants of this case. 

(a) J is not a good quantum number, and 
while ZL and S orient separately, there is never- 
theless marked interaction between the two. 
The original level is split apart through the 
action of the electric field on the orbital moment. 
The magnetic levels are of course not split by 
the electric field directly, but rather by the 
magnetic field set up by the orbital moment 
precessing about the electric axis. In this case, 
except for a few specially degenerate levels, all 
levels are completely nondegenerate and hence 
cannot split further in a magnetic field, but can 
merely change their positions.” 

It is unlikely that this case is the correct one 
because of the magnetic behavior of the chrome 
alum crystal, and because all the absorption lines 
observed split in a magnetic field," yet do not 


"Kramers (Proc. Amst. Acad. 33, 959 (1930)) has 
shown that for atoms containing an odd number of equiva- 
lent electrons there is in crystal fields no distinction 
between states with the same absolute value of m. Stated 
loosely this means that there is no distinction between 
electrons revolving clockwise and anticlockwise in their 
orbits. In fields of high symmetry two or more of these 
J+‘ levels may combine to give a level of still higher 
degeneracy. 

' Penney and Schlapp have treated such a case in their 
theory of the magnetic susceptibilities of rare earth salts. 

'* Mulliken, Rev. Mod. Phys. 2, 60 (1930). 

‘\f the g factors of the high excited states were zero, 
trip|ts would also be observed on account of the multiple 
char cter of the basic state, but the center of gravity 
woul! be shifted by 1w. In this instance all the dissym- 
metres in the patterns would have to be due to the lowest 
State. and since all the lines originate at the lowest state 
the cissymmetries would perforce be the same for all the 
patt: rns. It is observed that this is not the case. While the 
Separ:tion dissymmetries may be attributed to experi- 
ment.l error, the intensity dissymmetries cannot. Further, 
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show any shift of the center of the Zeeman 
pattern from the field-free line. 

(b) The theory of Van Vleck and Schlapp and 
Penney. The basic *F term is split by a cubic 
crystal field into three levels, two triply degener- 
ate and one, the lowest, a single level. The 
separation of the highest and lowest levels is 
about nine volts. A much weaker rhombic field 
causes the two triply degenerate levels to be- 
come nondegenerate. In a magnetic field all the 
levels split further into four components. Since 
the lowest level is far removed from the others, 
its interaction with them is small, and since it is 
already nondegenerate with respect to the orbital 
moment, only the spin moment orients in a 
magnetic field, thus giving a g factor of 2. 
The excited levels, on the other hand, interact 
with one another strongly and therefore do not 
have g factors of 2. They also show an appreci- 
able second-order Zeeman effect. Thus the 
magnetic susceptibility of chrome alum, whose 
Curie constant corresponds to a g factor of 2 
arises from a very special case of the energy 
levels of solids. 

This theory seems highly improbable for the 
following reasons. Firstly, there are actually 
several low-lying levels belonging to the quartet 
system, and the theory permits only one. And 
secondly, the g factor of the excited levels as 
well as the low levels must be very close to 2. 
In the theory a g factor of 2 can arise only in the 
special case in which there is no strong inter- 
action with a neighboring level, and it seems very 
unlikely that all the levels observed should be of 
this unusual type. Furthermore, although the 
levels occur close together, in some cases being 
so close that when split by a magnetic field their 
components actually cross, there is no inter- 
action. Any second-order Zeeman effect or inter- 
action that exists must certainly be less than 
10 percent of the first-order Zeeman splitting. 

(c) There is slight or no interaction between L 
and S. Here the electric field of the crystal in 
splitting the levels apart has so completely 


if these upper levels had g factors of zero because of non- 
degeneracy they should show a second-order Zeeman effect. 
None is observed. Finally it seems unreasonable to suppose 
that the excited states should be nondegenerate and the 
basic states not, particularly when one considers that in 
the ions of the iron group there is a variety of lowest 
levels, and in almost every case they are known to split. 
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quenched the orbital moment that it does not 
orient in attainable magnetic fields. Besides, the 
electric field has so tightly bound the electron 
in its orbital potential trough that magnetic 
fields are not able to shift the positions of the 
levels appreciably. In other words, there is 
either a very small second-order Zeeman effect 
or none at all. This case or Case IVb seems most 
probable in view of the experimental evidence. 


Case IV 


The very strong field case. The crystal electric 
fieids are strong enough to produce splitting 
which is large compared to the distance between 
multiplets, and they thus break down the coup- 
ling between the individual electrons, Y/=L. 
(a) Strong interaction between 7 and s. Fhis 
case is ruled out for the same reasons given in 
IIIa. (6) Little or no interaction. Except for 
different positions of the levels this case cannot 
be distinguished from IIIc. 

We feel that in the present state of our knowl- 
edge the experimental data are insufficient to 
distinguish between the strong field Case IIIc and 
the very strong field Case IVb. We feel also that 
we cannot yet determine whether the excited 
states belong to the doublet or to the quartet 
system. However, we might point out that in the 
cubic crystal of chrome alum the magnetic 
selection rule Am=0, +1 is obeyed, and that if 
transitions of the type for which Am=-+2, 
+3 occur they must be very rare. In the absorp- 
tion spectrum of gadolinium (Gd***) ion in a 
monoclinic crystal the selection rule Am= +4 is 
known to hold, and such transitions are very 
frequent. Therefore, when the selection rules 
for cubic fields are known (and we have work 
under way at present in which we hope to 
determine them) this distinction may perhaps be 
made; for if jumps for which Am= +3 are per- 
mitted, the upper states must be doublets to 
account for their absence. 

In the unit cell of chrome alum there are four 
Cr*++*+ ions, each surrounded by six molecules of 
water. After the water molecules the nearest 
neighbors of a chromium ion are two SO;- ions, 
placed symmetrically about it, which produce 
an electric field directed along the line joining the 
ions. The three ions lie on a trigonal axis of the 
cell, and the trigonal axes for the four chromium 
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ions point in four different directions. In zero 
field or so long as spin-orbit interaction is so 
small that the spin can orient freely, the absorp- 
tion lines arising from these four types of chro- 
mium ions are superposed. On the other hand, 
if the orbital moments orient in a magnetic 
field or if they influence the orientation of the 
spin moments, the positions of the lines in the 
presence of a magnetic field are different, for 
these positions are determined by the resultant 
of the magnetic and electric fields. The resultant 
will be different for each of the four types of 
chromium ions, and four Zeeman patterns which 
do not perfectly superpose but instead yield a 
very complex and difficultly resolvable spectrum 
will be produced. Thus the fact that we obtain 
fairly simple and sharp Zeeman patterns is in 
itself evidence that only the spin orients and 
that interaction terms must be small. 

Lehmann” has investigated the Zeeman effect 
of the absorption lines of ruby at liquid air 
temperature and he reports that the Zeeman 
pattern both as to the separation and intensity 
of its components is dependent upon the orienta- 
tion of the lattice in the magnetic field. If this 
effect is not dye to the presence of another level, 
which might or might not fade out when the 
temperature is further lowered, the orbital 
moment must contribute to some extent in 
ruby. It is interesting to note that the magnetic 
behavior of CrsO3, which is the coloring matter 
in ruby, is decidedly anomalous and that the 
susceptibility at low temperatures shows <e- 
partures from Curie’s law in the same direction 
as for the free Cr*++*+ ion, in which the orbit does 
contribute largely to the magnetic moment. 
One might think from the two cases of chromic 
oxide and chrome alum that the basic state of an 
ion in a crystal is not intrinsically different from 
the excited states; it is not a special case. -\c- 
cordingly, if the magnetic susceptibility shows 
that the orbit is contributing to the magn«tic 
moment, the Zeeman pattern of the absorp'ion 
lines will show the same, and if the magnetic 
susceptibility shows that the orbital momen: is 
completely frozen out, the Zeeman pattern will 
indicate spin only orientation. 


15 Lehmann, Ann. d. Physik 19, 99 (1934). 
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In concluding, we feel that the experimental 
evidence is not in agreement with the theories of 
Van Vleck and Schlapp and Penney and that 
their basic model will have to be changed to 
conform with the following experimental evi- 
dence from chrome alum. The energy levels 
appear in groups of closely spaced levels, the 
lowest group belonging to the quartet system. 
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These levels split in a magnetic field and possess 
a g factor very close to 2. Any interaction be- 
tween the levels, if it exists, must be very small, 
of the order of magnitude of 10 percent or less 
of the splitting. In addition, these levels cannot 
show a second-order Zeeman effect of more than 
5 percent of the normal Zeeman splitting in 
fields of 25,000 gauss. 
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Electronic Structures of Polyatomic Molecules and Valence 
VI. On the Method of Molecular Orbitals 


RoBert S. MULLIKEN, Ryerson Physical Laboratory, University of Chicago 
(Received April 18, 1935) 


The use of the method of nonlocalized molecular orbitals in building up a conceptual scheme 
or qualitative theoretical framework into which empirical data (chemical and spectroscopic) 
can be fitted, is emphasized. This should be distinguished from the use of the method, often 
with rough ‘‘LCAO” approximations (linear combinations of atomic orbitals), in semiquanti- 


tative calculations. 


N earlier papers of the present series and 
elsewhere! an attempt is made to describe and 


interpret electronic structures of molecules in 
terms of electron configurations composed, so far 
as reasonably possible, of molecular orbitals. 
Since the object of this work seems not to have 
been made entirely clear, it seems desirable to 
restate it briefly. 

The whole matter turns about the meaning 
and use of the idea of electron configurations. 
Ever since Bohr succeeded in explaining the 
main features of atomic structures and spectra 
and of the periodic system, electron configura- 
tions have been in constant use for describing 
atomic structures approximately. Each electron 
is assigned to a definite orbit (or orbital, in the 
quantum mechanics), having a fairly definite 
energy or term value, at least for an atom in a 


‘This series: R. S. Mulliken, Phys. Rev. 40, 55; 41, 
49, 751 (1932); 43, 279 (1933): I, II, III, IV; J. Chem. 
Phys. 1, 492 (1933): V. Also R. S. Mulliken, Chem. Rev. 
9, 347 (1931); Phys. Rev. 43, 765 (1933); 47, 413 (1935). 
Also J. E. Lennard-Jones, Trans. Faraday Soc. 30, 70 
(1934); in Table III and the text of this paper, the desig- 
nations 7; and J. have been interchanged as compared 
with IV of this series, evidently due to a misunderstanding 
caused by excess of brevity in the writer’s discussion of 
electric moment components in IV. Also papers of Hund, 
Herzberg, Hiickel and others. 


definite state of ionization. It must be agreed, to 
be sure, that this term value varies somewhat 
depending on the state of the atom-as-a-whole, 
e.g., 3p has somewhat different term values in 
-++353p, 3P and ---3s3p, 'P. 

The quantum mechanics, of course, shows that 
the concept of an electron configuration is in 
final analysis incompatible with an exact correct 
description of the electronic structure of an atom. 
Nevertheless it shows, particularly when Har- 
tree’s idea of the self-consistent field is intro- 
duced, that one can still in ordinary cases use a 
unique electron configuration as a basis for a 
simple and relatively good approximation. (Wave 
function of atom =antisymmetrical linear com- 
bination of permutations of a definite electron 
configuration assignment.) 

With molecules, matters become more com- 
plicated. Because of the less simple force field, 
the possibility of specifying a unique electron 
configuration tends more easily to be lost. In 
other words, the inexactness of the electron 
configuration idea is sufficiently great here to 
admit the possibility of two or more alternative 
electron configurations each of which may form 
the basis for a fairly good approximation to the 
complete wave function. In such a situation, as 





376 ROBERT S. 
Slater has pointed out,? we are forced to adopt 
views of molecular electronic structures which 
are more or less intermediate between those 
suggested by the alternative configurations. In 
this way we get fairly definite although often 
complicated pictures of molecular structures. 
Of course, we might abandon altogether the 
attempt to use electron configuration ideology; 
but the time seems not yet ripe for this. In 
fact, the idea of electron configurations will 
probably always be valuable as a conceptual aid 
and in classification. 

Methods of assigning molecular electron con- 
figurations may be divided into two types. In one 
of these atomic orbitals are used exclusively; in 
the other, molecular orbitals of some sort are used 
for shared electrons. The first method follows 
the ideology of chemistry and treats every 
molecule, so far as possible, as composed of 
definite atoms or ions. The electron configura- 
tion is then the sum of the configurations of these 
atoms or ions, e.g., the normal state of Lig may 
be represented as 


15°2s-2s1s?, 12+). 


(1) 


So far as homopolar binding is concerned, this 
method, which grew out of the work of Heitler 
and London on He, was developed by Pauling 
and Slater. It has had notable success as a 
qualitative conceptual scheme for interpreting 
and explaining empirical rules of valence and in 
semiquantitative, mostly semiempirical calcula- 
tions of energies of formation and other proper- 
ties of molecular electronic states, particularly 
the normal states. 

Departing from chemical ideology, the second 
method treats each molecule, so far as possible, 
as a unit (cf. I and II of this series'). This method 
was introduced by Hund for diatomic molecules, 
and was developed by Hund, the writer, and 
others. In its pure form, the method goes too far, 
as was pointed out by Lennard-Jones: it is best 
to use molecular orbitals only for shared electrons 
and atomic orbitals for the others. Even then, 
the type of orbitals to be used may often be 
chosen in various ways: one may use exclusively 
two-center (‘“‘localized”’) orbitals even for poly- 

2 J. C. Slater, Phys. Rev. 41, 255 (1932). 

3 J. E. Lennard-Jones, Trans. Faraday Soc. 25, 668 
(1929). The method could thus now more accurately be 
called the ‘‘method of atomic and molecular orbitals,’”’ but 


it is the use of molecular orbitals which constitutes the dis- 
tinguishing and vital feature of the method. 
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atomic molecules (Hund), or one may use many- 
center (“‘nonlocalized’’) orbitals for the latter 
(present method). These possibilities of course 
reduce to one in the case of diatomic molecules. 
As an example, the normal state of Lis appears as 


15?(0,2s)71s?, !Z+,, 
or 15?(2s+-2s, o,)"1s?, 12+,, 


where o,2s and (2s+2s, o,) are two types of 
shorthand for representing the bonding molec- 
ular orbital.4 

It is the writer’s belief that, of the various 
possible methods, the present one may be the 
best adapted to the construction of an explora- 
tory conceptual scheme within whose framework 
may be fitted both chemical knowledge and data 
on electron levels from molecular spectra. A 
procedure adapted to a broad survey and inter- 
pretation of observed relations is here aimed at, 
rather than (at first) one for quantitative calcu- 
lation, which logically would follow later. Given 
an observed molecule or ion of known shape and 
size, what is its electronic structure in terms of 
an electron configuration using, in general, non- 
localized orbitals for shared electrons? What is 
the relation of this structure to the molecule’s 
spectrum, to its ions, and to the structures of 
other similar and dissimilar molecules? This 
procedure has already yielded good results for 
diatomic molecules,‘ and it is now a question of 
extending it to polyatomic molecules.! To a 
considerable extent the present method and ob- 
jective are very analogous to those used by Bohr 
in developing the theory of atomic structure.’ 

Various limitations and qualifications must of 
course be recognized. In particular the usefulness 
of the present method is confined to cases where 
internuclear distances (r) are fairly small, as in 
stable molecules. In other cases, or in studying 
processes of dissociation, the present method 
must be supplemented by comparison with that 
using atomic orbitals alone. Even for stable 
molecules, comparisons should always be made 
with the latter method, if maximum insight is to 
be secured.’ 

4Cf., e.g., R.S. Mulliken, Rev. Mod. Phys. 4, 1 (1932); 
see sections A3, A8, D1 in regard to notation; for energy 
levels of Hz, cf. p. 46. Cf. also Phys. Rev. 46, 549 (1934), 
on new notation. P 

5 This was noted in the early papers on diatomic 
molecules: R. S. Mulliken, Phys. Rev. 32, 186 (1928);. 


F. Hund, Zeits. f. Physik 51, 579 (1928); cf. also G. Her2- 
berg, Zeits. f. Physik 57, 601 (1929); J. E. Lennard-Jones. 
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As a simple illustration of the difference be- 
tween the use of atomic and of molecular orbitals, 
we may consider the three lowest electron levels 
of the He molecule, ali taken for an inter- 
nuclear separation 0.7410-* cm equal to the 


equilibrium separation (r.) of the normal state ~ 


of He. In terms of the probable best simple 
electron configuration in each case, these three 
states appear, by using atomic orbitals, as 


1s-1s, 1Z+,; 1s-1s, *Z+,; 15%, 'E+,, (3) 


the third configuration corresponding to H~Ht+ 
or H+tH~-. Using molecular orbitals, we have 


(1s0,)?, 12+,; 1sog2pou, *Zty; 1sog2poy, 'Zty. (4) 


” 


Here ‘‘1sc,”’ is often written as (1s+1s, o,) or 
gis. The molecular ionization potential, for a 
process without change of the internuclear 
distance, has approximately the respective values 
16.3, 6.6 and 3.9 volts for the above states 
zt, Zt, and ‘Z*.. 

From the atomic orbital viewpoint (3), the 
values 16.3 and 6.6 must be interpreted as 1s 
term values for the H atom in two states of 
combination, and 4 volts as the term value of 1s 
in 1s? of H~ when combined with H*. The 1s 
term value in uncombined H is 13.53 volts, in 
H~ it is 0.715 volt. The great divergence between 
the term values 16.3 and 6.6 is an indication of 
the poorness of the atomic orbital approxima- 
tion in this case, since a good electron configura- 
tion should yield only moderate variations in 
term value when the state of coupling of the 
removed electron is altered. If we use the molec- 
ular orbital standpoint (4), we find the term 
value 16.3 volts for 1sc, from 'Z+,, and the 
values 6.6 and 3.9 volts for 2po, from *>+, and 
'Z*.,, respectively. It will be seen that according 
to the criterion of constancy of the term values of 
any orbital, the molecular orbital configurations 
represent (for small r values) the better approxi- 
mation in this case. 

Recent quantum-mechanical calculations of 
James and Coolidge® on Hz and Liz suggest that 
the use of carefully constructed molecular or- 
bitals for shared electrons may often give better 


°H. M. James and A. S. Coolidge, J. Chem. Phys. 1, 
825 (1933); H. M. James, J. Chem. Phys. 2, 794 (1934); 
3, 9 (1935). The authors for the most part did not use 
pure molecular orbital configurations, but their work 
nevertheless indicates that these would give rather good 
approximations. 
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quantitative results even in calculating dissocia- 
tion energies and other molecular constants than 
do methods using only atomic orbitals. Pre- 
viously it had been argued that the reverse is the 
case.’ This argument and its conclusion, how- 
ever, really never referred to the use of molec- 
ular orbitals as such, but only to a particular 
and decidedly rough type of approximation in 
which each molecular orbital is expressed as a 
very simple linear combination of atomic orbitals 
(“LCAO molecular orbitals,” let us call them). 
[ Example: the LCAO approximation for o,2s of 
Lig is a (254+2s2), where A, B refer to the two 
Li atoms. | 

A similar statement applies in part to the 
criticism’ that the use of molecular orbitals gives 
over-ionic wave functions (e.g., too much HtH- 
character for Hz). When a molecular orbital is 
occupied by a pair of electrons, overionicness 
makes the latter spend too much of their time 
together behind the nuclei which are connected by 
the orbital. It is probable that the overionicness 
resulting from the use of LCAO approximate 
molecular orbitals is in general much reduced 
when better approximations are used for the 
orbitals. This is for example strikingly brought 
out by the work of James on Lie, which shows 
that when a better approximation is used for 
the bonding molecular orbital, the fraction of 
time spent by an electron behind the Li atoms 
as compared with that spent between them is 
much less than for an LCAO approximation 
orbital. 

While the LCAO type of approximation is 


7 Cf., e.g., G. W. Wheland, J. Chem. Phys. 2, 474 (1934), 
who, however, favors the use of molecular orbitals in a 
_—_— way; and cf. some of the references given by 

im. 

8 In early papers,® Hund and the writer used for diatomic 
molecular orbitals a notation (2p0, 3sc, 3dz, etc.) corre- 
sponding to a correlation with the united atom. Lennard- 
Jones, on the other hand,’ emphasized a point of view and 
introduced a notation (2p0, 2po’, etc., later replaced by 
Hund or the writer with o,2p, o.2p or o2p, o*2p, later 
also the more explicit forms (¢+0,¢,), (¢—o, ¢.)) cor- 
responding to the LCAO approximation. Similar notations 
of both types have been introduced! for polyatomic cases. 
It has never been intended that any of the various nota- 
tions‘: ! should be interpreted too literally, since the real 
molecular orbitals are not in general closely approximated 
by either united-atom or LCAO approximations. If used 
cautiously with this understanding, however, the various 
notations can be more helpfully suggestive than mis- 
leading. In most cases, the LCAO notations are decidedly 
the more appropriate,? but in many cases, especially for 
excited orbitals, the united-atom or related notations are 
the better. 

( 033) especially J. H. Van Vleck, J. Chem. Phys. 1, 177 
1933). 





378 ROBERT S. 
very simple and convenient as a qualitative 
guide, it is in no way an essential part of the 
present method. Especially is it not essential to 
the qualitative conceptual scheme of the latter.’ 
Nevertheless, it was through the use of LCAO 
forms that a clear understanding of the bonding 
properties of molecular orbitals was first ob- 
tained. This was brought to the front especially 
by a paper of Lennard-Jones,’ led up to by work 
of Herzberg.® In subsequent work, LCAO forms 
have been found reliable guides to the bonding 
or antibonding properties of orbitals. 

The arguments as to the poor quantitative 
results obtained with molecular orbitals were 
based also on the idea that detailed interactions 
between electrons (analogous to those which 
give rise to atomic levels with different LZ and S 
values in atoms) cannot feasibly be taken care of 
in this method.’ Although it is true that the 
inclusion of such interactions would make 


quantitative calculations very complicated if one 
should attempt a completely theoretical calcula- 
tion from the beginning, their qualitative inclu- 
sion has always formed a vital part of the 
method of molecular orbitals used as a conceptual 


scheme for the interpretation of empirical data 
on electronic states of molecules. In this case, 
the term values of molecular orbitals are obtained 
largely empirically (from band spectra or ioniza- 
tion potential data), and it then remains as a 
theoretical problem only to estimate the de- 
tailed interactions between electrons. For in- 
stance, it is precisely because of the inclusion of 
such interactions that the present method was 
able to offer the first, and a very simple and 
satisfactory, qualitative explanation of the 
paramagnetism of oxygen gas. The theory of the 
electronic states, including the qualitative energy 
order, for specified electron configurations of 
molecular orbitals in diatomic molecules is rather 
well developed,’ comparing favorably with early 
stages in the theory of the energy states of elec- 
tron configurations of atomic orbitals in atoms. 

In the development of the theory of atomic 
structure and spectra, the concepts of orbits and 
later orbitals were long used freely in correlating 
empirical data without any necessity of being 
able always to obtain from pure theory the 
exact forms of the “‘real’”’ or ‘‘best’’ (Hartree 
self-consistent field type) orbits or orbitals, or the 
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energy order of the various states whose existence 
was predicted on the basis of combinations of 
quantum vectors. These various forms and 
interactions were, in fact, most easily semi- 
quantitatively determined by fitting the em- 
pirical data into the qualitative theoretical 
framework. 

In developing a similar qualitative theory for 
molecules, the use of LCAO forms for molecular 
orbitals, although not fundamental, is neverthe- 
less convenient, as has been mentioned, since 
these are capable of correctly representing sym- 
metry properties and other main qualitative 
features of the “real” or “best” orbitals. It is, 
however, these real or perfect, even if not ex- 
actly known, orbitals which are used concep- 
tually in the interpretation of empirical data. 

In the case of diatomic molecules, more or less 
nearly perfect “‘real’’ orbitals have long been 
known for the one-electron molecule H+; here 
elliptical coordinates are appropriate. It has long 
been realized that fairly good approximations 
might also be obtainable in terms of elliptic 
coordinates for the orbitals of many-electron 
diatomic molecules. It remained, however, for 
James and Coolidge® (cf. especially the last 
paper of James) to extend this method so as to 
get fairly good approximations. James makes use 
of power series with coefficients determined so 
as to minimize the energy. While the use of 
series in ordinary elliptic coordinates has thus 
given good results; one hopes that the use of other 
or more general types of coordinates may give 
even better molecular orbitals. 

It appears reasonable to expect that, as time 
goes on, analogous methods will be developed 
for polyatomic molecules, using expansions 
involving coordinates adapted to the symmetry 
of the molecule or of the particular orbital. 
For example, the nonlocalized molecular orbitals 
used for CH, in the present method! could 
probably best be developed as series using 
spherical polar coordinates. 

In connection with any electron configuration 
method, it must of course always be remembered 
that the orbitals used have a meaning only in 
relation to the electron configuration type of 
approximation, and only in relation to the 
particular way chosen for assigning an electron 
configuration. 
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The Rotational Entropy of Nonrigid Polyatomic Molecules 


MAXWELL L. ErpinorF ANpD J. G. Aston, School of Chemistry and Physics, The Pennsylvania State College 
(Received April 4, 1935) 


An expression is derived for the rotational entropy of nonrigid polyatomic molecules with 
zero potential and apart from nuclear spin and symmetry number corrections. Previous em- 
pirical formulae are derived as special simplified cases of the general formula. Applications are 
made to several important types of molecular models. 





I. INTRODUCTION 


HE theoretical calculation of the rotational 

entropy of molecules, apart from nuclear 
spin and symmetry number corrections, de- 
pends on the evaluation of the partition function! 
Sabg e+’ Leta, 8, -**,4/*” where the summation 
is extended over all values of all the quantum 
numbers a, 8,--+, \ and where é€q, g, ---, , are the 
corresponding energy levels, obtained as the 
characteristic values of the quantum-mechanical 
rotational problem. In the field of nonrigid 
molecules, these levels have been solved exactly 
for the ethane top? and for spherically sym- 
metric tops such as tetramethylmethane.’ No 
solutions of nonrigid asymmetric tops have been 
obtained. An empirical formula has been pre- 
sented and tested by Halford.‘ This formula, 
although useful, cannot be applied without 
uncertainty (Section III). Expressions will be 
derived for this sum, assuming limiting classical 
rotation. These formulae, together with vibra- 
tion spectrum data, make possible comparisons 
of theoretical entropies with those obtained from 
heat capacity measurements down to 11°K, 
which are being conducted in this laboratory and 
elsewhere. Even at very low temperatures these 
expressions agree very well with the summations 
over the quantum levels (Section IV). In the 
application of the expressions to be derived to a 
specific model, it is unnecessary to obtain the 
Hamiltonian form of the rotational energy, and 
to perform the laborious integration over the 
momentum space associated with the rotational 
assembly. 


1 Fowler, Statistical Mechanics, Chapter II. 

ass Brunauer and Mayer, J. Am. Chem. Soc. 55, 37 
* Lacoste, Phys. Rev. 46, 718 (1934). 
* Halford, J. Chem. Phys. 2, 694 (1934). 


II. EXPRESSIONS FOR THE ROTATIONAL PARTI- 
TION FUNCTION AND ROTATIONAL 
ENTROPY 


The molecule under consideration is an M 
particle body, consisting of a rigid skeletal 
frame to which are attached N rotating (zero 
potential) tops, which, apart from their spins 
relative to the frame, have angular velocity 
components due to the motion of the frame. The 
center of gravity O of the molecule is considered 
fixed. The kinetic energy of rotation about the 
center of gravity may be derived in different 
ways depending on the specific nature of the 
problem. A general one is: The usual Eulerian 
angles (see Fig. 1), 6, ¥, ¢, refer three | axes, 
origin at O, and fixed in the frame to three refer- 
ence axes fixed in space. Angular coordinates, 
a1, @2,***, a, are used to define the orientations 
of the tops relative to the frame. Components of 
angular velocity due to rotating frame and 
spinning tops are resolved along WN sets of three 
L axes (xj, yi, 2; for the 7th top), which are the 
principal axes for the N tops at O. If I2., Lyi, Ji 
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are the moments of inertia, and wi, wyi, w2i are Ai;=fi; (x1, x2, --* x,).* If the matrix is non- 
the corresponding angular velocities about these singular, the transformation to Hamiltonian 
axes, (and similarly for the other tops), then the form gives: 
rotational kinetic energy may be written: 

‘: @;=2K.E.= , B;P:P; (3) 
$,=K.E.=3(D (Tews? t+ Tywy? tLe?) }. (1) “ies 

i=l 

of matrix (B) and determinant [B] where 

The number of independent degrees of freedom Pi =0K.E./0%;; 
will be determined by the original constraints of 
the problem. If the number of these be 7, and Bij= (cofactor of Aij)/LA ]=fu' (x1 + +xr) § 


corresponding coordinates x;, %2,°** x,, and (B)X(A)=(1); [B]X[A]=1. 
velocities #1, %2, +++ &,, it is possible to bring ®, to 
the positive definite symmetric quadratic form, It is necessary to integrate over the momentum 


and configuration space the following definite 
integral, which is the limiting value approached 
by Dale +++ Dea, 8, °-+,d/*” with increasing 
temperature and especially rapidly for molecules 
of matix (A) and determinant [A] where’ with large moments of inertia: 


1 +00 +co Br Bi ‘s 
o=— oe f f ee f get 3 Se »-dpdx,:++dx, (4) 
¢ —oo —oo ar a 


where h is Planck’s constant, k is the Boltzmann constant, x; varies from a; to B;, p; from — » to+~. 
It will be proven that in general, the integral with respect to the momenta is: (2rkT)*"A ]!/h’. 
[A ] is readily obtained from (2). The problem is then reduced’ to the evaluation of: 


@,=2K.E.= Zz A i454; (2) 


i, j=1 





Br B1 
Q= [(ankryr3/ar} [ . f [A ]'dx,-+-dx,. (5) 


This expression may be integrated to furnish the complete solution of the problem. Examples will 
be given in IV. 

Proof: 

. 
Let $;= > B,;P,P; be a real positive definite quadratic form with nonsingular matrix (B) and 
i,j=1 

determinant [B]. If the P’s are subjected to a homogeneous linear transformation with the matrix 
(C) the new quadratic form in, say, the P’’s, has the matrix (C)’ (B) (C) where (C)’ is the transposed 
matrix of (C). In particular @; may be reduced by a real orthogonal transformation in r variables of 


Tr 
° 2 - . 
determinant 1 to the normal form: > D,P;’. Because of the positive definiteness of the form, none 
i=1 


of the D,’s may be identically zero.** The Jacobian of the transformation is 1. The integrand 1s 
positive and the limits of integration are from — # to + «. (4) becomes: 


1 +co +co Br B1 ‘ 
a f én f f on f e-12kT 2 DiPi*gP,!. --dP,'dx ++ +dx,. (4a) 
hr —o —a ar ay ‘ 


The integration with respect to the momenta gives at once: 


* The symbols for a matrix ( ), and determinant [ ], are stated explicitly here. Reference: Bocher’s Higher 
enclosing capital letters are adhered to throughout. Algebra, or Hilton’s Linear Substitutions. 
** The theorems from the theory of real quadratic forms 
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(2rkT)"!? be Bi 1 b 
ge ff (—.-—~) die, +-dx,. (4b) 
hr J ar eS ay D,Dsz: ? -D, 


D, De - 
becomes: 


-- D,=LC]’ [B] [C]=[B], since for a real orthogonal matrix [C ]’/=[[C }*=1. Thus (4b) 


Br Bi 
Q=[(2akT)"!2/h7] f ++ | [A]Mdxy-+-dx,. (5) 


If the coefficient of the integral in (5) is called 
\ and the integral itself A, then the rotational 
partition function 


A (6a) 


and the rotational contribution to the entropy 
(for the gaseous state and excluding nuclear 
spin and symmetry number corrections), is, by 
well-known methods: 


Spr=RIn {dAer!?}, (6b) 


III. A SpecraAL CAsE; DERIVATION OF A PRE- 
VIOUSLY SUGGESTED EMPIRICAL FORMULA 


In the special case that [A ] has a form closely 
Tr Tr 

related to sin? 611 J; where II J; is the product of 
i=1 i=l 


moments of inertia and where 6, the angle of 
inclination of the frame, varies from 0 to z, 
the other angular coordinates from 0 to 27, then 
(5) may be integrated at once to give* 


- } 
Of =D M-IETY IY itr} /hr (7) 
i=1 


and Sp = R In e"/?Q’ which is exactly the empirical 
formula proposed and tested by Halford.‘ It is 
now possible to discuss this formula theoretically: 
(a) There is an uncertainty in the correct mo- 
ments® to be substituted in (7), because of the 
nonrigidity of the molecule. The correct mo- 
ments are always in [A] of (6a, b). (b) Be- 
cause of the free rotation of the tops, the so-called 


* An example of the special case for which [A] has the 
Tr 
form sin? @]] J; is the asymmetric top of moments of 


apes t=1 
inertia A, B, C, for which [A ]=ABC sin? 6. Formula (7) 
will apply because formulae (5), (6a, 6b) reduce to (7) in 
this case. 

* Kassel, J. Chem. Phys. 3, 115 (1935). 





“reduced’’ moments* must be used. In certain 
cases, the magnitude of these reduced moments 
may be easily found. In general, however, the 
proper reduced moment will not appear in (2) 
but in its determinant [A ]. This is precisely 
what is necessary for the application of (6a, b). 
(See Section IV, (9) for example.) (c) In asym- 
metric molecules, where the coupling between 
spins does not vanish, [A] will contain the 
spin coordinates explicitly, and expressions of the 
form (7) are obtained only as a first approxima- 
tion, depending on the asymmetry of the mole- 
cule. A specific example of this is given in (11), 


(12). 


IV. APPLICATIONS 


(A) As a first illustration of the direct appli- 
cability of (5), (6a, b), the limiting rotational 
partition function and entropy of tetramethy]l- 
methane will be obtained and compared with the 
summation over the quantum levels. 

The model for this molecule consists of four 
freely spinning gyroscopes located at the corners 
of a rigid tetrahedral frame. Seven coordinates 
are used; three, the usual Eulerian angles, refer 
the frame to axes fixed in space, while the other 
four refer the spin of the tops relative to the 
frame. Details will be found in Lacoste’s paper.* 
The kinetic energy is given in form (2). Its 
determinant [A] is a symmetric 7th order 
determinant, which after rearrangement, is 
simply evaluated to be: JJ,‘ sin? 0 


where J, refers to the reduced moment considered in III, 
I, is the moment of inertia of a methyl top about the 
spin axis. 


Applying formulae (5), (6a, b) 
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ASTON 


_ akT Ml, 4)! 
f ‘ f f f- f f. sin 0d0dyd pda \dazda3da, 


and Sp=R In e7/2221/2719/2(RT)7/2( 193 T 4) /h?. (8) 


The summation over the quantum states: 


=>. , a : % y > a > efi, mg, mj, m2, m3, m4) 


1@ mp my m Mm, m3 mM 
has been carried out by Kassel,® whose result is: 
e7/2221/2719/2(RT)7/2( 7 37 ,4)3 Ro 
h? 4 
+terms negligible even at 100°K, (9) 


Sr’=R In 





where ¢ =0.0022 at 100°K and is inversely pro- 
portional to T. Neglecting Ro/4 which, even as 
low as 100°K produces an error of less than 0.005 
cal./deg. (9) becomes identical with (8). 

(B) As a second illustration, the method will 
be applied to a general class of molecules, for 
only Case V of which (see below) has a quantum 
mechanical solution been attempted. 

The model (Fig. 1) is a molecule consisting 





2K.E.=A(6 cos g:+y sin 6 sin ¢:)?+B(—6 sin 9,+y sin 6 cos ¢)? 





of two freely spinning asymmetric tops sub- 
ject to the restriction that a line from the center 
of gravity of one top through the center of gravity 
of the molecule be the common principal axis, 
and also spin axis of each top. Molecular ex- 
amples will be given below. The usual Eulerian 
angles refer three axes fixed in each top at the 
center of gravity of the molecule, O, toa reference 
frame fixed in space, origin at O. @ is the angle of 
inclination of Z,;Z2 to OZ»; y is the angle be- 
tween the fixed axis Xo and the nodal line ON; 
gi is the angle between ON and Xj; ¢» is the 
angle between ON and X2 where X,, Yi, Xo, Y> 
lie in the equatorial plane. Let A, B, C be the 
moments of inertia of top 1 about X,, Y; and Z;; 
and D, E, F, be the moments of inertia of top 2 
about Xe, Y2 and Zs. 

The kinetic energy of rotation is expressible in 
the form: 





+D(6 cos go+y sin @ sin g2)?+E(—6 sin g2+y sin 6 cos ¢2)? 


The determinant [A ] of (10) is readily found to be: 


Formulae (5), (6a, b) may now be applied, by using (11) to give the rotational partition function 


+C(¢i+wW cos 6)?+F(¢2—wW cos 6)”. (10) 





[A]=CF sin? ((-AB+DE+(AD+BE) sin? (¢:—¢2)+(AE+BD) cos? (yi—¢2)]. (11) 





and the rotational entropy. The integral required is: 


SOP SS. 


to the form:* 


* In order to bring (12) to the form (13a), it is required 
that A>B; E>D. Interchanging the magnitudes of A 
and B, or of D and E in (12) will interchange the mag- 


x sin ([-AB+DE+(AD+BE) sin? (g,— 






Integrating with respect to @ and y, and after a slight change of variables, (12) may be brought 





¢2) +(AE+BD) cos? (gi— ¢2) |id0dYdeidg2. (12) 





2a 
9=2nm { (1—k? sin? u)*du, (13a) 





nitudes of the coefficients of the sin? (gi—¢:) and 
cos? (¢ — ¢2) terms, but will not alter the value of the 
integral. 
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ENTROPY OF NONRIGID MOLECULES 


where uU= ¢i— $2, 
\=(CF)}(4r)(2akT)2/h', 


m=((A+D)(B+E))}, 


A>B; E>D, 


(A—B)(E—D)\ 
x-( ) <1. 
(A+D)(B+E) 


This is a well-known elliptic integral of the second kind, and the result may be put in the form: 


2°r*(kT)? 
2, =—— 


ht 


{(CF)(A+D)(B+E)} {1 a aoa acl 


BR? 1° 3 
| (13b) 


92 92 42 


where the series in [--- ] rapidly approaches 1 as AB or DE. 


CaseIl. A>B; E>D; C#¥F 

Example: Derivatives of diphenyl, where atoms 
have been substituted symmetrically in either 
one or both of the phenyl groups. The partition 
function is (13b) and the rotational entropy is 
R In e?Q. 


Case II. A= E; B=D; C=F 

Two freely spinning equivalent asymmetric 
tops. 

Example: Diphenyl. 
0 = [2°r®(kT)?/h* ]{(CF)(A+D)?}3 


K* 1*-3K* 
x{1-—- sel, 


22 22.42 
(A—D)? 
K*? =———- and_ Sr=R In e?Qzo. 
(A+D)? 
Case III]. A= B; E>D; C#F 
Top 1 symmetric. 
Example: Toluene 


3=[2%a*(kT)?/h*] {(CF)(A+D)(A+E)}}; 
K=0 and Sr=RIne’Qs. 
Case IV. A= B; D=E; C#F 
Two freely spinning nonequivalent symmetric 
tops. 


Example: All symmetrically substituted deriv- 
atives of ethane of the type: 


Lt 
y-6-6-x 
Y X 


Q;=[2%r5(kT)2/h*]{(CF)(A+D)?}}; 
K=0 and Sr=RiIne?Qy. 


(13c) 


(13d) 


(13e) 





Case V. A= B=D=E; C=F; A'=2A 
Two freely spinning equivalent symmetric 


tops. 
Example: Ethane. 


Q; = 2°r®(RT)?(A’C)/h*; 


K=0 and Sr=RiIne’Q;. (13f) 


SUMMARY 


(1) A general formula is derived for the limiting 
rotational partition function and entropy of non- 
rigid polyatomic molecules (apart from nuclear 
spin and symmetry number corrections and as- 
suming no potential). In order to apply the 
formula, it is unnecessary to obtain the Hamil- 
tonian form of the rotational energy, and to 
perform the integration over the momentum 
space. 

(2) A previous formula for the rotational en- 
tropy, suggested empirically, is derived as a 
special case of the general formula. 

(3) Illustrations of the use of the formulae 
derived are made to: 

(a) The spherically symmetric top with in- 
ternal gyroscopes, such as tetramethylmethane, 
and compared with the summation using quan- 
tum-mechanical energy levels. 

(b) The general case of a molecule consisting 
of two freely spinning asymmetric tops subject 
to the condition that a common principal axis 
through the center of gravity of one top and the 
center of gravity of the molecule be a spin axis 
for each top; rotational partition functions and 
entropies are obtained for five special cases of 
this model, and molecular examples given. 
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The Vibrations of Polyatomic Molecules in Excited Electronic States 


A. B. F. Duncan, Department of Chemistry, Brown University 
(Received April 29, 1935) 


The vibration frequencies found in the excited states of the molecules CO2, CSe, H2O, ClOz, 








SO2, Os, NH3, PCl3, CoH2, C2H4, CeHe and (CH3)2CO are considered. The frequencies appearing 
most strongly in all cases, and in most cases the only frequencies occurring, are totally sym- 
metrical. This is in accord with the extension of the Franck-Condon principle to such molecules 


by Herzberg and Teller. 








HE normal modes of vibration of a number 
of molecules are now well known. The 
numerical values of the frequencies, given by 
Raman and infrared spectra, have been assigned 
correctly to the proper modes of vibration in 
many cases. The frequencies which are present in 
the excited electronic states have not received 
so much attention, and the experimental data 
have been few until very recently. 

The selection rules for vibration-electronic 
transitions have been given in some detail by 
Herzberg and Teller.! It appears from their 
discussion that in absorption from the non- 
vibrating level of the ground state, only totally 
symmetrical vibrations (those whose coordinates 
do not change sign under any operations of the 
symmetry group) are strongly excited in the 
upper electronic state. The rules governing 
weaker transitions from this ground state as well 
as other lower states are also discussed. Herzberg 
and Teller considered the spectra of four mole- 
cules in relation to their rules and predictions. 
The spectrum of CH;I was found to be in com- 
plete accord with theory. The fluorescence 
spectrum of benzene was also in accord. The 
spectrum of ClO, was in partial agreement with 
the theory also.? The analysis of the SO: spec- 
trum was not in accord with theory.’ 

Since the publication of this important paper 
more data on other molecules have been pub- 
lished. The purpose of the present note is to 
show that totally symmetrical vibrations are 
indeed the ones occurring in excited states in all 


1G. Herzberg and E. Teller, Zeits. f. physik. Chemie 
B21, 410 (1933). 

2H. C. Urey and H. Johnston, Phys. Rev. 38, 2131 
(1931). The more recent analysis of Ku removes all dis- 
agreement with the theory. 

3A recent analysis of the spectrum of SO, by J. H. 
Clements is in complete accord with the ideas of Herzberg 
and Teller. 
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cases where spectra have been adequately 
studied and analyzed. The other vibrations are 
either rather weak or entirely absent. For the 
sake of brevity the results are summarized in 
Table I, which is largely self-explanatory. 
References to the normal state vibrations are 
omitted, if they can be found in Stuart’s Mole- 
kiilstruktur.4 The designation of these vibrations 
by subscripts is that used by Stuart rather than 
the original investigators, in order to avoid con- 
fusions. The totally symmetrical vibrations are 
starred. 


DISCUSSION 


In the examples given above the normal vibra- 
tions are known with some certainty. Three cases, 
however, call for comment. In benzene there are 
only two totally symmetrical frequencies, and 
there is little doubt that one of them has the 
value 991. The only observed frequency in the 
excited state, 920, is closer to this value than 
any other, totally symmetrical or otherwise. 
In ethylene the situation is less certain because 
there are a number of normal frequencies not 
differing greatly in value. It is difficult to decide 
to which normal frequencies those in the excited 
states correspond. The assignment of Price, 
however, is as reasonable as any. 

In acetone there seems little doubt that the 
excited state frequency 745 corresponds to 788, 
and 210 and 325 correspond to 376. 1190 was 
considered by Noyes, Duncan and Manning to be 
a modification of 1220, or some other C—H 
bending frequency in this region. One of which 
must be totally symmetrical. The possibility that 
this is the C—O frequency v2 (7) (1710 in normal 
state) must also be considered.”! 


4H. A. Stuart, Molekiilstruktur, Julius Springer, Berlin, 
1934, 
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VIBRATIONS OF POLYATOMIC MOLECULES 


TABLE I. 








Molecule Normal vibrations Excited state vibrations 





m* 1322 M1 
U3 2363 
V2 667.9 


COz 1330-11105 





1,* 655 YA 
v3 1523 
V2 396.8 


600+50° 


Excited state vibrations 


Molec ule Normal vibrations 





m* 257 
Ve 190 
v3* 510 V3 eens 
U4 480 V4 


v, 230-265" 


%-— 








u* (3600) 
ve* 1595.4 
U3 3756 


» 1300, 1600, (800)5 


1(s)* 1975 u(s) 1810, 176013 
v2(s)* 3370 v2(s) ————— 
v(a) 3277(3288) v/a) 

5(a) 729 5(a) 

5(s) 600 5(s) 300-0 








v,* 1105 
V* 529 
U3 954 





vu* 1146 
v2* 560 
U3 1340 


10992 
566.7 


v,* 1046 
P 700 
U3 1357 


1 878, 936, 902, 9541 





v,* 964. 
Ye 1630 
v3* 3336 
vq (3300) 


3, 933.8 





1% (as)* 3019 
v2(rs)* 1623 


v11(x%s) ————— 
v.(rs) 1340, 800, 
230, 200 
5(as)* 1342 5(ars) 480, 150 
(9 others not totally 
sy mmetrical) 


300,13 


M1 “900 tenes esc an 
920 (absorption)!* 


v,* 991 
v.* 3063 
(only totally sym- 
metrical vibra- 
tions) 
(CH,),CO". 19 oun) 788 
(r)* 1712 
a(n) 376 


v( ) 715° 20 

Vo( ) — 
v3( ) 210, 325 
? 1190 








Thus it appears that the symmetrical fre- 
quencies are most strongly excited by absorption 
from a nonvibrating ground state even in rela- 
tively complicated molecules. The same must be 
expected in more complicated molecules which 
give discrete spectra. So it may be concluded 
that a search for these frequencies should form 


°G. Rathenau, Zeits. f. Physik 87, 32 (1933). 

°G. Rodloff, Zeits. f. Physik 91, 511 (1934). 

7Z. W. Ku, Phys. Rev. 44, 376 (1933). 

8 J. H. Clements, Phys. Rev. 47, 224 (1935). 

°O. R. Wulf, Proc. Nat. Acad. Sci. 16, 508 (1930). 

10 A. B. F. Duncan, Phys. Rev. 47, 822 (1935). 

1K, W. F. Kohlrausch, Der Smekal-Raman- Effekt, Julius 
Springer, Berlin, 1931. 

BA, B. Duncan, unpublished observations of far 

ultraviolet spectrum. 

8 W. C. Price, Phys. Rev. 47, 444 (1935). 

“ Reference 4, p. 337. 

® E. B. Wilson, Jr., Phys. Rev. 45, 706 (1934). 


a logical starting point in the analysis of com- 
plicated spectra. 

Information on absorption from higher levels 
of the ground state is available in so few cases 
that it does not seem profitable to apply the rules 
of Herzberg and Teller which cover these cases. 


1D. E. Teets and D. H. Andrews, J. Chem. Phys. 3, 175 
(1935). 

039} B. Kistiakowsky and M. Nelles, Phys. Rev. 41, 595 
1932 

18 C. V. Shapiro, R. C. Gibbs and J. R. Johnson, Phys. 
Rev. 38, 1170 (1931). References to earlier work here. 

19 The C—C and C—O frequencies are those of HsCO 
(reference 4, p. 327). The nontotally symmetrical fre- 
quencies are probably 490, 526 and 1066. 

20W. A. Noyes, Jr., A. B. F. Duncan and W. M. 
Manning, J. Chem. Phys. 2, 717 (1934). 

*1 This was pointed out to the writer by Dr. W. C. Price. 
See also V. Henri, Comptes rendus 199, 849 (1934). 
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On the Stokes Phenomenon for the Differential Equations Which Arise in the Prob- 
lem of Inelastic Atomic Collisions 


O. K. Rick, Chemical Laboratory, Harvard University 
(Received March 7, 1935) 


The problem of inelastic atomic collisions, involving ex- 
change of energy between the atoms, depends upon the 
solution of a pair of coupled differential equations equiva- 
lent to a single differential equation of the fourth order, in 
which the independent variable is the distance, 7, between 
the two atoms. Asymptotig¢ forms for the solutions of this 
pair of equations may be found. The probability of energy 
exchange may be shown to depend upon the connections 
between the solutions for values of r which correspond to 
regions of space where the relative kinetic energy of the 


two atoms is positive and the solutions for values of r cor- 
responding to regions where the kinetic energy is negative. 
The finding of these connections involves a study of the 
asymptotic forms of the solutions when these forms are 
considered as functions of a complex variable. Such a 
study has been made in this paper, and the corresponding 
Stokes phenomenon investigated. The desired end was 
not achieved, but results of interest have nevertheless been 
obtained. A brief discussion is given of the previous work 
of Stueckelberg. 





$1. INTRODUCTION 


HE theoretical calculation of the probability 

of exchange of energy by two atoms on 
collision is a problem of considerable interest. 
It may be best treated by considering that the 
two atoms form an unstable molecule. The 
problem then resolves itself into the determina- 
tion of the probability of a radiationless tran- 
sition, whereby the molecule passes from one 
repulsive potential energy curve to another with- 
out change in total energy. 

It may be formulated mathematically as 
follows. We let ©; be the coniposite eigenfunction, 
representing the internal state of the system 
(that is, the electronic states of the two atoms and 
the rotational state of the unstable molecule 
which the two atoms are assumed to form), 
neglecting any interation between them. ®; obeys 
the wave equation 


(HH )— U;)$;=0, (1) 


where U; is the eigenvalue and [7 is the Hamil- 
tonian for the pair of atoms, neglecting any 
terms in the potential due to interaction of one 
atom with the other. H» depends upon the 
electronic coordinates and the angles giving 
the relative orientation of the two atoms, but the 
distance, 7, between the atoms does not enter 
Hy at all. U;, which gives the potential energy 
of the unstable molecule is, in this approximation, 
independent of 7. Only when the interaction 
between the two atoms is taken into account do 
we get the proper potential energy curves. 


We shall assume that only two states of the 
unstable molecule need to be taken into account. 
These two states we distinguish by letting 7 be 
replaced by 1 or 2, respectively. We shall further 
assume that a selection rule operates so that its 
rotational quantum numbers do not change in 
the transition. This would seem really not to 
be an additional assumption, but rather a 
necessary restriction required by the first as- 
sumption, since each set of rotational quantum 
numbers defines what is effectively a separate 
internal state of the system with its own po- 
tential-energy curve. If any other selection rule 
holds we should thus be compelled to consider 
transitions involving more than two internal 
states. 

The complete Hamiltonian for the system, 
which takes into account the mutual interaction 
of the two atoms, as well as their motion with 
respect to each other, may be written 


H=H)+V+H,, (2) 


where V consists of the previously neglected 
potential energy terms giving the mutual inter- 


action, and 

1/0 20 

H,= ——{ —+- ~), (3) 

x’ \dr? ror 
where x? =82zM/h?, where M is the reduced mass 
of the system and h is Planck’s constant. If we 
let ¥ be the eigenfunction corresponding to the 
Hamiltonian (2) we have 


(Ho + V+H,—E)y=0, (4) 
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ON THE STOKES 


where E is the energy of the system. Now, if, as 
we have assumed, only two states of the system 
need to be taken into account, then we may 
write 

y = F\®,/r+ F.2/r, (5) 


where F; and F, are functions of r. If we substi- 
tute (5) into (4), multiply through by #,; and 
integrate over all values of the internal co- 
ordinates, remembering (1) and the orthogonal 
properties of ©; and ®, and noting that //,(F;/r) 
=—K rd’ F;/dr’, we get 


U, F\ +01: F\+212Fe- xk a? F,/dr’ —EF, = 0, (6) 


where! 
Vui= fP, V®,dr (7) 


and Vi2= JP, VP, (8) 


where dr represents the volume element for all 
the coordinates except r and the integrals are 
taken over all allowable values of these co- 
ordinates. We now let V;=2::+U,. It will be 
seen that V,; is the quantity which gives the 
potential energy of the unstable molecule. By 
rearranging Eq. (6) slightly it becomes 


a? F\/dre+ r(E = V,)Fy = KV oF. (9a) 


Similarly 


d?F2/dr’+ «(E— V2) Fe= 012F), (9b) 


since V2:=Vi2. This pair of coupled differential 
equations, which are to be solved simultaneously, 
was originally obtained by London.? 

The most interesting case arises when the 
potential energy curves V; and V2 intersect, as 
shown in Fig. 1, and the point of intersection is 
of especial importance in the solution of the 
equations. We shall, in fact, find it convenient 
to set p=r—7, where 7 is the value of 7 at 
which the intersection occurs, and to let the 
arbitrary zero of energy coincide with the energy 
of the intersection. 

The pair of coupled differential equations, (9a) 
and (9b), is equivalent to a single differential 
equation of the fourth order. For large positive 
values of p, it may be shown that the four 
independent solutions represent, respectively, a 
stream of particles (really pairs of particles) 


‘We assume #, and 4 are real. If this is not the case 
the necessary modifications can readily be made. 
* London, Zeits. f. Physik 74, 143 (1932). 
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going out (toward larger p’s) on Vj, a stream 
coming in on V,, a stream going out on V2, anda 
stream coming in on V». Now these solutions 
connect with the solutions for negative values of 
p, two of which decrease exponentially as p goes 
to greater negative values, while two increase 
exponentially. The latter two solutions are ex- 
cluded as having no physical significance, which 
is equivalent to placing certain conditions on the 
coefficients of the solutions for positive values 
of p. This enables us to find the transition 
probability, as we shall see, provided we are 
able to connect the asymptotic solutions for p 
positive with those for p negative. Now Stueckel- 
berg* has attempted to find this connection by 
Stokes’ method of tracing the changes in the 
asymptotic solutions as they are carried over a 
large circle in the complex plane.‘ Stueckelberg’s 
work, however, was confined to the case in which 
the energy, E, was very different from the energy 
of intersection of the two curves. But the work- 
ings of the Franck-Condon principle, as well as 
the important role which the intersection of the 
two curves plays in Stueckelberg’s analysis, 
would lead us to expect that the most interesting 
case would be just that which Stueckelberg’s 
work excludes.® I have, therefore, attempted, by 
extending the method used by Stueckelberg, to 
solve the problem for the general case in which 
E is allowed to have any value whatsoever, 
specializing, however, the form of V; and V2, in 
fact setting 


Vi=—aip and V2=—azgp, (10) 


3 Stueckelberg, Helv. Phys. Acta 5, 369 (1932). For 
other approaches to the problem see Rice, Phys. Rev. 38, 
1943 (1931); Landau, Physik. Zeits. Sowjetunion 1, 46, 
88 (1932); Zener, Proc. Roy. Soc. A137, 696 (1932). 

4 Stokes, Mathematical and Physical Papers 4, 77, 283, 
5, 221, 283 (or Trans. Camb. Phil. Soc. 10, 105 (1857), 
11, 412 (1868), Proc. Camb. Phil. Soc. 6, 362 (1889); 
Acta Math. 26, 393 (1902)). 

5 Rice, Phys. Rev. 37, 1187, 1551 (1931); Jablonski, 
Zeits. f. Physik 70, 723 (1931). 
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where a; and a, are constants. It was believed 
that this specialization in the form of the 
potential energy curves would be a less im- 
portant restriction than that found necessary by 
Stueckelberg. It was found, however, that even 
with the restriction expressed by Eq. (10) it was 
not possible to solve the problem by means of 
the method used. In fact, during the course 
of the work, I have become convinced that 
Stueckelberg’s results, while they appear ex- 
tremely reasonable, are really open to grave 
objections. But in spite of the failure to attain 
the desired end and the generally negative 
character of the results, it has seemed to me 
that the considerations developed are of con- 
siderable interest mathematically, and should 
contribute to the understanding of this important 
problem, if not to its solution. 


§2. SOLUTION OF THE EQUATIONS 


The first step which must be taken is to find 
approximate forms of an asymptotic character 
for the solutions of Eqs. (9a) and (9b). Stueckel- 
berg has accomplished this by eliminating one 
of the F’s, let us say F2, and substituting a 
solution of the form 


Fy= el SothSith? Sot. -)/h (11) 
into the resulting fourth order equation, and 
equating separately to zero the coefficients of the 
various powers of h. As Stueckelberg gives the 
merest outline of this procedure, we shall show 
how it is done. There seems, however, to be 
some advantage to be gained by simply using 
the similar expression for Fo, 


Fo=elPothtith? Tat. -)/h (12) 
and substituting them together into Eqs. (9). 
Now if we are to set the coefficients of the various 
powers of h which occur on the two sides of 
each of the resulting equations equal to each 
other we must expand the exponentials. This at 
once leads us into difficulties on account of the 
negative power of h. We shall, therefore, sub- 
stitute the almost equivalent condition, that the 
solution given by Eqs. (11) and (12) must re- 
main valid for all values of h. By letting h go to 
zero we see that we must have Sp=7», a fact 
of some importance, as we shall see, which was 


RICE 


not brought out in Stueckelberg’s work. We now 
divide e%/*=e7/* out of the equations and 
expand the exponentials, letting primes indicate 
differentiation with respect to p. We obtain 
[(1/h?) (So +hSi! +h?So' +--+)? 

+(1/h) (So +hS1" +h?S2" +++ +) 

+ (82°M/h?)(E— V;) Je(1+h4S2+---) 


= (8m2M/h®)ve™(1+hT2+-++) (133) 


and a similar equation in which 7’s and S’s are 
interchanged and V2 occurs in place of Vj. 
This latter we shall denote as Eq. (132), and we 
shall use this subscript notation several times in 
the next few lines to designate equations related 
in this way in order to avoid having to write 
down both of them. The term in 1//? gives us 
(dividing (13) through by e*') 


So? +82?M(E— Vi) =82?Moye™-*, (144) 


Remembering now that 7)>=S» and letting V 
=(Vi+ V2)/2, the average of the two potentials, 
we see by adding (14;) and (142) 


So? /h?=T 9? /h? 
= x*[v12 cosh (T3—S1)—(E-—V)]. (15) 
Subtracting (142) from (14,) we get 
Vo— Vi= 212 sinh (T,—S)). (16) 


Using the relation between cosh and sinh we get 


from (15) and (16) 
So /h=To'/h 

= xl + {012° +3(V2—Vi)*}4—-(E—V)]}*. (17) 
Now, returning to Eq. (13,), the term in 1/h 
gives 
2S0'S1' +S" +.So"S2+822M (E— V1) Se 


= 89? M407 ve71- “1, (183) 


Using Eq. (14;) we find 
2S0' Si’ + So’ = 82? M42(T2—S2)e™- 1. (19) 


Adding (19:) and (192), remembering So= 7», 
and using (16) we get 


2So/ (Si +T1') +250" 


= 82?M(T2—S2)(V2— Vi). (20) 
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Subtracting (192) from (19;) we get 
So! (Si wy Ty’) 
= 82?Mvy2(T2—S2) cosh (T;—S;). (21) 


From (16) (Ve—V;)/2012=sinh (7,—5S;).. By 


differentiating both sides, 
[(V2—- V;)/2012] =(T'—S;') cosh (T,—S)). 


From this we can solve for S,'— 7)’. Substituting 
for S;/—T,' and cosh (7,—S,)(obtained from 
Eq. (16)) in Eq. (21) we get 


82?>M(T2—S2)= —2012S0'[(V2— Vi) /o12 |’ 
X (40.2?+(Ve—Vi)?]}-. (22) 

Substituting (22) into (20) and making some re- 

ductions we get 

(Sy' +71’) = —So""/So’ 

—3L(V2— Vi)?/40n? 'L(Ve— Vi)?/40n? +1 J. 


Integrating, 


(Si+ T)) = log So’ 


—log [(V2—V1)2/40122 +1} +C, (24) 


where C is the constant of integration. 
From (16) we have 


T;—S,=sinh™ { (V2—V1)/2012} 
=log {(V2—Vi1)/2012 


+((V2—Vi)?/4012? +1 ]'} =log w, (25) 
(where uw is defined by the equation) provided 
the sign of the inner square root of (17) is 


positive.® If the inner square root is negative’ 


T,—S,=—log ut+ir (26) 
where » has the same definition as before 
(positive sign in front of the square root). 

With Eqs. (17), (24), (25) and (26) we now 
have sufficient material to evaluate the ex- 
pressions (11) and (12) for Fy and Fo, respec- 


*See Peirce, A Short Table of Integrals, formula 679, 
which follows from the relation x =log (sinh x+cosh x). 

’ This inner square root is simply cosh (7;—S;), and if 
this is negative then 7;—S, is to be obtained by putting a 
minus sign before the square root in the logarithmic 
expression, which yields Eq. (26). 
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tively, as far as the second term in the expo- 
nential. Before doing this, it will repay us to 
establish a nomenclature which will distinguish 
the four different solutions and which we can 
use throughout the rest of the paper. In order to 
do this most effectively we shall proceed at 
once to use the approximation indicated in 
Eqs. (10), writing Vi=—ajp and V2=—azp. 
As in Fig. 1, we shall assume that V2>V; for p 
real and positive, hence we see that we must 
have 


(27) 


a) > Qo. 


Since we shall desire to consider the solutions, 
F, and Fs, as functions of a complex variable we 
shall write p= poe’, and using Eqs. (10) we may 
substitute into Eq. (17). We shall write 


1“ [— {v127-+4(a1— ar)" p97e?* } , 


—E—3(a;+a2) pee] (28) 


and 
vo=[+ {0122+} (a1— a2)? po7e?* } 3 


—E—}(ai+ a2) poe }'. (29) 


The reason for the particular choice of sub- 
scripts for the v’s is readily understood if we 
consider these expressions when pp is very great. 
Remembering (27) we readily find 


lim v1= ay! potei/2+ix/2 (30) 


po=a 


be id/2+in/2, (31) 


lim Vvo= as'po 


po=ao 


and 


In these expressions a;', ae’ and po! are to be 
taken as positive. It is thus seen that »;’ ap- 
proaches V; and 2? approaches V2, when 6-0, 
since E is negligible if po is large. 
Remembering now that Eqs. (25) and (26) 
are to be used with v2 or »;, respectively, we 
see that on substituting them, together with (17) 
and (24), into (11) and (12) we get four inde- 
pendent solutions of the differential equation. 
The four independent pairs of functions we 
shall write as F;,; and Fe, ; (¢=1, 2, 3, or 4). 
The expressions for the F;, ; are given below. 
(The Fs, ; are obtained from these merely by 
changing the sign of the exponent on u.) 
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F,, a= Agvy -tphe- "S40 = A gnu fy, 
F;, =A yy r~ther S140 = A nu fi, 
Fy, 3=Agve*h-t he" 240 = A gn gu} fo, 


F,, o= Aove ht teS "240 = A onou' fe 


(32) 
(33) 
(34) 
(35) 


where }\=[(V2—V1i)/2v:2 ? +1, where the A’s 
are arbitrary constants, and where the 7’s are 
abbreviations for the exponential parts of the 
solutions, and the f’s for the rest of the variable 
part, excluding the power of uw. The general 
solution may therefore be written 


F\=Ainu'fitAonu ‘fe 
+Aansu*fot+Asnu'fi, 
F.= Anim *fitAonou'fe 
+A snsu'feotA anu fi. 


§3. THE TRANSITION PROBABILITY 


It is readily seen that the limit of uf; as p 
becomes infinite is 24v,~}, that of p>f2 is 22v27}, 
while wf; and wf, become zero more rapidly 
than v;~? or vg}. Hence to a first approximation 
we may write for p real and positive 


= 24y,-¥(0)(A1(0) (0) +.44(0) n4(0)), 
lim F2(0) 


=D 3y—-1(0)(A.2(0)92(0) +A 3(0)93(0)). 


We introduce here a notation of which we shall 
hereafter make much use, whereby we insert in 
parentheses the value of 6. It is necessary to 
designate the value of @ even in the case of the 
A’s, for they are not the same for all values of @, 
this being, in fact, the essence of the Stokes 
phenomenon. 

Now A,(0)v:~*(0)9:(0) represents a stream 
of particles going from right to left on the 
curve V; with a velocity always proportional 
to |v,(0)| and A4(0)v:~-*(0)n4(0) represents a 
stream going from left to right on the same 
curve with the same velocity.’ In the same 


(37a) 


(37b) 


8It is really quite arbitrary which expression should 
represent which stream of particles, but we must have a 
given sign of the integral in the exponential correspond to 
a particular direction of motion in both (37a) and (37b). 
(By the term ‘‘stream of particles’”’ we really mean stream 
of pairs of particles. When we say the stream is ‘‘going to 
the right’’ we mean that the distance p is increasing.) 
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way A2(0)v2~3(0)92(0) and A3(0)v2~}(0)73(0) rep- 
resent streams of particles going respectively 
from right to left and from left to right 
on the other curve. The densities of these 
streams of particles will be found by multiplying 
these expressions by their conjugate com- 
plexes and will be, respectively, 2| A (0) |?»:(0)~', 
2|A4(0) {#»,(0)—, 2|A2(0) |2v2(0)—, and 2| A3(0) |? 
X v2(0)~!. Since »:(0) and »2(0) are proportional 
to the respective velocities we see that | A,(0) |’, 
|A4(0)\*, |A2(O) |? and | A;(0) |? are proportional 
to the currents (particles per unit time crossing 
a given point) in the various streams of particles. 
Now if we let A2(0) =0, then we see that all the 
incident particles coming from right to left are 
on the curve V;. After the collision they are 
traveling from left to right and a certain number 
have gone over to the curve V2. The proportion 
which go over, on collision, to the curve V2 is 
the transition probability ; we shall designate it 
as P (P, or P2 in the two cases considered below). 
We thus may write: 


If A.(0) =0, 
P,=|A3(0)|?/|A1(0) |’, 
| A (0) |?= | A3(0) |?+ |A4(0) “ 


(38) 


(39) 


and 


the last so that we may have conservation of the 
number of particles. On the other hand, we have: 


If A,(0) =0, 
P2=|A4(0) |?/|A2(0) |’, 
| A2(0)|?= |A3(0) |?+|A4(0) |?, 


in this case the transitions being from the curve 
V2 to the curve V,. 

Now it will be observed that of our solutions 
ni and nz become exponentially infinite and 73 
and 4, become exponentially zero as p goes to 
— . In order, therefore, that F; and F2 should 
satisfy the necessary conditions for a wave 
function we must have A;(7) and Ae(z) equal 
to zero. This is sufficient to determine the value 
of P, and the problem resolves itself into finding 
the relation between A,(0) and A;(z), etc. 


(40) 


and (41) 


§4. CHARACTER OF THE APPROXIMATE SOLUTIONS 


The approximate solutions which we have 
found have an asymptotic validity with respect 
to variation of p. This is roughly seen if we try 
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to find S2 and 7». If we take the coefficient of h® 
in Eq. (13:) neglecting S; and 7; where they 
occur and set it equal to zero, we get 


2So' So! + (So +2S0'S1')S2+S1°+5S," =0. (42) 


Now it will be observed that, as po goes to ~, 
the coefficient of S.’ changes as po’, the coefficient 
of S. changes as po~*, and the other term changes 
as po’. Therefore Se itself will change as po}. 
The same will be true of 7». Therefore, the 
approximate solution we have found, will, in the 
particular case we are considering, become a 
better and better approximation as po goes to ©. 
To make this rigorous, we should, of course, con- 
sider S3, 73, Sy, Ts, etc., but we may confidently 
accept the conclusion as correct. 

Now S2 and 7»: come into the expression for 
the solutions of the equation as exponentials, 
e*: or e”2, as the case may be, which multiply 
the rest of the expression. We see that this 
factor differs from 1 by an amount which goes 
to zero as 1/po? as po goes to ©, while @ remains 
constant, and the errors in a special solution, 
for example F;, ; will approach F;, 1 O(po~?). 


§5. THe SToKEs PHENOMENON 


The problem of the asymptotic representation 
of the solution of the uncoupled differential 
equation 


8G /dx2— «VG =0, (43) 


where V = —ax, has been considered in great 
detail by Stokes. The solution of this equation 
can be written in the form 


G= B,G,+ BG, 
= Byx-tetialz! 4 Boyletialesds (44) 


this expression has the same sort of asymptotic 
validity as the solutions (36) of Eqs. (9). How- 
ever, in this simple case, it is possible to obtain 
a complete asymptotic expansion for each of the 
solutions, which, of course, enables one to get a 
much better approximation than is given by (44). 
But even by using these asymptotic expansions, 
it is found that, in a certain region near those 
lines in the complex plane for which the ex- 
ponents of the expressions in Eq. (44) are real 
(the solid lines of Fig. 2), the error in the 
asymptotic representation of the large solution 
becomes greater than the entire magnitude of the 
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small solution. These regions become narrower 
and narrower as one goes out to |x| =, but 
along the lines themselves it is always true that 
the error of one solution becomes greater than 
the other solution. If now we cross this region, 
then, to use the words of Stokes, unless the 
coefficient of the large solution is zero the small 
solution becomes lost in a mist, and on emerging 
we will find that the arbitrary constant of the 
small solution will have changed. The change, as 
Stokes has shown, will be proportional to the 
coefficient of the large solution.? The constant of 
proportionality we shall call the Stokes coeffi- 
cient. 

Now if we do not use the asymptotic expansion 
but instead the much rougher solution (44) then 
the misty region becomes much larger. As |x| 
approaches infinity it is seen that one solution 
becomes infinitely large and the other infinitely 
small. The order of the errors is just the same as 
in the case of the coupled equations, discussed 
in §4. Thus the error of one solution becomes 

® Although implicitly contained in the work of Stokes, 
this proposition is not very clearly stated there. It may be 
easily seen as follows. If the coefficient of the large solution 
is zero, then the small solution does not get lost in a mist, 
and its coefficient remains unchanged as the critical region 
is crossed. Thus the small asymptotic solution alone 
represents one solution of the differential equation on both 
sides of the critical region. On the other hand, if the 
coefficient of the small asymptotic solution starts at zero 
and the coefficient of the large asymptotic solution is finite, 
then the coefficient of the small one changes by a definite 
amount as the critical region is crossed, and this gives us 
the representation of another solution of the differential 
equation on both sides of the critical region. A general 
solution is a linear combination of the solutions described, 
and it will be observed that the change in the coefficient 
of the small asymptotic solution in any such linear com- 


bination is proportional to the coefficient of the large 
solution. 
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infinitely great as |x| goes to infinity, compared 
to the entire other solution, over the whole 
region of the complex plane except certain 
regions near the dotted lines of Fig. 2, which are 
midway between the solid lines. As one goes 
from one of the dotted lines to the other the 
smaller solution will be lost in a mist nearly the 
whole way, but exactly the same argument used 
by Stokes? will suffice to show that it changes by 
an amount proportional to the coefficient of the 
large solution. We may formulate this as follows. 
Between 6=0 and @=27/3 the second solution 
is large. We therefore have 


B,(27/3) =B,(0) —aB,(0), (45) 
where —a is the Stokes coefficient. Similarly 

B2(42/3) = B2(27/3)+0B,(27/3), 
where 0 is the Stokes coefficient, and 


B,(27) = B,(47/3)+cB2(47/3), 


(46) 


(47) 


where c is the Stokes coefficient. 

Now if we carry the expression exp (3ia«x*) 
around the complex plane, letting 6 go from 0 to 
2r, it will change into exp (—3ia'xx*), and vice 


versa. And as 6 goes from 0 to 27 we see that x? 
changes to —ix~'. Therefore G,(27) = —iG2(0) 
and G2(27)=—iG,(0). But, as the differential 
equation has no singularities in the finite part 
of the complex plane, we must have G(27) =G(0). 
Hence B,(27) =71B2(0) and B2(27) =7B,(0). Suc- 
cessive use of (45), (46) and (47) together with 
the last two equations gives us!’ —a=b=c=1. 

It is thus possible to evaluate the Stokes 
coefficients for this simple problem. A similar 
method may be applied to the more complicated 
case of the coupled differential Eqs. (9). Before 
we can do this, however, we need to consider in 
more detail the nature and properties of the 
solutions of these equations. 


§6. PROPERTIES OF THE APPROXIMATE SOLUTIONS 
OF THE COUPLED DIFFERENTIAL EQUATIONS 


The properties of the approximate solutions of 
the coupled equations will be largely determined 
by the properties of v; and v2. To determine 


10 This method of treatment is due to Zwaan, Utrecht 
Diss. (1929), who considered it in connection with a dis- 
cussion of the Wentzel-Kramers-Brillouin approximation 
method. 
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these properties we examine Eqs. (28) and (29). 
The situation can be best visualized by con- 
sidering E+y,? and E+>»,* as functions of the 
real variable p(p=poe*®, the value of po being 
always positive, while @ is 0 if p is positive and 
a if p is negative). The straight lines of Fig. 3 
represent the original potential energy curves, 
while one of the dotted curves is E+y,;? and the 
other E+ v2”; these curves represent what may be 
called perturbed potential energy curves. Just 
which one of the dotted curves is to be associated 
with y; and which with v2 is indeterminate and 
depends, in fact, on p in a manner which may 
be arbitrarily fixed, except that, by (30) and 
(31), for very great positive values of p the upper 
dotted curve is E+», and the lower E+p,’, 
while if p has a very great negative value the 
reverse is true. At some intermediate point the 
identification changes. That this must be the case 
is obvious by inspection of (28) and (29). 
When p= +2%v12/(a1—az) the inner square roots 
of v; and v2 vanish. If we start from a point on 
the real axis, follow a path which encloses either 
one of these values of p, and return to the original 
point on the real axis »; will have changed into 
ve and vice versa. Now the two paths shown in 
Fig. 4 are obviously equivalent; therefore in 
going from A to B the identification of the dotted 
curves in Fig. 3 will have changed ; but the large 
semicircle is just the path implied in the case of 
Eqs. (30) and (31). However, if we take the path 
along the real axis, the point at which we depart 
and go around the branch point, and so the point 
at which the identification of the dotted curves 
of Fig. 3 changes, is quite arbitrary. This being 
the case, we shall find it convenient to let this 
point of change be far to the left of the point p; 
of Fig. 3. Then we may evaluate the constants of 
integration for the integrals fvidp and f vedp by 
writing them as definite integrals, /,,°vidp and 
J;.°vedp. In order to find the value of the first, 
let us say, of these integrals, at a point p whose 
po is large, we may break the integral in two 
parts, the first part being along the real axis 
from p; to po, the second part along the arc of a 
large circle with center at the origin and radius po. 
We may thus write 


p Po p 
f ridp= { rido+ [ vidp. 
Py Pa Po 


(48) 
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Energy avis 
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We may evaluate the second integral on the 
right of (48) if we use the binomial theorem first 
to expand the inner square root of Eq. (28) and 
then to expand the entire right-hand expression 
of (28). In this way we find 


E 
[rdo=iea! pf (4+ —— + +) de" 
* 60 2a; poe? 


= tas} po![ (2e*/?+ (E/arpo)e®!?+-++) 
— (3+(E/aipo)+---)] 
=ias[(3p!+(E/a1)p!+---) 
— (3 po!+ (E/a1) po! +---)]. 


Let us now set z=p+E/a,. Then in the limiting 
case where p is very large we see that 


(49) 


3E | 
= pi+- ee, 


2 a 


(50) 


Thus from (48), (49) and (50) 


Pp Po 
lim f ridp= { vd p — Zia $29! + 2ia,'s! 
Pg-0o py py 


=$§/2+2ia!z!, (51) 


where zo is the absolute value of z and 


Po 
(f ridp— Yias!a). (52) 
PL 


It is readily seen from (30) that the expression 
on the right-hand side of (52) approaches a 
finite limit. It is pure imaginary since »; is pure 
imaginary on the real axis. It would be zero were 
E+»; a straight line coinciding with V,. 
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L 








In the same way 


ep 
vol p= 6" / 


Pa 


lim 


Pyco 


2+ 2ias!yl, (53) 


where y= p+ E/ae and 


Po ; 
6’/2=lim (f rstp—Bias'y.). (54) 
Poco po 


It can be shown that 6’ = —6 in the following 
way. 6 is equal to £,v1dp where £ >, indicates a 
path starting at pi, going along the real axis to 
po, then going clear around the large circle to po 
again, and finally coming back to p. For if 7 is 
carried clear around the large circle from p» back 
to po it changes sign, by (30), that is to say 

vi(2e) = —»,(0); (55) 
hence, remembering that z=ze’* and that, when 
we have gone clear around the large circle, 
6=27, we see by (51) that 


f rdo=6/2+ Bias" tert v,(2m)dp 
Pl 


= §/2—2ia;}z 4 f v;(0)dp 


and from (52) we see that we may write 
F ovidp=s. (56) 
F 0,vedp=s'. (57) 


In Fig. 5 we show the paths indicated by fo, 
and £>»,. The path of £ >, is EM NH, while the 
path for £», is shown, rather distorted, by 


Similarly 





YPpreing, /(a,-a,) 


abcdefghijkl. Now the value of the integral taken 
about the closed path EM NHIJKLABCDE is 
zero, so the integral of v2 taken about JJ] J KLA- 
BCDE is —6’. But the value of the integral 1 
taken about abcdefghijki, that is 6, is easily 
seen to be equal to the integral of ve taken about 
HIJKLABCDE, or —é6’. This appears if we 
remember that carrying »; or v2 around one of 
the branch points on the imaginary axis changes 
the sign of the inner square root and so effectively 
converts one of these quantities into the other. 
If we start with v2 at E and »; at a, then the 
integrals along the segments marked with capital 
letters are equal to those along the segments 
marked with the corresponding small letters. 
Furthermore, the integral along ef cancels that 
along gh since p; is not a branch point after »; 
has been carried about the branch point on the 
imaginary axis. 

We shall find it necessary for our subsequent 
calculations to know how our different solutions 
change when we go around a large circle in the 
complex plane from a point where 6=0 back to 
the same point but with @=27. It is easy to 
see from (51), (55), (56), and the original 
definitions of y; and n, that 


m(2r)=ns(O)e® and 4(27)=m1(0)e, 
while, similarly, since 6’ = —6, 
no(2r)=ns(O)e; and 3(27)=n2(O)e*. (59) 
It is also easy to show that u}(27) = —yu3(0), and 


(58) 


RICE 


that fi(27) =7f:(0) and fo(27) =7f2(0). From these 
we get 
Fy, (27) = a te Fy 4(O), 


and Fe, (27) = —te" Fe, (0), 


(60) 
(61) 


etc. 


etc. 


§7. THE STOKES PHENOMENON FOR THE 
CouPLED EQUATIONS 


We now have sufficient information at hand to 
consider the Stokes phenomenon for the coupled 
differential equations. We consider again Fig. 2, 
letting now the origin of the figure, from which 
the various lines radiate, be the point p=0. 
It is readily seen from (51) and (53), remember- 
ing that 6 and 6’ are pure imaginary, that, in 
the limit of large po, the integrals J}, *vidp and 
J,,°vedp are pure imaginary along lines parallel 
to the dotted lines of Fig. 2. There will be certain 
small regions near the lines defined by the con- 
ditions 6=0, @=27/3, 0=47/3, ---, where the 
ratios of the four solutions of the coupled equa- 
tions will be finite, and where the errors in the 
asymptotic forms will be small. In these regions 
all four arbitrary coefficients, A1, A2, A3 and A,, 
are significant. In all other regions the error in 
the asymptotic form of any solution is larger 
than the entire value of all smaller solutions. 
Thus, unless the coefficient of the largest solution 
is zero, all other solutions become lost in a 
mist throughout nearly the whole region between 
two of the dotted lines. It can be shown by a 
simple extension of Stokes’ reasoning’ that on 
passing through the misty region each coefficient 
will have changed by an amount which can be 
separated into several parts, the number of 
parts being equal to the number of solutions 
whose magnitudes are greater than the solution 
in question, and each part being proportional to 
the coefficient of one of the solutions of greater 
magnitude. The order of size of the solutions is 
entirely determined by the relative sizes of 11, 
ne, n3 and 7, since the ratios of these quantities 
become exponentially infinite if one approaches 
po=® in any other direction than those indi- 
cated by the dotted lines in Fig. 2. Since 


ns>ns>ne>m if O0<@0<22/3 or if 44/3 <0<27 


and n> n2>s> 4 if 21/3<0<47/3 


we may write, assuming that we begin with 
6=27/3, 
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A i(0) =A (27/3) +a,A 2(2 1/3) +4 3A 3(2 1/3) +ay4A 4(2 1/3) 


A,(0)= 
A,(0)= 


A (42/3) =A (27/3) 


A,(27/3) +4d23A 3(2 1/3) +44A 4(27/3) 


A3(27/3) +434A 4(27/3) 
A,(27/3) 


A2(47/3) =A (27/3) +A o(27/3) 
A3(47/3) = b3,A 1(27/3)+b22A (27/3) +A 3(27/3) 
A,4(47/3) = byA (27/3) +b42A o(27/3) +by3A 3(2 1/3) +A 4(2 1/3) 


A (27) a A i(4 1/3)+c¢12A 2(47/3) +¢)3A 3(4 1/3) +¢4,A 4(47°/3) 


A,(27r)= 
A;(27)= 
A,(27)= 


where the a’s, b’s and c’s are the Stokes coeffi- 
cients. 

But since the coupled differential equations 
have no singularities in the finite part of the 
complex plane, we know, as in §5, that F,(27) 
=F,(0), which leads, in conjunction with (60) 


and (61) to the conditions 
A,(27r)=te—"A4(0), A3(27)=1e-“A2(0), 65) 
95 
A2(27)=ie"A3(0), Ag(27)=1e"A,(0). 


By use of (62), (63) and (64) it is possible to 
express all the A’s in (65) in terms of A,(27/3), 
A;(22/3), A3(27/3) and A4(27/3). Since the 
latter are completely arbitrary the multiplier of 
each one of them on the right-hand side of the 
resulting equations must be set equal to its 
multiplier on the other side. We thus obtain 
sixteen equations in the a’s, b’s and c’s, which 
are as follows: 


14+-Cy2be1+C13b31 + C14b41=0, 
CietCisdse+Cradae (66.2) 
Cis+Cisbas (66.3) 
C1 ie, (66.4) 
ba t-cobutcoda (66.5) 
1+ Cesbs2+Cosbae (66.6) 
Cos t+Coadas (66.7) 


(66.1) 


A (4 1/3) +603A 3(4 1/3) +¢94A 4(4 1/3) 


A;3(47/3) +63,A 4(42°/3) 
A 4(42°/3) 





(66.8) 

(66.9) 
(66.10) 
(66.11) 
(66.12) 
(66.13) 
(66.14) 
(66.15) 
(66.16) 


C24 =1e"ass, 
bsitCsaba1 =0, 
bs2+Csaba2 
1+ Csaba; 


= ie, 
=1e—"do3, 
=1e~" ae, 
bar =ie", 
bas =e" a0, 
bas = teas, 
1 =1e" 14. 
It turns out that of these 16 equations only 15 
are independent. They enable us to express the 
b’s and the c’s in terms of the a’s, and to get three 


relations between the a’s themselves. First of all 
we see that (66.16) gives us 


dys= —te—", (67.1) 


From (66.15), (66.12) and (66.11) we obtain 
(67.2) 


1—a 13294=> 1e~"a93. 


Substituting for c23 from (66.7) and for b32 from 
(66.10) in (66.6), using (66.12), (66.14), (66.8) 
and (66.15) to evaluate the rest of the quantities 
in terms of the a’s, and finally simplifying with 
the aid of (67.2) we get 


12024+13034 — 0120 34023= 0. (67.3) 
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It is obvious that further conditions will need 
to be found before we can evaluate the a’s and 
complete the solution of the Stokes problem. 
Before considering how other conditions may be 
discovered, however, we shall show how the 
transition probabilities may be evaluated in 
terms of the a’s. 


§8. EVALUATION OF THE TRANSITION PROBA- 
BILITIES IN TERMS OF THE STOKES 
COEFFICIENTS 


As we have noted at the end of §3 it is neces- 
sary to have A,(r) and A2(z) equal to zero. 
We have preferred in the later sections to con- 
sider the values of the A’s only for the directions 
of the dotted lines of Fig. 2; we can continue to 
do this, for it is readily seen that A,(7) and 
A2x(m) are zero if A,(27/3) and A2(27/3) are 
zero. This condition, it is seen, simplifies Eqs. 
(62) to a considerable extent. Now to calculate 
the transition probability P, we further set 
A,(0)=0, so that the second equation of (62) 
becomes 


0 = 193A 3(27/3) +44A 4(2 1/3) 


which enables us to express all the A(0)’s in 
terms of A3(27/3) and the a’s. We then get 
from (38) 


(68) 


| 1 — 34023/d24| , 
| dis — 014023/d24| . 





and (39) gives the equation 


| @13—@14023/d24| . 


= | 1 — 034023/d24| 24 | d23/d24| S, (70) 


Similarly from (40) and (41) 


|@13/a14|? 
P,= 


| 23 —13024/014| . 


(71) 





and 


| d23—o4013/A14| . 


= | 1 —13034/d14|?+ | @13/a14|*. (72) 


$9. OTHER RELATIONS BETWEEN THE STOKES 
COEFFICIENTS 


Other relations between the a’s may be ob- 
tained by considering the properties of the 
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solutions when only partial circuits are made. 
If we start on the positive real axis with any 
solution or combination of solutions which is 
real, then it is obvious that on carrying this 
solution to some particular @ the result will be 
the conjugate complex of that which would be 
obtained by carrying the same solution to —@ 
and the same po. Now suppose we set up the 
expressions , 
Fy(7)=Ai(m) Fi, 1(7)+A2() Fi, (7) 
+A3(7) Fi, 3(7) +A 4(2) Fi, 4(7) (73a) 
and 
F.(2)=Aj(2) Fo, 1(7) +Ao(7) Fo, o( 7) 
+A 3(7) Fs, 3(77) +A 4(2) Fe, a(7), (73b) 
where 


A;(r)=A,(37/2), 

Ao(m) = (b2;/2)A1(32/2)+A2(32/2), 

A3(m) = (b3:/2)A1(32/2) (74) 
+ (b32/2)A2(34/2)+A3(37/2), 

A 4(1) = (b41/2)A 1(34/2) + (ba2/2)A2(32/2) 
44(b43/2)A3(327/2) +A a(307/2). 


Now if and only if the expressions Fi(7) and 
F(a) are real (quite regardless of whether the 
A(z)’s for the small solutions have any sig- 
nificance when F)(7) and F2(7) are considered as 
approximate solutions of the differential equa- 
tions) we see that F\(47/3) = Fi*(22/3) (where 
the asterisk means conjugate complex) and 
F2(42/3) = F2*(27/3). Therefore, if Fi(0) and 
F,(0) are real it follows that F,(7) and F2(z) 
(as given by (73a) and (73b)) are also real. 

Now since y!, f; and fe start out by being real 
when @=0, and since 7:(0)=74*(0) and 72(0) 
=73*(0), it is readily seen that the conditions 
that F\(0) and F.(0) shall be real are 


A,*(0)=A,(0) and A2*(0)=A;(0). (75) 


To find the conditions which must be fulfilled 
in order for F(x) and F2(x) to be real we first 
note that, po being large and fixed, 


ulm) fun) =e" /4u3(0) fu(0) 
uA) f(a) =e *4u-4(0) fu(0) 
= —e-#4-¥(0) fx(0), 


(76a) 


and 
(76b) 





ON THE STOKES PHENOMENON 


while exactly the same relations hold with fe 
substituted for f;. From (51) the imaginary part 
of J5,°“ vidp is 6/2 and from (53) the imaginary 
part of J,,°“vedp is 6’/2= —6/2. We therefore 
see that in order for F;(7) and F;(7) to be real 
the following relations must hold. 


Ai(n)=e'7/4e-/q,, Ag (m) =e" !4e-*/°Q,, 


(77) 
Ao(r)=e'M4e!2q9, Ag(r) =e'*/4e8/2Q4, 


where 1, Q2, G3 and ay are real, the quantities 
which multiply these a’s being just what is 
necessary to counteract those factors of the 
particular solutions which are not real. 

Now, as we have said, if one of the sets of 
equations, (75) or (77), holds the other must 
hold also. This gives us further conditions on 
the a’s. We evaluate the A(0)’s in terms of the 
A(3x/2)’s from (62) and substitute into (75). 
We thus get a pair of equations we shall call 
(75A). Then we may evaluate the A(z)’s in 
terms of the A(3z/2)’s from (74) and substitute 
into (77) getting a set of four equations which 
we designate as (77A). Eliminating the A(37/2)’s 
from (75A) and (77A) we are finally left with 
two equations involving the a’s and also, of 
course, various combinations of the a’s and }’s. 
This final pair of equations we designate as 
(75B). Since the a’s are arbitrary, the quantity 
which multiplies any one of them on the left 
of one of the pair of equations (75B) must be 
equal to the corresponding multiplier on the 
right. We thus see that the two equations (75B) 
actually furnish eight relations between the a’s 
and b’s, which with the use of Eqs. (66) may be 
reduced to relations involving the a’s alone. 
Equating the multipliers of a, on both sides of 
one of the pair of equations (75B) gives Eq. 
(67.1) again; the same procedure on the other 
equation gives 
(78) 


a34>= — 1d24*e~**, 


Equating the multipliers of a; in the pair of 
equations (75B) gives 


(79) 


13 =13* 


and another relation which is not independent 
of Eqs. (67), Eq. (78) and Eq. (79). The coeffi- 
cients of a2 and q; give relations which are too 
complicated to handle directly. This being the 
case, further relations were sought, as shown in 
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the next paragraph, by invoking the law of 
conservation of particles, as expressed, in part, 
by Eqs. (70) and (72). This enables one to express 
five of the a’s in terms of the other one, and it 
may then readily be shown that the equations 
resulting from the coefficients of a2 and a; are 
satisfied. 

Multiplying (70) through by |ao,|? and (72) 
by |ai4|? and noting that, by (67.1) and (67.2) 
we have |d23414—24013|?=1, we see that 


(80) 
(81) 


12 | |= 
| d24— 34023| <4 | 223 | 2=1 


| {= 
| @44—@13034|?+ |Q13|?= - 


and 


Now by (67.1), (78) and (79), 


2 


! | ° ° 
| @14—@13034|*= | —1e-*° + d13024*7e" | 


ies * oa 
=1—413024 — 13425 +137d24 a24 


which, by (81) is equal to 1 —a;;”. It is thus seen 
that, since a;3 is not to be zero, we must have 


13 = (dog+G04*) (1+ ao4a24*)™. (82) 


It is now readily possible from (82), (67.2), (78) 
and (67.3) to get d23, @34 and dj. in terms of do3. 
We have 


do3= —1e"(1 —do4?)(1+deud24*) (83) 


and dj2= —A34= 124*e-", (84 ) 


whence 


P= P2= P= (do4t+Go4*)?(14+aeude4*)-*. (85) 

The expression on the right-hand side of (85) 
has all the characteristic properties of a proba- 
bility, inasmuch as it may be seen to be real and 
lie always between 0 and 1. Unfortunately, 
however, it has thus far proved impossible to 
find one more relation among the a’s which 
would permit the final evaluation of P. This 
failure, nevertheless, is in itself rather interesting, 
inasmuch as it shows what may happen when 
one attempts to apply Zwaan’s adaptation of 
Stokes’ method to more complicated differ- 
ential equations. In the case of the simple second 
order differential Eq. (43) the problem is solved 
by taking the asymptotic solutions once around 
the complex plane. However, in the case of 
the coupled differential equations, which are 
equivalent to a fourth order differential equation, 
even the more stringent conditions applied in 
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this section, which result from what may be 
roughly described as taking the solutions half- 
way around the plane, do not suffice to determine 
the Stokes coefficients. If we could find con- 
ditions to apply when the solutions are carried 
some smaller part of the way, say one-third of 
the way around, then we might hope to solve 
the problem. But the discovery and application 
of such conditions promises to be so complicated 
that nothing would be gained over solving the 
differential equations by series, a process which, 
theoretically at least, may always be carried out, 
but is actually forbiddingly involved. 


§10. CRITIQUE OF STUECKELBERG’S WORK 


It seems, at first, rather strange, in view of 
the above, that Stueckelberg should have been 
able to get definite results even for the special 
values of the total energy E to which he confined 
himself. It appears to me, however, that he 
was able to do this only by using a scheme which 
implicitly but quite arbitrarily introduces the 
extra condition which in the present investiga- 
tion it has not been possible to find. Stueckelberg 
assumes that the discontinuity in the coefficients 
of the solutions must always occur along lines 
where the ratio between two of the solutions is a 
maximum or minimum. For example, suppose we 
consider solution 1 and solution 2; then we find 
the locus in the complex plane where log (71/72) 
is real—such a locus will, at least in the limiting 
case, be the locus along which either 71/72 or 
n2/1, as the case may be, increases most rapidly 
as we go out in the complex plane away from 
the origin. Such a line naturally can be con- 
sidered to begin at a point where the logarithm 
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of the ratio of the 7’s in question is zero; the 
branch points considered above are such points. 
Stueckelberg’s analysis is based on the assump- 
tion that the ratio of two 7’s has acquired its 
asymptotic properties at a distance from its 
particular branch point which is small compared 
to the distance between branch points—this is 
what confines him to the case where E is far 
from the intersection of the potential energy 
curves. However, it seems to me that one is not 
allowed to make assumptions about the lines at 
which or the regions within which the jumps in 
the coefficients will take place, nor to decide 
whether a given solution causes the coefficients 
of one or more than one of the other solutions to 
jump, without a detailed investigation of the 
asymptotic series and the errors involved therein, 
and the limits of the regions in which these 
errors are large. In the simple case we con- 
sidered in §5 we showed that it was possible to 
do as much with the rough asymptotic forms as 
with the more exact asymptotic expansions, but 
this is by no means necessarily true for a more 
complicated case. Stueckelberg, furthermore, has 
not considered all possible pairs of solutions, but 
only those the Jogarithms of whose ratios become 
zero at one of the branch points. It is, again, 
not certain that this is justified. It is entirely 
possible that Stueckelberg’s analysis is correct; 
I do not believe, however, that it has been 
proved to be correct. It might not be impossible 
to find out whether it is correct, or if it is not 
exactly but nearly correct how good an approxi- 
mation it is under various circumstances, but 
this will require a detailed discussion of the 
asymptotic expansions which promises to be very 
tedious. 
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A possible error in Eyring’s recent calculations of absolute reaction rates due to the short life 
and consequent unsharp quantization of the activated complex is noted. The existence of this 
error is made more probable by a consideration of the target area required by Eyring’s equations 
at low temperatures. There is no doubt that his treatment becomes asymptotically correct at 


high temperatures. 





HE purpose of this paper is to comment on 

some aspects of the recent interesting at- 
tempt by Eyring* to compute absolute rates of 
chemical reaction from the geometry of the 
potential surface. Eyring considers that the 
activated complex has all the properties of a 
stable compound except for one degree of freedom 
which may be treated as a translation. He calcu- 
lates the partition function and hence the equilib- 
rium constant from the vibrational frequencies 
just as is done for stable molecules. It would 
seem that this procedure is open to question on 
the basis of the uncertainty principle, for with a 
velocity of 10° cm/sec. and a very generous 
reaction zone 5 X 10~* cm deep, the complex exists 
for only 5X10-™ sec. and has an uncertainty in 
energy of 2X10" sec.-! or 700 cm™. The vibra- 
tional quantization of the complex is thus very 
imperfect. 

If the quantum-mechanical transmission co- 
efficients could be calculated for a simplified 
potential surface the correctness of Eyring’s 
procedure could be decided definitely. The 
simplest surface which retains the essential 
features of the problem is 


E=0 for |x| >a, 


E=Eo+27°v’my? for |x| Sa. 


Here the barrier is a parabolic cylinder. Un- 
fortunately it has not been found possible to 
calculate the transmission of such a _ barrier. 
Some information can be gained, however, by 
studying the classical treatment of a similar 


1 Published by permission of the Director, U. S. Bureau 
of Mines. (Not subject to copyright.) 

* Associate physical chemist, U. S. Bureau of Mines, 
Pittsburgh Experiment Station, Pittsburgh, Pa. 

*H. Eyring, J. Chem. Phys. 3, 107 (1935). 


barrier. To make the problem correspond more 
closely to that solved by Eyring we suppose a 
reaction 


A+BC=ABC=AB+C 
with the potential energy surface 


E=27n?v97ud? for x>a, 


E=2rvpud?+Eo+27°v?mr*? for xSa. 


Here d is the separation between B and C, yp the 
reduced mass of BC and vy its vibration fre- 
quency; x is the distance between the center of 
gravity of BC and the projection of A on the BC 
axis; r?= y’+2", where y and z are perpendicular 
to the BC axis; m is the mass of A and »v the 
doubly degenerate bending frequency of the 
complex. For convenience we assume that uw is 
very much larger than m and that the moment of 
inertia of ABC equals that of BC. We are not con- 
cerned with the nature of the potential surface 
beyond the (degenerate) saddle point, since this 
will not affect the rate. For calculating the num- 
ber of collisions the molecule BC may now be 
replaced by the plane x=a. The number of 
collisions per cm? on this plane when there is one 
atom of A per cc is (kT /27m)*. The number of 
collisions on the annulus 7 to r+dr with a 
perpendicular component of velocity R to R+dR 
is 


dZ = (2am/kT)'e-" ®!2*? RrdRdr 


= (24/mkT )*e—®0/*7e—“!k Trdudr, 


where u=mR?/2—Ep. If every atom penetrates 
the potential barrier whose kinetic energy 
perpendicular to the barrier at least equals the 
height of the barrier where it strikes, the rate of 
reaction will be 
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(u/2x2v2m) $ 


b= f dZ 
u=0 “ r=0 


= -*(mm)HRT herb ue—"/*Tdy 


0 


= y~*(kRT/22m)ie—F/*?, 


For this same potential energy surface Eyring’s 
method gives 


ki=(h/kT)2(RT/22m)3e-20/ #7 (1 — eh kT)—2 


which reduces to the classically calculated value 
when hv<kT, but in general exceeds it by the 
factor (hv/kT)?(1—e—"”/*7)-. It is readily seen 
that if the classical reactivity of atoms with ex- 
cess kinetic energy wu is increased by a factor f(x), 
the condition necessary to give Eyring’s result is 


if uf(u)e—“!*du = (hv)?(1 —emhlk?)—2 
“0 


for all T. It has not been found possible to solve 
this equation for f(w), and it seems probable that 
no solution exists. Some light is cast upon this 
difficulty by a very simple calculation. We may 
certainly assume that no atom can react for 
which u<0, since then the energy is insufficient 
to penetrate the barrier at its lowest point. The 
number of collisions per cm? with u=0 is 
(RT /2mm)'e—*0!*?, From Eyring’s equation the 
number of successful collisions at very low 
temperatures is (h/kT)*?(kRT/2xm)%e—*0/*?, The 
ratio of these gives a target area of h?/2xmkT, 
which is of the order of 10~—* cm? and replaces o? 
of the usual theory, as Eyring has pointed out. 
Now for low temperatures the Maxwell distribu- 
tion curve is very steep, and the atoms with 
u=0 become a nearly homogeneous group with 
u=0, thus permitting easy calculation of the low 
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temperature target area for any assumed f(z). 
We need merely note that the target area from 
the generalized classical treatment is (2rmkT vy?) 
JS, uf(u)e-“!*7du, and that for small values of T 
the integrand is appreciable only near u=0. 
If uf(u) vanishes with u and behaves near u=0 
as u”, the target area for low temperatures will be 
proportional to 7”. If uf(u) is finite at 7=0, the 
target area becomes constant as 7° approaches 
zero. If uf(u) increases as u~™ when u goes to 
zero, the low temperature target area varies as 
T-"; the integral diverges, however unless 
m<1 and hence can never quite represent the 
Eyring target area, which varies as 7~'. It seems 
less important to discuss this merely annoying 
inability to find an f(u) exactly equivalent to 
Eyring’s result than it is to note that his equa- 
tion unquestionably requires an infinite target 
area for atoms with energy exactly Eo. Such a 
requirement seems physically unsatisfactory. 

It is recognized that arguments of this char- 
acter are not absolutely decisive, but they do 
at least suggest that the quantization of the 
activated complex is not sharp enough to justify 
the statistical methods used by Eyring, except 
when hv<kT; in that case classical theory gives 
the correct partition function, and unsharp 
quantization cannot harm the result. 


Note added in proof: Private correspondence and con- 
versation with Dr. Eyring after this paper was written has 
revealed that he agrees with the preceding discussion. His 
statement that when more accurate values for any part of 
the partition function were known from a source other 
than the Eyring-Polanyi energy surfaces, they could of 
course be used, was intended to include refinements of the 
type discussed here. There may be some question as to the 
wisdom of stretching the term partition function so far, 
and some disagreement in estimating partition functions 
in particular cases, but there is no disagreement with regard 
to the general program outlined by Eyring. 





JULY, 1935 


JOURNAL OF CHEMICAL PHYSICS 


VOLUME 3 


The Orientation of Unsymmetrical Molecules at Interfaces 


S. H. BAUER AND A. PoLLack, George Herbert Jones Laboratory, The University of Chicago 
(Received April 22, 1935) 


A function giving the distribution of the angles of orientation of the molecular axes from the 
normal to the interface can be obtained if one makes use of a number of special assumptions 
regarding the forces of orientation. This is compared with the distribution calculated on assum- 
ing an extreme form of Langmuir’s principle of independent surface action. The importance of 


having such a function is pointed out. 





INCE the publication of the often quoted 

paragraph by W. B. Hardy in 1912! on the 
orienting effects of unsymmetrical fields of 
molecules at interfaces, much work has been 
done in this direction which has led many to 
state that in cases of saturated films close to one 
hundred percent orientation exists. Such esti- 
mates were made on the basis of qualitative 
reasoning? as well as on quantitative calcula- 
tions. An obvious objection against the theory 
of surface orientation, sometimes presented as 
being of great importance while at other times 
minimized, is that the random motion of the 
molecules due to their thermal agitation will 
serve to overcome the orienting forces.‘ 

The purpose of this paper is to present a more 
or less quantitative estimate of this effect based 
on the classical principles of statistics and 
some of the present day concepts of oriented 
dipole layers, it being expressly understood 
that the picture used is a highly artificial one. 

From a kinetic theory point of view, one must 
consider a distribution in orientation so that an 
“oriented layer’’ really stands for a more or less 
defined layer of molecules whose molecular axes 
make an average angle with the surface which is 
not zero. In the case where the distribution in the 
inclination from the normal to the surface alone 
is to be considered, the generalized Boltzmann 


' Hardy, Proc. Roy. Soc. A86, 634 (1912). 

“ Harkins, Alexander’s Colloid Chemistry; Chem. Catalog 
Co., 1926, Vol. 1, p. 192. 

§’Langmuir, Colloid Symposium Monograph I11; Chem. 
Catalog Co., 1925, p. 48ff., who makes use of the principle 
of independent surface action, found the probability that 
a molecule like methyl alcohol should have its electrical 
axis oriented normally to the alcohol-vapor interface with 
its OH group in the alcohol phase is 300 to 10° times as 
great as the probability of orientation in the opposite 
direction. 

‘Compare: Langmuir, reference 3, p. 54 and Adam, The 
‘ea and Chemistry of Surfaces; Oxford, 1930, pp. 146, 
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function reduces to the well-known expression 
dn= Ce~V/*T sin 6d8, (1) 


where @ is the inclination of the molecular axis 
from the Z axis (taken normal to the surface); 
dn is the number of molecules between the 
angles @ and 6+d6; V(@) is the potential energy 
of orientation of the molecule at 6; and C is a 
constant, containing among other factors, the 
a priori probability of a given orientation. 

If it is assumed that C is independent of @ 
(which may not be true for interfaces where the 
molecules are crowded closely together so that 
the a priori probability will depend on geometric 
factors) it becomes obvious from (1) that 


62 
—1) e~Y@lKT sin 0d0 
M1 


61 7 

=| e-YIET sin Od0+ | e-Y/*T sin 6d0, (2) 
0 62 

where 1/s is the fraction of molecules oriented 

between 6, and 42. On the other hand, assuming 

random orientations, that fraction can be cal- 

culated from 


(3) 


In (2), 1/s can readily be found as soon as V(@) 
is known. By extending the method of Lang- 
muir®> this function can be obtained directly; 
whence the distribution function. Further, on 
differentiating V(@) one can get some idea of the 
type of force which one must imagine to exist 
near and at the interface in order to explain 
the orienting effects. This method will be il- 
Justrated in the latter part of the paper. 

One can, however, approximate V(@) by as- 
suming a field of force as a function of @ of a 


1/s=4(cos 6,;—Ccos 62). 


5 Reference 3, p. 50. 
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form which would be expected under these 
circumstances. If this were to result in an ex- 
pression for V(@) containing one arbitrary con- 
stant, its value can be estimated from the average 
potential energy of orientation per molecule 
which can be obtained experimentally; thus 


of e~V@/kTV(@) sin 040 
7 0 


j= : (4) 
of e~V@IKT sin 6d0 
0 





It is true that differences of opinion exist as to 
what experimental quantity should be used as a 
measure of 6 but fundamentally they are very 
much alike. 

Harkins considers, as will be done in this 
discussion, the energy of adhesion, defined as 


Oy Ove OY: 
ban (nt) +(n-12") (1-12) 
oT oT oT 


= E+ F2—E,;, (5) 





qT 


cy ce 
Ea qpotar) — Ea (nonpolar) “ee 


p= 
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where y represents the free surface energy and 
Ei, Es, E;, the total surface energies of the 
liquids 1, 2, and of the interface between them, 
respectively. E4, then, is the total energy neces- 
sary to tear apart a column one square centi- 
meter in cross section of liquid 1 from a similar 
column of liquid 2 when the ends of the columns 
are in contact. 

It can be readily seen from the available data 
that E, for organic liquids, particularly for those 
of aliphatic nature, differ very little from one 
another even when some are of the “polar” 
type, whereas such a difference is markedly 
shown, when the second liquid is water, in E;. 
This must be due to the fact that in cases where 
E; is small (i.e., when Ey, is large) the polar 
groups of the organic molecules in the interface 
are solvated in the water phase.’ Hence the 
difference between the E, for a hydrocarbon on 
water and that for the corresponding acid, 
aldehyde, alcohol, etc., on water is a measure of 
the potential energy of orientation per cm? of 
surface.’ We can now write 


Y@O/kT (rl sin 0+ xr? cos 8) sin 6d 0 





3 


where the second factor of the right member is 
the average area of a cylindrical molecule of 
length / and radius r projected on the interface 
while the first factor gives the potential energy 
of orientation per cm? of interface for the polar 
molecules. The factor 3 in the denominator is 
arbitrary but can be justified on the ground 
that (Ea(potar)y —EAqonpolar)) is the potential en- 
ergy per cm? of three processes of orientation, 
v1Z., 


(a) Change in the orientation of the water molecules in the 
interface as the bounding field changes from that due to the 
polar ends of organic molecules to that produced by water 
vapor; 

(b) Disorientation of the organic molecules in the inter- 
face—from the polar ends in the water phase to that of a 
random distribution; 

(c) Finally, as the water surface is sufficiently removed, 
there takes place another orientation in the organic liquid- 
vapor interface so that the nonpolar ends extend toward 
the vapor phase. Since we are interested in the energy 


’ (6) 


at 


cf eV MET sin 6d0 


en 





required for step (b) only, and further, since it is generally 
conceded that as much as 94 percent of the free surface 
energy is localized in the first layer of molecules,* the factor 
3 cannot be greatly in error. The form of the resulting 
distribution function will certainly be little affected. 


6 Harkins presents this view very clearly: reference 2, 
p. 208. “Ifa bar of . . . octyl alcohol of unit cross section 
is pulled in two at a plane, the molecules should first 
orient themselves on both sides of the plane where the 
break is about to occur. . . . The final work of separation 
will be effected most largely by the attraction between the 
nonpolar ends of the alcohol molecules. . . . Thus it 
would be expected that the work of separation of octyl 
alcohol from itself would be not greatly higher than the 
work of separation of octane from octane. . . . If, how- 
ever, the octyl alcohol is pulled away from water, the rela- 
tions are different. The hydroxyl groups should be turned 
toward the water. . . . The final break must come be- 
tween the polar OH groups of the alcohol and the water. 
Therefore the work of adhesion should be high.” 

7 For data which may be used in these calculations see 
Harkins and Ewing, J. Am. Chem. Soc. 42, 2539 (1920); 
Harkins and Cheng, ibid. 43, 35 (1921); Silberstein, Zeits. f. 
physik. Chemie A143, 157 (1929); Satterly and Colling- 
wood, Trans. Roy. Soc. Can. 25 (III), 205 (1931); Int. 
Crit. Tab. 4, p. 450ff. compiled by Harkins and Young. 

8 Edser, Fourth Colloid Report, 1922, p. 58. 
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ORIENTATION OF 





Langmuir employs in his calculations differ- 
ences in the E£; term only.® It is obvious, how- 
ever, that this method is a close approximation 
to that of Harkins since the F,’s differ very little 
for most organic substances. Adam,!° on the 
other hand, suggests that the latent heat term 
(—Tdy/dT) should not be taken into account 
in these considerations. 

At any rate, on equating expressions (4) and 
(6), one can solve for V(@) and for the unknown 
constant whatever value one uses for the average 
energy of orientation per unit area. To illustrate 
this procedure for obtaining the distribution 
function and to compare with similar results 
derived by extending Langmuir’s method as 
indicated above, we shall outline the calculations 
that were carried through for the case of octyl 
alcohol-water interface. The molecules are as- 
sumed to remain linear! and the force function 
is taken to be independent of 8, i.e., is a constant 
Fy) in a direction normal to the interface. Then 
obviously 





5a® a” 


MOLECULES AT 





INTERFACES 
V(0)= —M cos 8, (7) 


where M, depending on Fy and the length of the 
molecule, is the unknown constant to be de- 
termined by equating (4) and (6); V(7/2)=0 is 
here taken as the reference of potential energy. 
On substituting in the resulting equation cos @= x 
and M/kT=<a as is usually done, one gets 


akT f 


744 


1 —1 
edx= -pf e**(27rl(1—x?)} 
+1 
+arx)dx, (8) 


where p= (Eacp)—Eawp))/3 and the unknown 
to be solved for is now a. 

The right member of (8) can be integrated 
only when e” is expanded into a series. On 
performing the integration and simplifying, and 
at the same time also expanding the integral of 
the left member into a similar series for compu- 
tational purposes, there results the relationship 








1 a® a‘ 
prie(—-+— 4+ — 
2 16 384 9216 


na” 


an 





az a’ a‘ .a® 
=1(—4+—4— 4+ 
3 30 840 45360 


where » is an even number only, and m=n-—1. 
These series, although absolutely convergent for 
all values of a, do not converge very rapidly in 
this particular case, so that it was necessary to 
obtain the needed root of an eighth degree equa- 
tion. This was done by a method of successive 
approximations and we proved satisfactorily 
that further terms in the series will affect the 
result by less than 2 percent. 





*Values used by that writer are given in Alexander, 
reference 2, p. 536, and Colloid Symposium Monograph, 
reference 3, p. 53. See also work of Smyth, J. Am. Chem. 
Soc. 51, 2646, 2660 (1929). 

” His book, reference 4, p. 151. 

1 This is a point of departure from the ideas expressed 
by Langmuir (Colloid Symposium, 111) where he considers 
the molecules of long chained fatty acids and alcohols in 
such a position in the interface as to keep their polar ends 
anchored in the interface where they can behave as a 
two-dimensional gas while the rest of the chain twists 
and flips about in the body of the liquid. For octyl alcohol, 
this effect is probably not very important. 


3-5-7+++(n—1) ) 
n! 2-4-6-++(n—2)n(n+2) 


aa a’ a’ 
-) +part(<+—4 + + 
3 30 840 45360 
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The actual data used were:” 


E,.p) = 164.6 ergs/cm? for n-octyl alcohol at 20°C on 
water 
Eawp) = 107.4 ergs/cm? 
water 
ar? = 21.6 10-'* cm? 
2rrl =6.29 X 107'* cm? 
1/kT =2.49 x 10% 


for n-octane alcohol at 20°C on 


4, (10) 


wn 


giving a=2. 


from which the distribution function (1) was 
calculated. Curve 1 of Fig. 1 is a plot of the 
values of 1/s as calculated from (3) to be expected 
assuming random orientation, while curve 2 is 
the plot of the above calculated distribution 
(Eq. (2)) assuming a constant force acting normal 

12 The value of 7 was obtained from the area at zero 
compression (of film balance) for the aliphatic alcohols. 
The length was estimated from the molecular volume and 


value of r, as calculated from the density (=0.827) at 
yg Oe 
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to the interface and of magnitude consistent with 
the measured potential energies of orientation. 

We now extended Langmuir’s principle of 
independent surface action to an exaggerated 
degree by calculating the interfactial energy an 
octyl alcohol molecule has when it is oriented at 
different angles from the normal to the inter- 
face. Fig. 2 best illustrates the method used. 
The XY plane represents the water-alcohol 
interface; the cylinder represents the molecule 
oriented at an angle @ from the Z axis; the sec- 
tion of the cylinder above the plane CH FG is 
supposedly the polar end while that below the 
plane, the nonpolar end. It was then a simple 
geometrical problem to calculate the areas and 
their projections on the XY plane of the polar 
section in contact with the water and alcohol 
phases; similarly for the nonpolar section. The 
values of the areas at a given angle were then 
multiplied by their respective interfacial ener- 
gies and the products added to give the total 
energy of orientation at that angle. Again, 
V(/2) was taken as the zero of potential energy. 
From the resulting function V(@), curve 3 in 
Fig. 1 was obtained using Eq. (2). 





The geometrical constants needed in the cal- 
culations wer¢g rather arbitrarily estimated from 
the ratios of the lengths of the OH and CH, 
groups as obtained from. x-ray data™ and ou 
previous values of / and r. Thus 


r=2.62X10-* cm, lon =2.5X10-% cm, 
1=12.0X10-* cm, logHy, =9.5 X 1078 cm. 


z ¥ Er- 
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18 Clark, Applied X-Rays; McGraw-Hill Co., New York, 
1932, p. 319ff. 


ORIENTATION OF MOLECULES AT INTERFACES 


The four interfaces to be considered and their 
respective energies as estimated by Langmuir 
from such data as the heats of evaporation and 
the partial pressures of binary mixtures are:“ 


E; 
—12.6 ergs/cm? 
+ 6.7 
+24.4 
+37.4 


OH — water 
R—alcohol 
OH —alcohol 
R—water 


To say the least, the authors are convinced that 
little weight can be placed on such an analysis 
until many calculations of a similar nature have 
been carried through. It is to be noticed that 
although in cases 2 and 3 the V(@) functions 
differ markedly from one another and even more 
so do their corresponding force functions, the 
actual distributions in orientation are similar. 
Further, it is to be remembered that the con- 
clusions reached apply only to gaseous films 
where the a priori probability of any orientation 
is independent of @. As has been indicated above, 
in cases where the film approaches the condensed 
state when various crystallizing forces tend to 
bring in such a dependence, the constant C of (1) 
and of the subsequent equations would have to 
be expressed as an explicit function of the angle 
and made part of the integrand. The authors 
considered using a function which had maximal 
values at 6=0 and 7 to study the effect on 
the form of the distribution obtained. The 


“He does not give the values for Egy-alcohol and 
Er-alcohol. These were estimated from the energies 
Er-on, Er-r, and Egon whose magnitudes he gives 
as 33.7, 0 and —12.6, respectively, by taking into con- 
sideration the relative areas of the OH and R parts of the 
alcohol phase presented to the molecule in the interface. 
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arbitrariness introduced by selecting such a 
function at random is evident, particularly so 
when Rideal and his co-workers are inclined to 
believe that the maximal values really exist at 
other angles than 0 and z. The estimation of the 
ratio of the magnitudes of the maxima to those of 
the minima would be difficult if at all possible. 

We hesitate to discuss the significance of these 
particular distribution functions or of the ap- 
proximate validity of the method used to cal- 
culate them. We do wish to indicate, however, 
that did valid methods exist by means of which 
such distributions could be calculated, the result- 
ing functions would be used in all surface- 
chemistry discussions wherein ‘“‘orientation’’ is 
assumed or implied. Outstanding experiments 
which would be particularly affected are the 
contact potential measurements due to uni- 
molecular insoluble films spread on the surfaces 
of very dilute acid or alkali solutions. If the 
equation relating the potential difference due to 
the film with the average dipole moment of the 
molecules is approximately true, 


AV=2rnm 


the distribution functions obtained not only 
permit one to calculate the actual dipole mo- 
ment but also predict exactly the effect of tem- 
perature on AV (provided the nature of the 
film itself is not changed by the alteration of 
temperature). 

In conclusion, the authors wish to express the 
sincerest thanks to Professor W. D. Harkins for 
the numerous helpful discussions relating to the 


material of this paper. 


% Harkins and Fischer, J. Chem. Phys. 1, 852 (1933). 
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The Conductance of Potassium Chloride and of Hydrochloric-Deuterochloric Acid in 
H,O-D.O Mixtures. The Viscosity of H,O—D.O. 
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The conductances of 0.01 normal solutions of potassium 
chloride and of hydrochloric-deuterochloric acid in HxO-D.O0 
mixtures, and the viscosities of HxO-D.O mixtures, have 
been investigated as a function of D2O concentration 
at 25°C. The viscosity of D.O is 23 percent greater, the 
conductance of potassium chloride in D.O is 17 percent 
less, and the conductance of deuterochloric acid in D»O 
is 26 percent less, than the corresponding values in H.O. 
The viscosity, and the conductance of potassium chloride 


solutions and of hydrochloric-deuterochloric acid solutions 
show a maximum deviation of 0.4, 0.8 and 5.4 percent of 
the total values, respectively, from linearity. The decreased 
conductance of potassium chloride solutions in D.O is due 
mainly to the increased viscosity of the solvent. The 
marked deviation from linearity of the hydrochloric- 
deuterochloric acid curve is interpreted on the basis of a 
statistical correction to the proton exchange mechanism 
for the mobilities of H30*+ and D;0+ in H2O—-D.O mixtures. 





N. LEWIS and his co-workers,': ? using a 
e fraction of a cc of heavy water (90-97 
percent D.O, depending on the experiment), 
have shown that the mobilities of ions in D,O 
are less than, and that the viscosity of D,O is 
greater than, the corresponding values for HO. 
Because of the small quantities employed, the 
accuracy is open to question, and a more exact 
study of the problem is highly desirable. Further- 
more their use of linear extrapolation to obtain 
the values for pure D.O from the experimental 
values in H,O—-D,O mixtures is not always a 
correct procedure as will be emphasized below. 
The conductances of 0.01 normal solutions of 
potassium chloride and of hydrochloric-deutero- 
chloric acid in H,O-D.O mixtures, and the 
viscosities of HxO-D.O mixtures, were investi- 
gated as a function of D.O concentration at a 
temperature of 25°C. 


MATERIALS 


Potassium chloride was twice recrystallized 
from water and dried at 600°C. The constant 
boiling hydrochloric acid was prepared according 
to the usual method.* After atmospheric distilla- 
tions from an alkaline permanganate solution 
and from an acid chromate solution, the water 
samples were distilled in vacuum. The water in 
the vacuum still was frozen, the apparatus 
evacuated, the ice melted completely to allow 
the dissolved gases to escape, and the freezing, 
evacuation and melting repeated. The distillation 


1 Lewis and Doody, J. Am. Chem. Soc. 55, 3504 (1933). 

2 Lewis and MacDonald, J. Am. Chem. Soc. 55, 4730 
(1933). 

* Bonner and Titus, J. Am. Chem. Soc. 52, 633 (1930). 


was conducted at 40°C with the receiving flask 
at 0°C. The waters averaged about 1 10~-° mho 
in specific conductivity. The D.O concentration 
was determined from the increase in specific 
gravity, AS, by using the corrected Lewis 


formula, Np,o =9.377AS — 1.01AS*.4 


APPARATUS 


A vacuum tube oscillator of the beat frequency 
type® constructed in this laboratory® produced 
a current free from harmonics, and from any 
direct current component. The frequency was 
1740 cycles per second. The Leeds and Northrup 
conductivity bridge is described in detail else- 
where.’ An L. and N. impedance coupled ampli- 
fier amplified the bridge output sufficiently that 
a cell resistance of about a thousand ohms could 
be determined to within +0.001 percent. 

A small cell, shown in Fig. 1, having a capacity 
of 4.5 cc, was constructed of Jena 161! glass. 
The cell constant was determined according to 
the directions of Parker and Parker.* As a 
change of the frequency from one thousand to 
five thousand cycles produced a decrease in the 
measured resistance of only 0.01 percent, the cell 
was sensibly free from the Parker effect. 

The ordinary Ostwald U-shaped type of vis- 
cometer was employed. The time of outflow at 
25.00°C was 75 sec. for five cc ‘of H,O. 

4 La Mer and Baker, J. Am. Chem. Soc. 56, 2641 (1934). 

5 Lampkin, Radio Broadcast, July, 1928, p. 156. 

‘The authors wish to express their thanks to Mr. 
Samuel Korman under whose direction the oscillator was 
constructed. 


7 Paul H. Dike, Rev. Sci. Inst. 2, 379 (1931). 
8 Parker and Parker, J. Am. Chem. Soc. 46, 312 (1924). 
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ELECTROLYTIC CONDUCTANCES 


EXPERIMENTAL PROCEDURE 


Preparations of solutions 


The limiting factor in the precision of the 
measurements is the accuracy with which the 
concentration of the solution is known. In 
the preparation of potassium chloride solutions, 
between 20 and 35 mg of potassium chloride was 
transferred from a weighing bottle to a 50 cc 
flask. If the weights of the potassium chloride, 
determined by the difference in the weight of 
both the flask and the weighing bottle, did not 
check within 0.06 mg, the sample was discarded.° 
The average of the two weights was taken as the 
weight of potassium chloride used. Enough water 
was now added to make an approximately 0.01 
normal solution. 

The method of preparing the hydrochloric 
acid solutions differed from that employed for 
potassium chloride solutions, since the principal 
error arises from the loss by evaporation during 
transfer. The water for the solution (specific 
conductivity about 1x10-* mho) was placed in 
a 25 or 50 cc volumetric flask and about 45 mg 
of hydrochloric acid of constant boiling point 
added from a small weight burette, the amount 
of acid added being checked by the increase in 
the weight of the solution in the flask. In two- 
thirds of the cases the results agreed to within 
0.1 mg, and the maximum variation was 0.2 mg. 
As this loss in weight is more likely to be water, 
the only opportunity for evaporation of the 
hydrochloric acid solution being from the tip of 
the burette while the water has an opening about 
one (cm)? in area from which to evaporate, the 
decrease in the weight of the burette was taken 
as the weight added. 


Determination of conductance 


The conductivity cell was filled with the 
solution until the electrodes were in equilibrium 
with the solution, after which the cell was never 
washed with pure water. When ready for use, 
it was drained, dried, filled with the solution 
being tested and placed in the constant tempera- 
ture bath until its conductance became constant 
which required about an hour. It was now rinsed 
three times and refilled. After the first filling 
the difference in conductance between successive 


* The weighings were reproducible to +0.03 mg. 
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Fic. 1. Conductivity cell. 


fillings was rarely more than 0.02 percent and 
usually about 0.005 percent. 


Viscosity determinations 


After the viscometer, containing five cc of 
water, had been in the 25°C bath for ten minutes, 
the time of outflow was measured six or eight 
times and the results averaged. The water was 
removed, the viscometer dried, and the experi- 
ment repeated. Chromic acid “‘cleaning”’ solution 
was used but once as it has a solvent effect on 
the glass, the apparatus being merely air-dried 
between runs. At the conclusion of the experi- 
ments, the HO value was rechecked, and the 
time of outflow for H,O at 40° determined in 
order to evaluate the constants of the vis- 
cometer, C and C’, in the equation for relative 
viscosity (n= p(Ct—C’/t)).!% 1 
Calculation of results 


The calculation of the conductance is evident 
from Table I, where the data of a typical run 
are given for both potassium chloride solutions 
and hydrochloric acid solutions. The data for 
the potassium chloride solutions were corrected 
for the conductance of the solvent. This cor- 
rection is of no importance in HCl solutions, 
where the conductance of traces of H2CQ; is 
completely repressed.!” 


Experimental data 


The conductances of 0.01 normal potassium 
chloride solutions in HyO—D,O mixtures is given 
in the second column of Table II and curve I of 


10 Hatschek, The Viscosity of Liquids, D. Van Nostrand 
Co., 1928, p. 41. 

1 Reilly and Rae, Physico-Chemical Measurements, 2nd 
ed., D. Van Nostrand Co., 1932, Chap. 37. 

12 Davies, The Conductivity of Solutions, 2nd ed., John 
Wiley and Son, 1933, Chap. IV. 
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TABLE I. Typical conductance data. 








KCI Solutions 
10° x 


Sp.gr. KCl Sol’n Sp.gr. Eq./cc R _ 4.902 108 10° «x 
Solvent (g) (g) Sol’n (<10°) (ohms) R-10-* water KCl Np.o 


1.02202 0.02740 37.720 1.02249 0.9933 3640.0 1346.7 2.1 1344.6 
0.02263 32.029 1.02247 0.9661 3724.0 1316.3 1 1314.2 0.2060 





2. 
02678 34.532 1.02252 1.0605 3406.6 1439.0 2.1 1436.9 





HCI Solutions 
20.25% 


HCl A=4.902/R-C 
1.05923 0.04458 1.05930 0.9907 1484.8 333.2 

.04906 1.05931 1.0884 1351.8 333.2 

04368 1.05930 0.9693 1514.6 333.9 








TABLE II. Summary of results. 
All measurements at 25°C; electrolyte concentrations, 0.01 normal; values in parentheses are extrapolated. 








Relative 
Np:o Axcl Acer Acacl+peb \cH30++D30+)* viscosity Axei(n)? 8 


0.000 141.4 ta4 411.9 339.8 1.000 141.4 
206 135.7 69.2 378.1 308.9 1.046 
368 131.6 67.1 354.9 287.8 1.082 
.550 127.0 64.8 333.5 268.7 1.123 
.709 123.5 63.0 318.7 255.7 1.161 
.848 120.3 61.3 309.2 247.9 1.195 
-.980 117.3 59.8 303.0 243.2 1.228 
1.000 (117.0) (59.7) (302.1) (242.4) (1.232) 











1Tc]- in KCl equals 0.5098 (MacInnes, Shedlovsky and Longsworth, J. Am. Chem. Soc. 54, 2758 (1932). 
*\(HsO++D;0+) =A(HCI+DCl) ~ACI-- 


‘ 


Fig. 2. The value of Ao.o. in H,O, 141.4 Kohl- The results for 0.01 N HCI-DCI solutions are 
rausch units, agrees well with Shedlovsky’s given in Table II and curve II of Fig. 2. The 
value of 141.32. The value in 100 percent D,O, 
extrapolated from 98.34 percent D,O, is 117.0, a 
decrease of 17.2 percent in the conductance 
compared to the value in H,O. Lewis and 
Doody,' using an 0.086 molar potassium chloride 
solution, obtained a decrease in conductance of 
about 14 percent. The change in conductance 
with D,O concentration is not quite linear. The 
maximum deviation (50 percent D,O) is 0.8 
percent. 

As the ratio of the conductances of K* and 
Cl- is very nearly independent of viscosity,'®: 4 
the conductance of Cl- may be calculated if 
the assumption is made that the relative degrees 
of hydration remain the same. The values of 
Aci- thus calculated, with 0.510 as the trans- 
ference number of Cl- ion in KCl, are given 
in the third column of Table II. 

13 
» Re Tie I Soe ow gee oe Nae ae Solu- a Ate a a a 
tions, The Chemical Catalog Co., 1922, Chap. 5. ¥j 


46 MacInnes, Shedlovsky and Longsworth, J.Am. Chem. Fic. 2. Curve I, Axci. Curve II, Acaci+pcn- Curve III, 
Soc. 54, 2758 (1932). A\(H30++D;07*)- Curve IV, (Axci°7)- 
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ELECTROLYTIC CONDUCTANCES 


value of Ao.o: for HCl in HO, 411.9 Kohlrausch 
units agrees very well with Shedlovsky’s value 
of 411.88. The value of Ao.o: for DCI in D.O 
extrapolated from the three values between 
97.50 percent DO and 98.75 percent D,O is 
302.1, a decrease of 26.6 percent from the con- 
ductivity of HCl in H,O. Lewis and Doody,' 
using 0.017 M HCl obtained a decrease of about 
28 percent. The conductance at 50 percent D,O 
is less by eighteen Kohlrausch units (5.4 percent 
of the total conductivity at that D,O concentra- 
tion) than a linear relation would give. This de- 
parture from linearity is shown in Table III and 


TABLE III. Deviation of conductance of (HCI+DCl) from 
linearity. 








Alinear— Aexp. x 102 


Alinear 
_— Aexp. 


11.2 
16.6 


Avexp.) 


378.1 
354.9 
333.5 17.9 
318.7 15.3 
309.2 9.5 


Adinear) 


389.3 
371.5 
351.4 
334.0 
318.7 


Np, Oo 


0.2060 
3682 
511 
7095 
8488 


4texp. 











Fig. 3 where 
Alinear — Aexp. 
(Atinear — Acxp.) and a 4 10? 


Aexp. 


are given as functions of the D.O concentration. 

The viscosity of HyO—D.O mixtures is given 
in column 6 of Table II. The maximum deviation 
from linearity is 0.4 percent of the total at 50 
percent D,O. The relative viscosity extrapolates 
to 1.232 for pure D,O, and the absolute viscosity 
to 1.103 centipoises,!® in exact agreement with 
Lewis and MacDonald’s value.’ 


POTASSIUM CHLORIDE SOLUTIONS 


Numerous studies of the relationships between 
conductance of ions in aqueous solution and the 
viscosity of the medium": '4 have shown that 
the experimental results may be expressed satis- 
factorily by the equation, A-(7)? equals a con- 
stant where is usually less than one and, in 
general, approaches unity as the size of the ion 
increases. 

To investigate the influence of viscosity, the 
addition of a second solvent or a change in 
temperature has been employed. Unfortunately 
these procedures introduce a variation in other 


* Viscosity of HzO at 25°C =0.895 centipoise, I.C.T. 
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Deviation of Curve I 
2 S$ 8 & & 
« % 9 


a 9 
yn» & 
Deviation of Curve 


% 
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°o 10 20 WO AO 50 60 70 BO 30 100 
No.0 


Fic. 3. Deviation of conductance of (HCI+DCl) from 
linearity. Curve I, (Alinear— Aexp.)- Curve i, C(Atinear— 
Mexp.)/Aexp. ]X 10°. 


important properties of the medium, so that the 
effects cannot be ascribed exclusively to vis- 
cosity. By varying the proportions of H;O and 
D,O, the viscosity of the medium can be changed 
without producing these complications ; the tem- 
perature, the electrolyte concentration, the size 
of the solvent molecules, and even the dielectric 
constant!’ remains practically unchanged. 

The product (Axci:()°-*!*) gives a constant 
value within the limits of experimental error as 
shown in the last column of Table II. 7 is the 
relative viscosity of the pure solvent, which may 
be used since the viscosity of 0.01 N KCI solu- 
tions in H,O is only 0.04 percent greater than 
pure H,O.'8 Even on the assumption of a linear 
dependence upon the viscosity 90 percent of the 
decrease in conductance is accounted for by the 
increase in viscosity of the D.O. 


HYDROCHLORIC-DEUTEROCHLORIC ACID 
SOLUTIONS 


It has been recognized for many years that 
the conducting ions in acid solutions are not 
free protons, but hydrated protons (H,O*+)— 
oxonium ions.'* 2° When an H;O+t ion is in 
sufficiently close contact with a water molecule, 
as in aqueous solutions, it need not retain its 
extra proton since there is another configuration 
of equal energy in which the extra proton has 
changed molecules, thus: 

17 Lewis, Olson and Maroney, J. Am. Chem. Soc. 55, 4731 
(1933) find D for D2O is 99 percent of that of HO. 

18 Jones and Talley, J. Am. Chem. Soc. 55, 624 (1933). 

19 Bernal and Fowler, J. Chem. Phys. 1, 515 (1933). 


20 Eucken, Jette and La Mer, Fundamentals of Physical 
Chemistry, McGraw-Hill, 1925, p. 445. 
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H.0+H;0t2H,0++H.0. (1) 


At any given time there is an equal probability 
of the proton being on either of the water 
molecules. However, if there is an applied field, 
F, the probability of the proton being found on 
one or the other of the water molecules is altered 
by an amount proportional to F, the migration 
of the proton being directed toward the region 
of lower potential energy. 

The ions Cst+, Rbt, Tht, K+, NH,*, Cl-, Br- 
and I~ have about the same mobility in aqueous 
solutions, which is the maximum mobility of any 
ions except H;O+ and OH-. The ordinary 
mobility of ions of the same charge is primarily 
a function of their effective ionic volume and 
not of their mass. The conductance of H;O+ is 
about five times that of the ions just mentioned. 
If the mechanism of conduction of H;O+ were 
the same as that of other ions, the radius, calcu- 
lated on the basis of Stokes’ law, would be 
absurdly small (2.610-® cm). According to 
the present view, the conductance of H;O* must 
be considered as composed of two parts, (1) 
ordinary conduction similar to that shown by 
other ions and (2) proton exchanges leading to a 
Grotthus chain type of conduction. In Hy,O—-D.O 
solutions, in addition to proton exchanges be- 
tween H:O molecules, several other types of 
proton and deuteron exchanges are possible. 


D,0+D;0+2D;0++D,0 (2) 
HDO+H,DO+@H,DO++HDO (3) 
HDO+HD.0+tZHD.0t+HDO —_ (4) 
HDO+H,DO+@HD,0t+H.0 (5) 
HDO+HD,.0t@2H,DO+t+D.0 (6) 

H.0+D,0+2H:DO++D,.0 (7) 

D,0+H,0+@HD,0++H,0 (8) 


The exchanges represented by Eqs. (2), (3) 
and (4) are similar to that represented by (1) in 
having the same energy after the proton or 
deuteron transfer as before, while if the ex- 
changes represented by (5), (6), (7), or (8) occur, 
an absorption of energy from, or an evolution 
of energy to, the surrounding medium must 
occur. The frequency of the exchanges illustrated 
by (1), (2), (3) and (4) are all of the same order 
of magnitude, while the frequency of those 
illustrated by (5), (6), (7) and (8) will be less 


BAKER AND V. K. 


LA MER 


because of the necessity of exchanging energy 
simultaneously with the surrounding medium. 
At a concentration of 50 percent D.O, we have 
the maximum probability that an acid ion will 
not be next to a water molecule with which its 
extra proton or deuteron can exchange readily. 
We should accordingly expect a lower con- 
ductance than that calculated according to the 
additivity law. The conductance of (H;0+ 
+D,;0*) ; 1.€., of Aqici+pcp —Aci-> is given in 
the fifth column of Table II, and in curve III 
of Fig. 2, on the assumption that the conductance 
of Cl- is the same in 0.01 normal solutions of 
both KCI and HCI-DCI."® The equivalent con- 
ductance of H;0+ in HO is 339.8 Kohlrausch 
units; that for D,;0+ in D,.O is 242.4, a decrease 
in conductance of 28.6 percent. The equivalent 
conductance of (H;0++D,;0*) in 50 percent 
D.O is 17.6 units less (6.5 percent of the total 
conductance) than linear extrapolation would 
give.2! This marked deviation from linearity is 
exhibited by the empirical equation, \(4,0++p,0*) 
= 339.8 (H,O)+255.8 (HDO)+242.4 (D,0), 
which fits the experimental points with an 
average deviation of 0.2 unit. The conductance 
constant for HDO, instead of being the average 
of the constants for HO and D,O, is only slightly 
larger than that for DO. 

It is important to point out that the deviation 
of the conductance curve for (H;0++D,;07*) 
from linearity discussed above is not due to the 
effect considered by Onsager and Fuoss.”* They 
show that the conductance of a solution of two 
electrolytes with a common ion but with non- 
common ions of widely differing mobilities should 
not follow the additivity law. For NaCl and 
HCl at a concentration of 0.005 molar each 
(total concentration 0.01 molar), the deviation 
from linearity is only 0.76 percent, although the 
mobilities of Na+ and H;O* stand in the ratio 
of one to seven. The phenomenon we have ob- 
served is of considerably greater magnitude. 

We are indebted to the Rockefeller Foundation 
for the loan of the heavy water samples used in 
these measurements, and wish to extend our 
thanks for this loan. 


*1 This deviation from linearity was predicted qualita- 
tively by Dr. Otto Halpern of New York University in a 
private conversation in April, 1934. 

22 Onsager and Fuoss, J. Phys. Chem. 36, 2758 (1932). 
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The Vapor Pressures of the Isotopic Forms of Water 


MiLton H. Want AND Haron C. Urey, Department of Chemistry, Columbia University 
(Received May 3, 1935) 


The ratio of vapor pressures of H2O" and HDO" and of 
H,O"* and H,O'*8 have been determined at four tempera- 
tures between 11.25°C and 46.35°C. The ratio of vapor 
pressures of HDO"™ are very nearly the geometric mean of 
the vapor pressure of H.O"* and DO". The vapor pressure 
of H,O'8 varies from approximately 0.8 to 1.2 percent less 
than the vapor pressures of HO" for temperatures be- 


tween 46.35°C and 11.25°C. The application of these to the 
separation of the hydrogen and oxygen isotopes by distilla- 
tion is discussed. It seems possible that the hydrogen iso- 
topes can be separated by the fractional distillation of 
water, but that the separation of the oxygen isotopes by 
this method will be considerably more difficult. 





T has previously been observed!:? that both 
the hydrogen and oxygen isotopes are con- 
centrated by distillation of water, indicating 
that the lighter isotopic molecules have the 
higher vapor pressures. In particular the vapor 
pressure of D,O has been investigated by Lewis 
and Cornish,” and at the boiling point, the H2O 
has approximately 5 percent greater vapor 
pressure, with larger values at lower tempera- 
tures. It has been assumed that the vapor pres- 
sure of HDO should be the geometric mean of 
the vapor pressures of H2,O and D,O.* It was 
found that an experimental determination of the 
HDO vapor pressures and also the ratio of vapor 
pressures of H,O" and H,O!8 would be advisable 
to see whether the distillation of water was 
feasible for the concentration of either or both 
the hydrogen and oxygen isotopes. Since this 
work has begun, Wynne-Jones** has published a 
note investigating this same question by a 
method similar to that described in this paper, 
and finds that his results are inconsistent with 
the assumption that the vapor pressure of HDO 
is the geometric mean of H.O and D.O. 

It is impossible to determine directly the vapor 
pressure of HDO since it probably can never be 
secured in the absence of H,O and D.O. Hence 
some indirect method is necessary for the de- 
termination of its vapor pressure. The fractional 
distillation of dilute solutions of deuterium 
water can be used for the determination of the 
vapor pressure of HDO if it is assumed that the 


1E. W. Washburn and E. R. Smith, J. Chem. Phys. 1, 
426 (1933). 

*G. N. Lewis and R. E. Cornish, J. Am. Chem. Soc. 55, 
2616 (1933). 

*B. Topley and H. Eyring, J. Chem. Phys. 2, 217 (1934). 

* W. F. K. Wynne-Jones, J. Chem. Phys. 3, 197 (1935). 


ideal solution laws are obeyed. In order that this 
method can be used, it is necessary that the 
distillation be carried out under such conditions 
that the Rayleigh distillation formula holds, 
which means that there must be no refluxing, and 
spattering of the liquid surface must be avoided. 
Refluxing would give a fractionation factor too 
high, and spattering of the liquid to the walls 
would have the effect of making the fractiona- 
tion factor too low. The apparatus used was 
designed to avoid these difficulties. 


EXPERIMENTAL 


The apparatus was constructed so that the 
distillation of water took place under its own 
vapor pressure at any desired temperature. 
Three liters of half-percent heavy water were 
placed in a five-liter flask which was fitted with a 
mercury-sealed stirrer to insure a uniform tem- 
perature and composition throughout the liquid 
during distillation. The vapor was condensed in 
another large flask immersed in an ice bath. The 
entire apparatus was barostated by means of a 
manometer regulator which actuated a relay, 
which in turn operated a vacuum pump. With 
this arrangement, the pressure variation was less 
than one millimeter. The vessel was heated by 
means of a lamp, arranged so that the evapora- 
tion took place as fast as possible without ebulli- 
tion, in order to avoid spray to the walls of the 


‘flask. At the lower temperatures (11.25°C) no 


refluxing could take place, since the entire 
apparatus was at room temperature. At the 
higher temperatures, refluxing was avoided by 
encasing the distilling flask and stirrer in an air 
bath maintained at a temperature about 10° 
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higher than the temperature of distillation. 
After the water had been degassed and brought 
to temperature, the distillation proceeded very 
smoothly ;—there was no bubbling or spraying 
from the surface. The temperatures were calcu- 
lated from the vapor pressure of water and were 
checked approximately with a thermometer 
immersed in the liquid. 


MeETHODsS OF ANALYSIS 


The quantity of water used for the distilla- 
tion was determined volumetrically, while the 
residue, which usually approximated 50 cc, was 
distilled under vacuum into a trap and accu- 
rately weighed. The density of the original 
materials and of the residues was determined by 
the pycnometer method. The pycnometers had a 
capacity of approximately 30 cc and the densities 
were reproducible to within several parts per 
million. For each sample the density was the 
average of three or more determinations. These 
measurements gave the total increase in density 
due to both the hydrogen and oxygen isotopes. 

The samples were then made alkaline by the 
addition of sodium, and electrolyzed. The oxygen 
of the sample was burned with an excess of 
deuterium-free hydrogen obtained by the elec- 
trolysis of light water,‘ to avoid any fractiona- 
tion of the hydrogen isotopes which would cause 
unknown variations in the densities. The density 
of these samples was determined by a pressure- 
controlled float similar to that described by 
Gilfillan.® The constant of the float used was 0.51 
p.p.m./cm Hg and the results were reproducible 
to several tenths of a part per million. 

The ratio of vapor pressures was calculated 
from the density measurements by use of the 
Rayleigh distillation formula, namely, 


(N, ‘Ny)*/@-D= W,/ W, 


which in this form is applicable to dilute solu- 
tions. 


RESULTS 


Table I lists the temperature, the densities 
of the initial and final materials, the densities of 


*L. A. Weber, M. H. Wahl and H. C. Urey, J. Chem. 
Phys. 3, 129 (1935). 
5E. S. Gilfillan, J. Am. Chem. Soc. 56, 406 (1934). 
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Fic. 1. Plots of the logarithm of the ratio of vapor 
pressure against the reciprocal of the temperature. The 
scale at the left applies to the HO-HDO plot and that 
at the right, the H.O'-H,0!8 plot. The cross represents 
data of Lewis and Cornish. 


the initial and final materials after removal of all 
deuterium, the factor by which the sample was 
reduced during distillation, together with the 
ratio of vapor pressures calculated from the 
Rayleigh distillation formula. 

In the first experiment, all the analyses were 
made with pycnometers. The heavy oxygen 
isotope was removed by passing COs» through the 
sample and the difference in density determined 
with pycnometers. This may cause an error in the 
ratio H,O*/H,O' since the total increase in 
density is only several times the precision of 
pycnometer measurements. A control sample was 
treated with CO, to avoid any possible errors 
due to variations in the ratio of oxygen isotopes 
in the COz used. 

The mean deviations of all pycnometer meas- 
urements is 1.7 p.p.m., which changes the vapor 
pressure ratio of H2O to HDO by 0.2 percent. 
Similarly, the mean deviation of all float meas- 
urements is 0.27 p.p.m., which changes the vapor 
pressure of H,O" to HO" by 0.05 percent. The 
discrepancy between the results of Wynne-Jones 
and those recorded here is larger than could be 
caused by errors in our analyses. 

In Fig. 1 the logarithms of the ratios of the 
vapor pressures are plotted against the reciprocal 
of the temperature. The units to the left of the 
diagram represent the HXO-HDO ratio and those 
to the right the H2O"*-H,O' ratio. The solid lines 
represent a least-squares solution of the points, 
which was deemed necessary because of the 
scattering of the points, particularly in the case 
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VAPOR PRESSURES OF WATERS 


TABLE I. 











Sp. Gr. after removal 


Specific gravity 


of deuterium 


Moles initial material 








Initial 
material 


Final 
material 


Initial 


Temp. material 


16 
Puon Px.o 


Px20" 





Final 


material 


Moles final material Prop 





M25°C 
11.25 
23.00 
35.6 
46.35 


1.000549 
1.000522 
1.000540 
1.000512 
1.000524 


1.000759 
1.000733 
1.000714 
1.000661 
1.000677 


1.000024 

1.000028, 
1.000028, 
1.000028, 
1.000028; 


1.000038 

1.000040, 
1.000037 ¢ 
1.0000365 
1.000036. 


58.917 
63.208 
52.002 
63.039 
82.081 


1.084 
1.089 
1.076 
1.066 
1.062 


1.014 
1.012 
1.0083 
1.008. 
1.0077 








The specific gravities in column 4 are greater than unity because of an increased concentration of the oxygen isotopes in the electrolytic water 
used. This was approximately 4 percent heavy water supplied by the Ohio Chemical and Manufacturing Company. 


of the H,O'—-H,O'* values. The curves when 
extrapolated to 100°C give values of 1.026 for 
H.O/HDO and 1.003 for H,O'*/H,O0'8. 

The extrapolation of the ratio of the vapor 
pressure of H,O' and H.O'* to 100°C cannot be 
expected to give a precise value, since the error 


in the individual measurements is so great and 


the number of measurements too small. From 
Lewis and Cornish” results, together with their 
statement that their column contained about 40 
theoretical plates, it is possible to calculate the 
ratio of vapor pressures, which is found to lie 
between 1.004 and 1.006, the average being 
about 1.005. It is plotted in the diagram, and as 
can easily be seen, it is not inconsistent with 
our ratio of vapor pressures at lower tempera- 
tures and indeed may be a better value for this 
temperature than our extrapolated value. 

From the slopes of the curves, the following 
values of A// were obtained : AJ7= AIT yp0 — Al n,0 
=138 cal. and A/J=AHy,o'%—AlTy,0'°= 26.4 
cal. The former value is approximately half the 
value given by Topley and Eyring,’ for Al/p,o 
—AlTy,o. If the vapor pressure of HDO is the 
geometric mean of the light and heavy molecule, 
the difference in the heats of vaporization of 
HDO and H.O should probably be one-half the 
difference in heats of vaporization of HO and 
D.O. The difference in heats of vaporization of 
HO" and H,O'8 cannot be regarded as very 


TABLE II. 


precise, and in fact may be in error by as much 
as a factor of two. Using Lewis and Cornish? 
value again, we would estimate this difference 
in heats of vaporization as considerably smaller 
than 26.4 secured from our experiments. 

Using the formula given by Topley and 
Eyring,* 


Pp,o/Pu,0= 1.35¢259/#7 


one can calculate the ratio of vapor pressures of 
the heavy and light molecules. In Table II are 
recorded this ratio of vapor pressures, both the 
geometrical and arithmetical means of the 
heavy and light molecules, and the experi- 
mentally determined value of the vapor pressure 
of HDO. 

It can be seen that the experimental values of 
the vapor pressure of HDO are in better agree- 
ment with the geometric mean than with the 
arithmetical mean of the light and heavy mole- 
cules. The measured values are consistently 
smaller than the corresponding geometrical 
mean, which is somewhat less than the arith- 
metical mean. However, the precision of the 
results is hardly sufficient to be certain that the 
disagreement of experimental values with either 
the geometric mean or arithmetical mean is out- 
side the limits of experimental error. 

In Table III are tabulated the vapor pres- 
sures of the various isotopic species of water at 


TABLE III. 








Pp.o 
Po 


0.853 
853 
.869 
885 
.897 


P20 +Pp.0 
2 


Temp. 


£.25°C 
11.25 
23.00 
35.6 
46.35 


(Px.0Pp.0)* 


9.24 
9.24 
19.58 
40.92 
72.93 


Pupo 


9.23 
9.18 
19.51 
40.80 
72.83 





9.27 
9.27 
19.63 
41.00 
73.04 


Pp,o% Pupo" Py,0% 





8.53 
8.53 
18.25 
38.50 
69.07 
722.8 


9.23 
9.18 
19.51 
40.80 
72.53 
740.7 


9.86 
9.89 
20.82 
43.15 
76.41 
757.7 
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the temperatures at which the experiments were 
conducted. The values for 100°C were extra- 
polated from the graph. 


DIsCUSSION 

As first pointed out by Lewis and Cornish, the 
fractionation factor for the separation of deu- 
terium by the distillation of water should be 
equal to the square root of the ratios of the 
vapor pressures of H.O and D.O provided the 
vapor pressure of HDO is the geometric mean of 
the vapor pressures of the other two waters, and 
provided the equilibrium constant for the 
exchange reaction, HO0+ D.2O=2HDO is equal 
to 4. The first of these conditions, as shown 
by our data, is very closely true. The second, 
according to the calculations of Topley and 
Eyring,® is approximately 3.4 instead of 4.0 at 
the boiling point. Our calculation of the ratio 
of vapor pressures of H2O and HDO is made on 
the assumption that there are no D.O mole- 
cules present. If we assume the value of the 
equilibrium constant for the reaction, HXO+D.O 
= 2HDO to be 3.26 at 25°C, and carry through 
an exact calculation for the fractionation factor, 
assuming that the vapor pressure of HDO is the 
geometric mean of the vapor pressure of HO 
and D,O, we find that the ratio of vapor pres- 
sures which we have listed may be in error by 
about 1 part in 10,000, an amount which is 
negligible considering other sources of error. Thus 
the fractionation factor should be very nearly 
that assumed by Lewis and Cornish.? With frac- 
tionation columns of the type devised by a num- 
ber of workers, particularly by Fenske® and his 
co-workers, it should be possible to make a very 
effective separation of the hydrogen isotopes, 
providing no unforeseen difficulties develop. 
100 theoretical plates should give a fractionation 
factor of approximately 10 for the separation of 
hydrogen isotopes by the distillation of water, 
and it thus appears that the distillation method 
may be a serious competitor of the electrolytic 
method. Recently Cryder, Fenske and _ co- 
workers’ have succeeded in securing a twofold 
concentration of deuterium by the distillation 
of water in a packed column containing about 
50 theoretical plates. The concentration was not 


6M. R. Fenske, C. O. Tongberg and D. Quiggle, Ind. 
and Eng. Chem. 26, 1169, 1213 (1934). 

7D. S. Cryder, M. R. Fenske, D. Quiggle, C. O. Tong- 
berg and A. R. Lux. Presented Spring Meeting of Am. 
Chem. Soc., 1935, New York City. 
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as large as anticipated for the height of column 
used ; however, the results are very encouraging. 

Our results are in disagreement with the ex- 
periments of Wynne-Jones** who investigated 
the fractionation ratio in distillation, using 
approximately 10 and 50 percent heavy water. 
He finds a fractionation ratio at 25° of 1.10, in 
disagreement with the fractionation ratio of 
1.076 which we find at 23°C. The method which 
he used is similar in principle to that used here, 
and it would appear that the precautions we 
have taken were also taken by him. However, we 
are unable to see any source of error in our experi- 
ments. His report of his experiments is so brief 
that it is impossible to judge any possible 
sources of error in his work. It seems possible 
that the temperature of his water was not at the 


temperature of his bath but that it might be 


somewhat lower due to the cooling effect of 
evaporation. It does not seem likely that the 
temperature could have been sufficiently low to 
have accounted for the discrepancy. If his value 
should prove to be better than ours, we might 
expect that a similar error in the ratio of the 
vapor pressures of H,O" and H,O!* might occur 
in our results, and that therefore this ratio 
might be larger than that recorded. 

The separation of the oxygen isotopes by 
distillation is not so promising unless consider- 
ably more effective fractionation columns can be 
devised. 100 theoretical plates in this case would 
give a fractionation factor of approximately 1.6, 
not a particularly favorable fractionation factor. 
At the lowest temperatures investigated, the 
fractionation factor would be approximately 
3.3 for a 100-plate column. This, of course, is 
quite an effective fractionation factor, but it 
seems probable that the capacity of a still 
column would decrease rather seriously at these 
lower vapor pressures, namely approximately 10 
mm of mercury. In this case it may be necessary 
to make use of the exchange equilibria discussed 
by Urey and Greiff.® 
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The absorption spectra and vapor pressures of the hydrogen and deuterium iodides have been 
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measured. The absorption of the hydrogen compound occurs at longer wavelengths than does 
the deuterium iodide. The vapor pressures are only slightly different, the deuterium iodide 
having a greater value at a given temperature. The theoretical reasons for these differences are 
discussed, and equations for the vapor pressures are given as determined by this and other 


work. 





HE absorption spectra of the hydrogen 
halides have been the subject of consider- 
able work.' The results have led uniformly to the 
belief that the absorption occurs from the zero 
(v=4) vibrational level of the normal !Z state to 
a higher, unstable electronic level whose vibra- 
tional states are not quantized, at least in the 
region to which absorption occurs. The result 
of such a process would be the continuous ab- 
sorption with a more or less vague short wave 
limit which has been, observed. If this picture 
be the true one, the corresponding deuterium 
halide should show marked deviation in its 
absorption spectrum, due both to.the difference 
in zero point energy'* of the ground states and to 
the different eigenfunctions characteristic of 
these levels. The heavier deuterium atom should 
penetrate less into those regions which are not 
permitted on classical theory. We have therefore 
compared hydrogen and deuterium iodides in an 
effort to discover if the results agree with the 
commonly accepted theory. 

Pure hydrogen iodide was prepared by passing 
hydrogen and iodine over a heated platinum 
catalyst, and then distilling the resulting product 
from a trap in a carbon dioxide snow-acetone 
mixture to one in liquid air until, after 3 distilla- 
tions, a colorless product was obtained. In 
order to avoid possible contamination from 
hydrogen adsorbed on a catalyst, the deuterium 
iodide was obtained by heating an excess of 
iodine with deuterium gas in a 2-liter bulb which 





! Tingey and Gerke, J. Am. Chem. Soc. 48, 1838 (1926); 
Dutta, Zeits. f. Physik 77, 404 (1932). 

1 This effect has been found by Franck and Wood in the 
absorption of the two water vapors in the Schumann region 
(Phys. Rev. 45, 667 (1934)). The difficulties of quantitative 
absorption measurements in this region made it impossible 
for them to find more than a long wave limit shift, which 
Was qualitatively in agreement with the predicted differ- 
ence in zero point energies. 
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had been previously baked out for some hours in 
vacuo at 450°C. The bulb was then sealed on and 
opened into the distillation train in the manner 
already described by Rittenberg and Urey.’ 

After the distillations the deuterium iodide 
which was frozen in a liquid air trap, the absorp- 
tion cell, and the quartz spiral manometer used 
to measure the pressure were sealed off from the 
pumps. The gas under measurement thus came 
into contact only with glass and quartz in a 
dark room and showed no measurable decomposi- 
tion even after the illumination necessary to 
obtain the absorption spectra. 

The deuterium gas was more than 99 percent 
pure as measured by combining a sample with 
oxygen and determining the refractive index of 
the water so-formed in an interferometer. 

The quantitative absorption spectra were ob- 
tained by the use of a Judd-Lewis sector photom- 
eter, a Hilger quartz spectrograph (E-3) and the 
continuous spectrum from a condensed spark 
between tungsten electrodes under distilled 
water. The fused quartz absorption cell had a 
12.41 cm optical path between the inside faces 
of the plane parallel ends. Several different 
pressures of gas were used and at each pressure 
eight pairs of exposures were made with eight 
different settings of the comparison beam. The 
frequency number where each pair of exposures 
possessed equal intensity was ascertained under 
a low power comparator. The molecular absorp- 
tion coefficient was calculated from the equation 
a,=(760/PXL)-e where P was the observed 
pressure in mm, L the length of the gas path in 
cm, and e, the setting of the comparison beam, 
was the logarithm of the ratio of the intensity of 
transmitted light to the intensity of the incident 


* Rittenberg and Urey, J. Am. Chem. Soc. 56, 1885 
(1934). 
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TABLE I. 


Mol. abs. Mol. abs. 
coeff. coeff. 








Frequency X 10-2 Frequency X 1072 





305.6 0.012 
309.7 024 


322.5 0.109 
334.7 .219 
341.9 .328 
346.2 437 
349.5 546 
352.0 656 
355.5 -765 
357.6 .874 


311.6 .036 
315.2 048 
316.9 -061 
318.6 O85 
321.7 -097 


0.008 334.8 .125 
016 336.8 
023 338.4 
.030 340.4 
.038 341.2 
.046 346.5 
.053 350.8 
.050 352.9 
061 355.3 
075 358.2 
.100 360.3 











light. The absorption of the empty cell was 
ascertained and a correction for its absorption 
was applied to the values obtained for the gases. 
The temperature of the gas in the absorption 
cell was approximately 25°. 

The value of the absorption coefficient a, is 
given in Table I for different frequencies and 
plotted in Fig. 1. The value of a, is seen always 
to be higher for hydrogen iodide than for 
deuterium iodide at a given frequency as may be 
predicted. 

In order to obtain a more quantitative com- 
parison between the absorption curves reference 
should be made to Fig. 2. Here are plotted the 
potential energy curve of the 'Y state which is 
the same for HI and DI, treated as an harmonic 
oscillator, and the two eigenfunctions of the zero 
vibrational states, as well as a hypothetical 
completely repulsive upper state. In absorbing 
radiation of any frequency v the molecule in 
the zero vibrational state is raised to some point 
on the upper potential energy curve, for which 
there will be an eigenfunction for every value 
of the potential energy. A typical example is 
shown by the dotted line. According to the 
theory suggested by Franck and worked out in 
terms of the wave theory by Condon,’ the prob- 
ability of a transition will be governed by the 
probability of a molecule in the ground state 
being in such a position as to be able to jump 
to that portion of the upper curve required by 
the frequency of the radiation, without change 


® Condon, Phys. Rev, 28, 1182 (1926); 32, 858 (1928). 
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in nuclear separation, r. Since the eigenfunction 
of HI has larger values at greater values of r than 
has that of DI the absorption of the former should 
extend to lower frequencies. Further, the absorp- 
tion curves should at the point of maximum ab- 
sorption differ by the difference in zero point 
energies (330 cm~'). This can readily be seen to be 
indicated by the results of Fig. 1. 

Stueckelberg and Gibson, Rice and Bayliss’ 
have worked out theoretical expressions to indi- 
cate the value of a, for cases similar to that under 
consideration. The resulting expressions have 
been applied to the continuous absorption of 
oxygen and iodine with considerable success. 

In the present case, however, two important 
constants are lacking which are necessary for 
the use of the equations of Gibson, Rice and 
Bayliss. We do not know the frequency of 
maximum absorption, nor the total width of the 
absorption bands. The band of absorption of the 
halides apparently extends far into the ultra- 
violet with no indication on our plates of any 
decreasing absorption on the short wavelength 
side. This is probably due to factors which make 
the absorption under consideration different 
from the case of iodine, which they studied. The 
reduced mass is much smaller and the upper 


4 Stueckelberg, Phys. Rev. 42, 518 (1932); Gibson, Rice 
and Bayliss, Phys. Rev. 44, 193 (1933). 
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curve possibly steeper in the region of strong 
absorption, giving rise to an unusually broad 
absorption band. 

It is possible to obtain a rough interpretation of 
our data using the expression of Stueckelberg, 
though here, too, the lack of the complete curve 
makes a fit somewhat dubious. The theoretical 
expression also requires that the upper curve 
have a constant slope, which is probably not 
true for an absorption which is occurring to an 
upper curve when it is approaching the dissocia- 
tion energy level. A fit of the data for HI gives 
as a slope of the upper curve 2.3 10° cm A-. 
The same constants applied to DI lead to a value 
approximately the same, but the requirement 
that »/c: y should give a straight line* is no 
longer fulfilled. This probably indicates a lack of 
linearity in the upper potential energy curve. 
The value of the slope, however, is a very reason- 
able one, being about 1/5 as large as that found 
for oxygen. 


* Reference 4a, p. 522. 
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All of the recent work has also supported the 
picture of the absorption process as given, but 
there are some points which should perhaps be 
clarified. Rollefson and Booher® have shown that 
the absorption extends to longer wavelengths 
than can be accounted for if the decomposition 
products are an excited iodine atom (?P;) and a 
normal hydrogen atom. They suggest that two 
separate processes are occurring, the one just 
mentioned and another giving the two normal 
2P3/2 and *S, states. It may be shown that the 
absorption observed by these authors at long 
wavelengths (a, =10-', assuming that the value 
of e is about 10 where blackening of the plate is 
no longer visible) is of the same order of magni- 
tude as the relative populations of the first and 
the zero vibration levels. Consequently it might 
appear that absorption might still be due to a 
single electronic transition giving a *P,; iodine 
atom with the long wave part accounted for by 
transitions from the level v=1 in the lower state. 
However the absorption occurs to a wavelength 
(3900A) which is shorter than could be produced 
by such a process and the occurrence of a transi- 
tion to a state giving normal atoms seems to be a 
necessary conclusion. That the latter is the sole 
process, and involves a completely repulsive 
upper state is quite plausible. 

In a more recent publication, Dutta and Deb*® 
have obtained continuous emission spectra with 
long wave limits from the hydrogen halides. 
As has already been pointed out by Urey and 
Bates’ this is not the type of spectrum to be ex- 
pected from such a process, which should have 
no real long wave limit. The latter investigators 
looked for such spectra in the hydrogen halogen 
flames, but decided that those actually observed, 
having no long wave limit and extending beyond 
the red, were characteristic of a recombination 
of two halogen atoms rather than of hydrogen and 
halogen. This is confirmed by the thermal 
studies of Kondratjew and Leipunsky.* The 
hydrogen halogen recombination spectra should 
obviously be similar to these, but extend to 
shorter wavelengths. 


— and Booher, J. Am. Chem. Soc. 53, 1728 
1931). 

6 Dutta and Deb, Zeits. f. Physik 93, 127 (1934). 

7 Urey and Bates, Phys. Rev. 34, 1541 (1929). 
8 Kondratjew and Leipunsky, Zeits. f. Physik 50, 366 


(1928) 
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THE VAPOR PRESSURE OF HYDROGEN AND 
DEUTERIUM IODIDE 


The pressures were measured with the ap- 
paratus described above. A five junction copper- 
constantan thermocouple, accurate to +0.1°, 
was employed for the temperature measure- 
ments. In the first series of experiments, with the 
small amounts of the iodides originally prepared 
for the study of the absorption spectra, repro- 
ducible vapor pressures for both solid and liquid 
states were obtained, but the triple points could 
not be determined with sufficient accuracy. With 
the liquids equilibrium was established very 
rapidly, but the solids produced reproducible 
pressures only after standing for considerable 
lengths of time at constant temperature. At the 
lowest temperatures, produced by carbon di- 
oxide-acetone mixtures, at least two hours were 
required. Constant temperatures at other points 
up to the boiling points of the iodides were 
maintained by means of a series of freezing 
baths. In the entire series of measurements the 
tube which served as the thermocouple well 
was fixed to the trap containing the hydrogen 
or deuterium iodide in such a way as to measure 
the temperature of the bath in the region nearest 
to the sample. 

The triple points were determined by means 
of heating curves. Larger quantities (about 20 g) 
of the iodides were prepared and the apparatus 
was rebuilt so as to effect the direct immersion 
of the thermocouple well into the hydrogen 
or deuterium iodide. With this arrangement, the 
triple points were readily determined, but satis- 
factory vapor pressures could no longer be ob- 
tained. 

The new data presented here for liquid hydro- 
gen iodide are in good agreement with values 
available in the literature.° For the solid the 
agreement is less satisfactory. The results of 
this and other investigations are perhaps best 
discussed by comparison with vapor pressures 
calculated from the boiling point, triple point, 
and calorimetric data presented by Giauque and 

Steele, McIntosh and Archibald, Zeits. f. physik. 
Chemie 55, 129 (1906); Maas and McIntosh, Proc. and 
Trans. Roy. Soc. Canada 13, 65 (14); Steele and Bayster, 
J. Chem. Soc. 97, 2607 (1910); Henglein, Zeits. f. Physik 


18, 64 (1924); Miravalles and Moles, Anales soc. espaii fis. 
quim. 23, 509 (1925). : 
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Wiebe" for hydrogen iodide. If the various ex- 
perimental points are placed on a large scale 
graph showing the logarithm of the pressure as a 
function of the reciprocal of the temperature, 
those for the liquid are found to lie on a line 
slightly above the calculated values, while 
those for the solid lie almost without exception 
below the calculated-curve. A line drawn through 
the average of the experimental points for the 
solid intersects the corresponding line for the 
liquid at some distance above the triple point. 
Furthermore, in the 15° temperature range im- 
mediately below the triple point, the ‘“experi- 
mental curve” shows a curvature opposite to that 
required by the heat capacity data. Closer agree- 
ment between measured and calculated pres- 
sures is observed in the neighborhood of the 
carbon dioxide sublimation point. It is evident, 
then, that most of the available data for solid 
hydrogen iodide represent pressures below the 
true equilibrium values. This is consistent with 
our experience that reproducible vapor pressures 
could be obtained with the solid only when con- 
stant temperatures were maintained for con- 
siderable lengths of time. 

We have chosen to express the vapor pressure 
of hydrogen iodide as a function of the tempera- 
ture by means of curves with slopes determined 
by the calorimetric data of Giauque and Wiebe 
and passing through our measured value of the 
pressure (log p=2.579) at the triple point 
(222.3°). The equation thus obtained for the 
liquid represents the direct experimental data 
quite accurately. In addition to the assumption 
that the heat capacity of the vapor is given by 
5/2 R, the following calorimetric data, taken 
from the measurements of Giauque and Wiebe, 
were employed in setting up the equations: 
AF1937.75 (vaporization) 


= 4724 cal. (at the boiling point), 
AHo22.3 (fusion) = 686 cal. (at the triple point), 


Cp(l) =19.11—0.021T; Cp(s) =5.72+0.0297. 
The equations follow: 


log p(HI, 2) = 26.119+0.002293T 
—7.111 log T—1636/T, (1) 


1 Giauque and Wiebe, J. Am. Chem. Soc. 51, 1441 
(1929). 
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log p(HI, s) = 10.493 —0.003167T 
— 0.377 log T—1406/T. (2) 


The corresponding relations for deuterium 
iodide are obtained by comparing, as shown in 
Table II, the measured vapor pressures for 


Takle II. Vapor pressure of solid and liquid hydrogen and 
deuterium todides. 








log p(HI1)(calc.) log p(H1I)(exp.) log p(DI)(exp.) 





(a) Liquid 

813 2.81 2.823 

.675 2.685 

58 2.586 

.52 supercooled liquid 2.535 
(b) Triple point 

(2.579) (2.57 


(3so2uh 
2.562(1) 


9) 
2.584 
2.560, 2.567 


(c) Solid 











deuterium iodide with those calculated for 
hydrogen iodide at the same temperatures. The 
table also shows the comparison between meas- 
ured and calculated values for hydrogen iodide. 

In addition to the tabulated results, a number 
of experiments were performed at temperatures 
slightly below 200° in which equilibrium was ap- 
proached from either side. In every case the 
pressure could be shown to be approaching a 
value consistent with the above results. For the 
solid, only the results obtained where constant 
temperature was maintained for sufficient peri- 
ods of time are recorded in the table. The 
measured pressures for deuterium iodide at the 
triple point, obtained with-the 20 g batch of 
material are obviously low. The true value, by 
interpolation from the other measurements on the 
liquid, is 2.572. 

Examination of the table shows that the 
logarithm of the vapor pressure of liquid deu- 
terium iodide may be obtained by adding 0.01 
to the calculated value for hydrogen iodide at 
the same temperature. The corresponding addi- 
tion for the solid is not so evident from the data, 
but is set at 0.012 to take account of the differ- 
ence between solid and liquid hydrogen iodide 
at the triple point of deuterium iodide. There- 
fore, for deuterium iodide, 
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log p(DI, 2) =26.129+0.002293T 
—7.111 log T—1636/T, (3) 
log p(DI, s) = 10.505 —0.003167T 
—0.377 log T—1406/T. (4) 


Substituting for p= 760 mm in Eggs. (1) and (3) 
gives for the boiling points most consistent with 
our experimental data, 237.5° for hydrogen 
iodide and 237.0° for deuterium iodide. We are 
at a loss to account for the deviation of our 
boiling point for hydrogen iodide from that of 
Giauque and Wiebe (237.75°). Our probable 
error in temperature measurement is of the order 
of +0.1°, the deviation of hydrogen iodide from 
the ideal gas laws might be expected to operate 
in the reverse direction, so as to widen the 
discrepancy, and the correction of our pressures 
to mm of mercury at 0°C would raise our boiling 
point only 0.1°. 

The conclusion expressed in Eqs. (3) and (4), 
to the effect that the ratio of the vapor pressures 
of the two hydrogen iodides is a constant over the 
entire range for either the solid or the liquid 
states, can be true only if the heats of vaporiza- 
tion, fusion and sublimation of the two sub- 
stances are nearly identical at the same tempera- 
tures. The difference in heats of vaporization, 
which play the most important role of the three 
mentioned, can be estimated in the neighbor- 
hood of the boiling point by an application of 
Trouton’s Rule. 

Trouton’s rule can be applied to the hydrogen 
iodides with greater confidence than to water or 
to the other hydrogen halides. The lower electric 
moment and the numerical value of Trouton’s 
constant, 19.9 for HI, suggest that the iodides 
should behave as ‘“‘normal’”’ liquids. The rule is 
applied in the form suggested by the Hildebrand 
modification" in which the heat temperature 
quotients at temperatures of equal equilibrium 
molal volume for the vapors are set equal. 
From Eq. (4).it is found by approximation that 
at 269.93° the molal volume of vapor at equilib- 
rium with deuterium iodide is the same as that 
for hydrogen iodide at its boiling point. Whence, 
for deuterium iodide, 


AFT o36.93 (vap.) =4713 cal. 


1 Hildebrand, J. Am. Chem. Soc. 37, 970 (1915). 
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and at the boiling point of hydrogen iodide, the 
value may be estimated from the heat capacity 
of hydrogen iodide to be 4709 cal. as against 
4724 cal. for hydrogen iodide. The effect of this 
difference upon the vapor pressures is within 
the limit of accuracy of our measurements. 
Vapor pressures at the boiling point calculated 
from a single pressure at the triple point and 
using, respectively, the two heats of vaporiza- 
tion differing by fifteen calories would be differ- 
ent by 0.7 mm, or the calculated boiling points 
would differ by 0.05°. 

By reference to the summary presented by 
Claussen and Hildebrand” in connection with 
their vapor pressures of the hydrogen and 


12 Claussen and Hildebrand, J. Am. Chem. Soc. 56, 1820 
(1934). 
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deuterium fluorides, it is seen that deuterium 
iodide differs from hydrogen iodide in the same 
direction as deuterium fluoride from hydrogen 
fluoride, but to a considerably smaller extent. 
The greater vapor pressure of the deuterium 
compounds in these two cases cannot be as- 
cribed to the same cause. The fluoride is known 
to have properties differing in general from the 
other halides. This has been supposed to be due 
to the more ionic nature of the hydrogen fluorine 
bond. The iodide should according to recent 
theories be at the other extreme.” A possible 
explanation seems to be that the polarity of the 
deuterium iodide should be smaller, because of 
the smaller anharmonicity factor operative in 
the region of zero point energy. 


13 Mulliken, J. Chem. Phys. 2, 782 (1934). 
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The Crystal Structure of Lepidocrocite 


F, J. Ewinc, Gates Chemical Laboratory, California Institute of Technology 
(Received February 27, 1935) 


The unit of structure, space group symmetry and de- 
tailed atomic arrangement of lepidocrocite are determined 
from x-ray data, and the existence and location of hydro- 
gen bonds in the structure are established from considera- 
tions of interatomic distances. The orthorhombic lattice is 
end centered on (100), and has the axes ao=3.87A, bo 
12.51A, and co=3.06A. The space group is V;!’—Amuam. 
The parameters are determined as upe= —0.332, uo 
= +0.282, and u9q = +0.075. The structure is described 


INTRODUCTION 


EPIDOCROCITE, FeO(OH), is described 
by Posnjak and Merwin! as forming or- 
thorhombic crystals with the axial ratios 
0.64 : 1: 0.43. The crystals are red and trans- 
parent, and form thin plates tabular on (010), 
slightly striated in the direction of the c-axis. 
The monohydrates of iron and aluminum 
oxide, with the general formula MHOz, exist in 
two modifications. The structure of the modifica- 


as consisting of iron-centered oxygen octahedra joined by 
the sharing of edges into two-dimensionally infinite layers, 
with the successive layers held together by hydrogen 
bonds. The relationship between the diaspore-goethite 
structure and the béhmite-lepidocrocite structure is dis- 
cussed from the standpoint of the coordination theory, and 
found to depend on two alternative ways of satisfying the 
electrostatic valence rule. 


tion corresponding to diaspore, AlHO:, and 
goethite, FeHOs, is known,? but no previous 
investigation of the structure of the second 
modification corresponding to lepidocrocite, 


1E. Posnjak and H. E. Merwin, Am. J. Sci. 47, 311 
(1919). 

2 M. Deflandre, Bull. Soc. Franc. Mineral. 55, 140 (1932); 
K. Takane, Proc. Imp. Acad. Tokyo 9, 113 (1933); 5. 
Goldsztaub, Comptes rendus 195, 964 (1932); F. J. Ewing, 
J. Chem. Phys. 3, 203 (1935). 
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FeO(OH), and béhmite,? AIO(OH), has been 
made. 

In the present investigation the unit of 
structure, space group symmetry, and detailed 
atomic arrangement of lepidocrocite are de- 
termined from x-ray data, and existence and 
location of hydrogen bonds in the structure are 
established from considerations of interatomic 
distances. The lepidocrocite structure and its 
relation to the diaspore-goethite structure is 
discussed from the standpoint of the coordination 
theory. 


THE UNIT OF STRUCTURE AND SPACE GROUP 
SYMMETRY 


Oscillation and Laue photographs were pre- 
pared with crystals of lepidocrocite from Eiser- 
feld, Westerfald, Germany. The crystals were 
tabular on b, forming thin plates about 1 mm? 
in area. Oscillation photographs of Mo Ka 
radiation reflected from or transmitted through 
(010) with [001] as the axis gave the value 
3.06A for co on the basis of the intervals between 
side spectra; and gave dp as 1.935, and bo as 
6.2552 on the basis of the pinacoidal reflections. 
Calculated values of md for observed Laue re- 
flections require that m;=n2=2. The correct unit 
cell is therefore given by the orthorhombic axes 


ao=3.87A, bo=12.51A and co=3.06A. 


No first-order Laue reflections were observed 
in which (k+/) was odd, and the lattice is there- 
fore based on To’, end centered on (100). As- 
suming holohedral symmetry for the crystal, the 
possible space groups* are V;,!" to V;”. Observed 
first order Laue reflections such as (520) and 
(580) eliminate V;2° and V,”*. Numerous first 
order prism reflections on the oscillation photo- 
graphs, such as (011), (031), (051), etc., also 
eliminate the space groups V;'* and V,”!. Of the 
two remaining space groups, V,!’ and V,!° 
the former allows only two prism zones to reflect, 
while the latter permits reflections from all 
three. No first-order reflections from planes of 
the type (hol) were observed in any of the 





*M. E. Nahmias, Zeits. f. Kryst. 85, 355 (1933). 

* No evidence was found in support of the statement by 
S. Goldsztaub, Comptes rendus 193, 533 (1931), that 
lepidocrocite possesses the space group symmetry V>}. 
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photographs, although a number of such planes 
were in a position favorable for reflection. The 
space group of lepidocrocite is thus indicated to 
be V;,!7—Amam. 

The assignment of four formulas of FeO(OH) 
to the unit cell gives a calculated density of 
3.96 g/cm*, which is in agreement with the 
observed density of 4.07. 

The eight oxygen atoms in the unit cell must 
be placed on the intersections of the symmetry 
planes to avoid O—O distances of 1.93A or less. 
The four iron atoms cannot be at centers of 
symmetry, because the intensity of the reflection 
(200) is observed much stronger than that of 
(400). The general positions of the atoms, re- 
ferred to a center of symmetry as origin, are 
therefore ; 
4Fe ati, u,0; 3,u,0; 41, u+3,4; 3, 4-u,}4: 
40, at }, m4, 0; etc.: 

40, at 4, ue, 0; etc. 


PREDICTION OF THE ATOMIC ARRANGEMENT 


A structure may be predicted for lepidocrocite 
from consideration of the following arguments, 
which are derived from the rules given in 
Pauling’s* coordination theory of the structure 
of ionic crystals. 


1. The elements of the structure are iron-centered oxy- 
gen octahedra, in which the Fe—O distances are close to 
2.00A, and the O—O distances have values in the neigh- 
borhood of 2.70A. 

2. The electrostatic valence rule is satisfied, whereby the 
anions are effectually neutralized by adjacent cations. 
The two formal ways of distributing the hydrogen ions in 
FeHO, are indicated by writing (OH;)~*’? anions in one 
instance, and distinct O- and (OH7~) anions in the other. 
It has been shown’ that the assumption that all the anions 
are of the form (OH;)~*/? leads to the diaspore-goethite 
structure. It seems logical to assume, therefore, that 
lepidocrocite exists as a second modification of FeHO: by 
virtue of having the alternative distribution of hydrogens. 
The O= and (OH)~ anions corresponding to this second 
alternative are required by the electrostatic valence rule 
to be shared between four and two octahedra. 

3. Edges shared between octahedra are contracted to 
about 2.50A, with other edges being compensatorily 
lengthened. 


The length of the a and ¢ axes suggest the 
probable orientation of the octahedra, since the 
value 3.87A for do is close to the length 3.82A 


¢ Linus Pauling, J. Am. Chem. Soc. 51, 1010 (1929). 
5F, J. Ewing, J. Chem. Phys. 3, 203 (1935). 
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of the diagonal of an undistorted octahedron, 
and since the c axis, 3.06A, can readily be 
associated only with a lengthened edge of an 
octahedron. 

With the character of the sharing and the 
orientation of the octahedra defined, and with 
the additional restriction that an identity period 
of 12.51A be maintained in the direction of the 
b axis, the layer structure indicated in Fig. 1 
immediately suggests itself. The layers are 
formed by the sharing of edges between octa- 
bedra in such a way that the oxygens O,; near 
the middle of the layers are common to four 
octahedra and correspond to O=, while the 
oxygens Oy, on the outer ridges of the layers 
are common to two octahedra and correspond 
to (OH)-. The layers can be regarded as formed 
by the sharing of edges between “double rutile 
strings” that were found as characteristic aggre- 
gates in diaspore. 

The a and c¢ identity distances are contained 
in the plane of the layer and correspond to a 
body diagonal and an edge of an octahedron, 
respectively. The } axis is perpendicular to the 
layers, and its value of 12.51A suggests that it 
represents two layer intervals. The most ap- 
parent way to satisfy this periodicity is to 
translate every other layer by c/2. 

The structure predicted in this manner has 
the same space group symmetry JV,!’ that is 
found for lepidocrocite, which is a strong indica- 
tion that the correct structure has been obtained. 

In this structure the parameters for the 
oxygens are fixed at “;= +0.267 and u2= +0.067 


Fic. 1. A portion of two layers of lepidocrocite, showing 
the iron-centered oxygen octahedra joined to form layers, 
and the layers tied together with hydrogen bonds (drawn 
as tubes). 
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by the requirement that shared edges be 2.50A 
long, while the iron parameter zw is calculated as 
—0.332 on the assumption that the Fe—O 
distances are as nearly equal as possible. These 
parameter values differ only by 0.008 to 0.015 
from the values to be derived from intensity 
data in the next section. 


VERIFICATION OF THE ATOMIC ARRANGEMENT 


Spectra up to the twentieth order were photo- 
graphed by reflecting Mo Ka radiation from the 
pinacoidal plane (010) of a crystal oscillating 
about [001 ] as an axis. The observed intensities 
of reflection were related to the atomic arrange- 
ment by means of the formula 


1+ cos? 20 
"Need 
2 sin 20 


w F2e—28 (sin Od)? 


in which C is a constant, (1+cos? 20)/2 sin 26 is 
the polarization factor, w is the Ott factor for 
varying time of reflection, e~?2(si= ”/? jis the 
temperature factor, and F has the usual form 


F= > fje2tiGzitkujtte;), 
‘ 7 


The value of B was estimated as 1.00. Values 
of f were taken from Pauling and Sherman’s‘ 
calculated f values for Fe and O. 

The f values for Fe are so much larger than for 
O at large angles of reflection that it is possible 
to evaluate u, the Fe parameter, independently 
of the oxygen parameters u; and u2. This was 
done by writing the intensity formulas concerned 
for a given observed inequality, assigning values 
to the oxygen parameters most favorable to the 
observed inequality, and solving graphically for 
permitted values of the iron parameter. The 
solution of ten inequalities over the range’ 
—0.50<u<—0.250 served to limit u to the 
region —0.3235+0.004. 

The oxygen parameters were varied over the 
ranges 0.065 to 0.085 and 0.265 to 0.285, which 
ranges were suggested by the predicted values, 
while the iron parameter was varied within its 
above limits; the intensities therewith calculated 


6 Linus Pauling and J. Sherman, Zeits. f. Kryst. 81, 1 
(1932). 
7 The space group requires one Fe in this interval, 
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TABLE I. Observed and calculated intensities of reflection on oscillation photograph from (010) with [001] as axis. 
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a Estimated visually. 


were found to agree with observed intensities for 
the values 


u= —0.322+0.001; u,=+0.282+0.005; 
uo = +0.075+0.005. 


The probable errors as given correspond to the 
minimum displacements giving qualitative dis- 
agreement between observed and calculated in- 
tensities. The extent of the agreement obtained 
between observed and calculated intensities using 
these values of the parameters is shown in 


Table I. 


DISCUSSION OF THE STRUCTURE 


The atomic arrangement as calculated from 
intensities differs from the predicted arrangement 
only by shifts of 0.008 to 0.012 in the parameter 
values, and the predicted structure for lepidocro- 
cite is thereby verified. As was mentioned 
previously the layers in this structure have the 
(OH)~- anions on their outer surfaces, and it is 
now possible to show that these anions give rise 
to OHO groups, or hydrogen bonds, holding the 
layers together. 

Hydrogen bonds are postulated in lepidocro- 
cite, as in diaspore,* for the explanation of certain 
O—O distances. In Table II are given the inter- 


atomic distances in lepidocrocite. In each layer 
the O—O distances along shared edges are 2.59A, 
and the other O—O distances, along the lengthened 
octahedral edges, are 3.05A and 3.06A. Between 
the layers the O—O distances are 2.70A, 3.10A, 
and greater. The O—O distance of 2.70A between 
the layers is not only the shortest unshared O—O 
distance in the structure, but is inexplicably so 
unless it is regarded as a hydrogen bond. Alter- 
natively the only forces holding the layers to- 
gether would be the relatively weak ones arising 
from dipole interactions, against which the re- 
pulsive force between the oxygens would be 
effective in raising the distance above 2.70A. 


The value of 2.70A for the OHO distance in 


TABLE II. Interatomic distances in lepidocrocite. 








Atom 


Neighbors 
in same 
octahedron 


Neighbors in 
other octahedra 
of same layer 


Neighbors 
in next 
layer 





Fe 


2 Oy, at 2.00A 
2 O; at 2.00A 
2 O; at 2.01A 


2 O; at 3.06A 
4 On at 3.05A 
4 QO; at 2.59A 
1 Oy; at 2.59A 


2 Oy at 3.06A 
40; at 3.05A 
1 Oy; at 2.59A 


3 Fe at 3.06A 


2 Oy at 3.97A 


20; at 3.97A 
4 Oy at 3.10A 
2 Oy at 2.70A 
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lepidocrocite is, moreover, equal within experi- 
mental error to the value 2.71A found in diaspore. 

The association of one hydrogen with each 
2.70A distance gives the correct stoichiometric 
composition to the crystal. If the hydrogens are 
placed exactly half-way between the bonded 
oxygens they lie at centers of symmetry, with 
the coordinates 000; 300; 033; 333. Fig. l isa 
pictorial representation of the lepidocrocite 
structure, with the hydrogen bonds drawn as 
small tubes between the layers. 

The octahedral units in the layers are con- 
siderably distorted from the shape of a regular 
octahedron, this distortion expressing itself in a 
flattening of that portion of the octahedron that 
is in the interior of the layer, and in a projection 
of the edge formed by the two (OH)~ anions 
away from the layer. The distortion is not as 
marked, however, as in certain other structures, 
notably those of bixbyite and the C modification 
of the sesquioxides.*® 

The O—Fe distances of 2.00 to 2.01A are 
equal to the sum of the radii concerned and do 
not exhibit the effect of Fe—Fe repulsion shown 
in hematite. The value of 2.59A for the length 


of the shared edges in lepidocrocite is in good 
correspondence with the value 2.55A found in 
hematite. 


8 Linus Pauling and M. D. Shappell, Zeits. f. Kryst. 75, 
128 (1930). 

® Linus Pauling and S. B. Hendricks, J. Am. Chem. Soc. 
47, 781 (1925). 
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The four Fe—O bonds around O, have direc- 
tions differing by 5° to 20° from regular tetra- 
hedral directions ; the directions of the four bonds 
around Oy, two to Fe and two to H, differ by 
only 3.5° to 9° from tetrahedral directions. 
While this approximate orientation of the 
hydrogen bond along a tetrahedral direction is 
probably dependent on the specific structure of 
the crystal, it nevertheless admits the possibility 
that the hydrogen bond is in part covalent. 

In connection with the problem of correlation 
between the microscopic structure and macro- 
scopic properties of crystals it may be pointed 
out that the tabular habit of lepidocrocite is such 
that the structural layers run parallel to the 
tabular face, that the excellent cleavage on (010) 
corresponds to the breaking only of the relatively 
weak hydrogen bonds, and that the striations on 
the tabular face run parallel to the “double 
rutile strings” in the layers. The parallel relation 
between the double rutile string and striations on 
developed faces is likewise noted in diaspore, 
and in rutile the striations run parallel to the 
single rutile strings. 

The structural formula FeO(OH) has been 
chosen for lepidocrocite as reflecting in some 
degree the character of the structure, although 
the OH~ anions do not actually exist as such, 
but rather as strings of (—HOHOHOH — ). 

I wish to express my indebtedness to Professor 
Linus Pauling for the interest and helpful advice 
that he has extended to this work. 
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THE DETONATION WAVE 


N 1881, Berthelot and Vieille! and inde- 

pendently Mallard and Le Chatelier? dis- 
covered that when an explosive gas mixture is 
ignited at the end of a uniform tube, the initial 
slow progress of the combustion process is 
quickly accelerated to a high velocity which 
remains constant regardless of the length of the 
tube. This constant velocity combustion process 
is called the “‘detonation wave.’ These early 
studies, and those of others’ which followed, 
show that, once established, the detonation 
wave proceeds at a constant velocity which 
depends upon the constituent gases and their 
proportions in the explosive mixture. On the 
other hand, the speed of detonation is inde- 
pendent of (1) the method of ignition; (2) the 
material forming the tube; (3) the diameter of 
the tube, so long as a certain diameter (6-8 mm) 
is exceeded; (4) the tube end construction which 
may be closed or open; (5) the courtse of the tube 
which may be straight, coiled, or even zig- 
zagged. Further, the influence of moderate varia- 
tions of the initial temperature and pressure is 
small. It appears, then, that the detonation 
velocity is a physical constant of the particular 
explosive gas mixture and as such should be 
amenable to thermodynamic analysis. 


THEORIES OF THE DETONATION WAVE 


In 1881, Berthelot and Vieille’ suggested that 
the speed of the detonation wave was approxi- 
mated by the mean translational velocity of the 
molecules in the burned gases, as calculated from 
the temperature attained in combustion. For 
some gas mixtures, it happens that the velocity 
thus predicted is verified by experiment but 
for many others, the agreement is unsatisfactory. 
In 1893, Dixon‘ combined this idea with a com- 
plicated conception of sound waves but again, the 
1 Berthelot and Vieille, Comptes rendus 93, 18 (1881). 
assy and Le Chatelier, Comptes rendus 93, 145 

3 For references, see Bone and Townend, Flame and 


Combustion in Gases, 1927. 
* Dixon, Phil. Trans. Roy. Soc. A184, 97 (1893). 


computed velocity agreed with experiment for 
some mixtures and not for others. Recently, 
Lewis> has advanced a chain reaction theory 
in which the detonation velocity for the optimum 
mixture of the constituent gases is compared to 
the final velocity of the “‘carrier’”’ in the reaction 
chain. This theory has not been extended to vari- 
ous mixtures of the same gases and so, at present, 
it must be regarded more as an attempt to 
visualize the reaction occurring within the wave. 
In 1899, Chapman® approached this problem 
from the point of view of thermodynamics 
(including the flow of compressible fluids) and 
was able to compute detonation velocities in good 
agreement with those observed for the various 
detonating mixtures of the same gases. Ap- 
parently without knowledge of this earlier work, 
Jouguet,’ in 1905-6, gave an extensive develop- 
ment of this method and achieved a similar 
success. Using the Chapman-Jouguet theory, 
Lewis* has recently recalculated the detonation 
velocity for various mixtures of hydrogen and 
oxygen using modern thermal data and con- 
sidering the dissociation of the combustion 
products. His results, in general, agree very well 
with experiment. 

This Chapman-Jouguet theory depends es- 
sentially on the thermodynamic relations for 
one-dimensional shock waves which were first 
correctly deduced by Hugoniot® in 1887. It is 
important to note that these Hugoniot relations 
fail to define any particular detonation velocity 
but permit velocities ranging from a certain 
minimum to infinity. So far as is generally known, 
no rigorous argument has been made for selecting 
the velocity at which detonation becomes 
stable. To meet this difficulty, Chapman simply 
assumed that. the detonation wave would pro- 
ceed at the minimum velocity and it so happens 
that the velocities thus computed do agree very 
well with experiment. Jouguet also postulated 

5 Lewis, J. Am. Chem. Soc. 52, 3120 (1930). 

6 Chapman, Phil. Mag. 47, 90 (1899). 

7 Jouguet, J. de mathematique 1905, 347; 1906, 6. 

8 Lewis, J. Am. Chem. Soc. 52, 3905 (1930). 


® Hugoniot, J. de l’ecole polytech. 57, 3 (1887); 58, 1 
(1889). 
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stability at this minimum velocity, basing his 
arguments for doing so on a physical interpreta- 
tion of his deduction that, at this speed, the 
burned gases move from the wave front at their 
acoustic velocity. Somewhat similar arguments 
have been advanced by Becker.’® It is the object 
of this paper to approach the problem of defining 
the speed at which detonation is stable by con- 
siderations of available energy and, for this 
purpose, the thermodynamic theory will be 
briefly outlined. 


THE HUGONIOT RELATIONS 


Suppose that a plane detonation wave is 
traveling at its characteristic high and constant 
velocity through a gas mixture that is confined 
within a uniform tube. At any plane 0 ahead of 
the wave, the gases are at rest, entirely undis- 
turbed by the approaching wave front which is 
located between two moving planes, 1 just ahead 
of the wave front and 2 just behind. The notation 
used is as follows: velocity u, pressure P, specific 
volume V, moles of gas mixture NV, temperature 
T, internal energy E, entropy S, available energy 
A. The unburned gases are fed into the wave 
front with properties identified by subscript 1 
and the burned gases are discharged from the 
wave front with properties identified by the 
subscript 2. At any point ahead of the wave, the 
unburned gases have properties identified by the 
subscript 0 and are at rest with respect to the 
tube so that u)=0. Aside from the velocity u 
and the available energy A, the properties at 
plane 0 are assumed to have the same values as at 
plane 1. Since the velocity of detonation is con- 
stant, it is assumed that the conditions im- 
mediately ahead and behind the wave front do 
not vary with time, and further, since the speed 
is high, it is assumed that the effect of heat losses 
is negligible. The principles of continuity of mass, 
conservation of energy, and the impulse law, 
then, provide the following relations: 

Continuity, 


ui/Vi=u2/ V2. (1) 
Momentum, 


P\+u,?/ Vi =P2+u,?/ Vo, (2) 


1” Becker, Zeits. f. Physik 8, 321 (1922); N. A. C. A. 
Tech. Memo. 505, 506. 
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Fic. 1. The variation of properties along the Hugoniot 
curve. 


Energy, 
E,+P; Vitus?/2=E2+P2Vet+ue?/2. (3) 


The simultaneous solution of Eqs. (1) and (2) 
gives Hugoniot’s expressions for the velocities 
u, and ue: 


U1/V1=U2/V2=[(P2—P:)/(Vi-— V2) ]!. (4) 


These expressions (4) can be used to eliminate the 
velocity terms in (3) and the result is Hugoniot’s 
equation, 


3(Pi+P2)(Vi-— V2) =(E2—E)). (5) 


Because of the high temperature attained by 
combustion within the wave front, the discharged 
material may be represented as a mixture of 
perfect gases. Further, the internal energy of 
these burned gases is assumed to be a function 
of temperature only. It follows that 


P2V2=N:2RT?2 (6) 
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T2 
and E,-Ex= [ C,dT —AE, (7) 
T1 


where C, is the instantaneous specific heat at 
constant volume of the burned gases and AE 
is the energy released during constant volume 
combustion. The values of Ne, C, and AE depend 
on the progress made by the combustion reac- 
tions which may be fixed either by assumption 
or by the requirements of thermodynamic 
equilibrium. For a given initial mixture, the 
four independent equations above (exclusive of 
equilibrium equations) establish relations be- 
tween the five variables “;, u2, P2, V2 and T>. 
These relations for compression waves, in which 
Ve< Vi, are shown in Fig. 1. To reach a solution 
for the stable detonation velocity u;, another 
independent relation is required. It will be seen 
from Fig. 1 that the Hugoniot relations admit 
velocities ranging from a certain minimum to 
infinity. Chapman assumed this minimum value. 


THE JOUGUET RELATIONS 


Jouguet has shown that further relations may 
be obtained from Hugoniot’s Eq. (5), the more 
important of which are conveniently derived in 
this somewhat different way. Since the condi- 
tions ahead of the wave front are regarded as 
constant, differentiation along the Hugoniot 
curve provides the equation, 


(Vi— V2)dP2—(Pi1+P2)d V2—2dE2=0. 


Inasmuch as the composition of the burned gases 
is known, either by assumption or by the re- 
quirements of thermodynamic equilibrium, the 
entropy S2 may be introduced by the relation 
T1dS,=dE2+P2d V2 so that the above equation 
becomes , 


dS_/d Vo = [( Vi —_ V2)/2 T? | 
X[(P2—P:1)/(Vi- V2) +dP2/dV2 }. 


Now at a point where dS2/dV2=0, the deriva- 
tive dP./dV_ acquires the subscript s and the 
second derivative is 


d*?S2/d V2? = [( Vi = V2)/2 T2 |L (d?P2/d V2?) Js. 
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For a mixture of perfect gases, in which y is the 
ratio of specific heats, it follows that 


[ d?P2/d V2? |, = (d/dV2)(—y2P2/V2)s 
= 72(y2+1)P2/ V2’, 


which is always positive. Consequently, in 
waves of compression (where V2<Vj), the 
second derivative d*S2/dV.? is always positive 
and the equation 


(P2—P;)/(Vi— V2) = —(dP2/dV2)s=y2P2/V2 (8) 


defines a point of minimum entropy. By multi- 
plying both sides by V2”, this equation becomes 


V2?(P2 —P,)/( Vi-—V2)=y2P2V2 or 


U2= C2 


where C2 is the velocity of sound in the burned 
gases at temperature 72. The physical signifi- 
cance of this result is not at once clear. Jouguet 
ingeniously postulated that a sound wave, 
evidently, pursued the wave front and, by con- 
trolling the pressure P2, governed the velocity of 
detonation. He argued that the detonation 
velocity would be stable when this sound wave 
was just able to overtake the wave front, that is, 
when u2=C2. So far as is generally known, this 
physical hypothesis has never been verified. 
Acting on these ideas, Jouguet chose Eq. (8) as 
the additional independent relation required to 
effect a solution for the detonation velocity 1. 
The variation of entropy and acoustic velocity 
along the Hugoniot curve is shown in Fig. 1. 
It will be seen that the detonation velocity thus 
postulated is exactly that assumed by Chapman, 
namely, the minimum value. 

For computation purposes, Jouguet combined 
Eq. (7) with Eq. (5) and used the equation of 
state to eliminate the pressure terms in both 
(5) and (8). Then, if w= Vi/ Ve, Eq. (5) becomes 


f CT =AE—(R/2)(u—1)(NaT2+ NiTi/u), 


and Eq. (8) becomes you?—(y2+1)u+MiTi/ 
N2T2=0, which are solved simultaneously for u 
and 72. When combustion is considered to pro- 
ceed to thermodynamic equilibrium, the addi- 
tional composition variables may be eliminated 
by the like number of equilibrium equations. 
Having found yu and 7:2, the velocity of detonation 
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can be computed from u;= pu(y2N2RT>2)', which 
results from Eqs. (4) and (8). 


THE AVAILABLE ENERGY OF STEADY FLOW 


Suppose that steady flow is established within 
a rigid channel and that the temperature datum 
for available energy is To. If this flow is ideal, 
that is frictionless and reversible, the energy 
equation accounting for any changes of state at 
thermodynamic equilibrium is TdS=dE+PdV, 
and the energy equation accounting for the 
mechanical equilibrium of the forces producing 
motion is 0= VdP+udu+dL, where dL repre- 
sents high grade energy derived from the flow 
such as work derived by mechanical means."! 
The sum of these equations is 


TdS—dL =dE+d(PV)+udu. (9) 


In general, the actual flow is passing from state 1 
at section 1 to state 2 at section 2 involves a loss 
of available energy. It is possible, however, to 
measure the actual loss of available energy by 
substituting for the actual flow process, a series 
of reversible transformations from which the 
equivalent change of available energy can be 
estimated. Following the usual convention of 
signs for the energy equation, (1—7)/7T)TdS 
=dA-+dL, where (1—7)/T)TdS is the available 
portion of any heat 7dS added reversibly to 
the system, dA is the change of the available 
energy of the system, and dL is high grade energy 
as definec above. The change of the available 
energy of the system is, then dA = TdS—T dS 
—dL. This equation becomes integrable by 
virtue of (9), that is 


dA =dE+d(PV)+udu—TdS 
so 


A,—A2= (AE, + Pi Vitu,?/2—T oS) 


— (Fo+ P2V2t+ue2?/2—ToS2). (10) 


These above considerations have been studied by 
J. Willard Gibbs and others.” 


11 For the equations of flow of compressible fluids see, 

among others, Stodola, Dampf und Gas-turbinen, 1924, 
p. 31. 
12 Gibbs, Trans. Conn. Acad. 2, 382 (1873); Scientific 
Papers, Vol. 1, p. 49; Gouy, J. de physique, November 
(1889); Stodola, Zeits. d. Ver. deutsch. Ing., p. 1088 (1898); 
Jouguet, Revue de mecanique, March (1906); Darrieus, 
Revue gen. de l’elect. 27, 963 (1930); Keenan, Mech. Eng. 
54, 195 (1932). 
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AVAILABLE ENERGY CHANGES IN THE DETONA- 
TION WAVE 


Since the conditions immediately ahead and 
behind the wave front are assumed not to 
change with time, the above expression may be 
used to measure the decrease of available energy 
suffered by the gases as they stream through the 
wave front. In this case, Eq. (10), by virtue of 
Eq. (3), becomes A1—A2=T7 o(S2—S;), the de- 
crease of available energy in calories per mole of 
gas mixture passed through the wave front. The 
initial gas mixture is fed into the wave front at 
the rate of u;/V; moles per second, so that the 
time rate at which degradation proceeds is 
D=T (S2—S1)ui/ Vi calories per second. When 
the velocity “, attains its minimum (the value 
verified by experiment) Se: is likewise at its 
minimum so that D must also be minimum, that 
is, normal detonation proceeds at that speed for 
which the time rate of degradation of available 
energy is least. 

The available energy entering the wave front 
as measured with respect to the stationary gases 
ahead of the oncoming wave, that is, the avail- 
able energy passing the moving plane 1 as 
measured from the datum state represented by 
the undisturbed gas mixture at some stationary 
plane 0 ahead of the wave front, is given by 
Eq. (10) as follows 


A,—Aop=(£i+ PiVitu,?/2— ToS) 
a (Eot+Po Votuo?/2 eal T So). 


These parentheses are identical except for the 
velocity terms, namely, 7, which is the velocity 
of the wave front, and uo, which is zero since the 
gases at plane 0 are at rest; so the above expres- 
sion reduces to A,;=u;?/2+Ao, the available 
energy fed into the wave front in calories per 
mole. The time rate at which available energy 
enters the wave front, in calories per second, is, 
then, F=(u,2/2+Ao)u:/Vi, from which it fol- 
lows that normal detonation proceeds at that 
speed for which the available energy supplied to 
the wave front (as measured with respect to the 
available energy of the stationary gases ahead 
of the wave) is least. 

At any given point along the Hugoniot curve, 
the extent of the degradation (D) within the 
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wave front in comparison to the available energy 
(F) supplied may be estimated from the ratio 
of these two quantities, that is, by 

p=D/F=T )(S2—S;)/(u1?/2+Ao0). (11) 
The derivative of p with respect to V2, along the 
Hugoniot curve, is 





d*p 


To d*S» [ 1 
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dp To dS» To V; (S2—S}) dP» ———] 
dV, AidV2 A,;? 2 (V;—V2)ldV2 Vi-Ve 


At a point where dS:/dV2=0, the derivative 
dP:/dV»2 acquires the subscript s so that, when 
equated to zero, the above equation reduces to 
Jouguet’s postulate, Eq. (8). The second deriva- 
tive at this same point is 








dV2 A, dV 


For waves of compression, the several terms of 
this expression are intrinsically positive. Refer- 
ring to the denominator of the fraction within 
the large brackets, it will be observed that the 
left-hand part is, by virtue of Eq. (5), less than 
(Ex—E,). The value of the available energy 
Ao, to have significance here, should be referred 
to a datum state that is stationary with respect 
to plane 0 and thermodynamically consistent 
with the datum used in assigning a value to F. 
Under such conditions, Ao will be some frac- 
tion of E,, and thus, the denominator of the 
above fraction will be less than E2 which, on the 
T-S plane, is a smaller area than the rectangle 
represented by the numerator 72(.S.—.S,) and 
so, the fraction is greater than 1. Accordingly, 
the second derivative is negative for waves of 
compression, and p isa maximum at that point on 
the Hugoniot curve defined by Eq. (8). This is 
the point defined by Jouguet’s postulate, and the 
corresponding detonation velocity is the mini- 
mum allowed by the Hugoniot relations, the 
same that Chapman had assumed and which ex- 
periment has verified. The fact that p is a maxi- 
mum at this point on the Hugoniot curve im- 
plies that the normal detonation wave proceeds 
at that speed which allows the greatest degrada- 
tion of each unit of energy supplied to the wave. 

This result is in accord with the fundamental 
principle of degradation of energy in all natural 
processes—a principle that is so universally 
verified that when the extent of degradation 
brought about by a natural process is discovered 
to vary with the speed of the process, it would 
seem only consistent that the stable speed would 
be that which allowed the greatest degradation. 


T'o(S2—S;) | 


2(P2—P;)(Vi-— V2) +((Vi-— V2)/V1)?Ao 





In the case of detonation, such a statement of the 
principle of degradation would not only explain 
the success of Chapman’s outright assumption 
and of Jouguet’s physical postulate but would 
also provide a fundamental starting point from 
which to derive such stability relations. The 
velocity of the normal detonation wave may, 
therefore, be regarded as fixed by arguments 
within the realm of thermodynamics. 


CALCULATED RESULTS 


In order to provide a specific illustration of the 
way in which the function p identifies the stable 
speed of normal detonation from the range of 
velocities defined by the Hugoniot relations 
without recourse to either (1) the physical con- 
siderations required by Jouguet’s postulate or 
(2) the outright assumption of Chapman, the 
following calculations have been carried out with 
a mixture of 1 mole of hydrogen and 3 mole of 
oxygen; a mixture for which Dixon‘ found the 
experimental detonation velocity to be 2821 
meters per second when the initial pressure and 
temperature were 1 atmosphere and 10°C. 
The leading calculated items are listed in Table I. 


TABLE I. 








p 
(Ao =0) 


0.00000 
.10997 
.25525 
.27725 max. 
.27465 
.24196 
.20006 


uy) 
(m/sec.) 


10.09 ra) 

11.27 4516.3 
15.41 2946.1 
19.69 2820.9 min. 
21.97 2836.2 
31.48 3046.8 
44.31 3403.8 


P; 
u=Vi/V2 ~~ (atmos.) 
1.0000 
1.1097 
1.4837 
1.8553 
2.0480 
2.8129 
3.7563 
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These computations consider the dissociation 
He=2H and HxO=H2+30z; all thermal data 
being taken from Bulletin 139 of the University 
of Illinois Engineering Experiment Station.” It 
will be observed that when p is maximum, the 
corresponding velocity “; is 2820.9 meters per 


‘8 The use of the recent spectroscopic thermal data will 
entail small modifications in Table I. The merits of the 
function p, of course, are independent of the correctness of 
the underlying data. 
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second which agrees with Dixon’s experimental 
value. Referring to Eq. (11), it will be seen that 
the value of Ao will influence the magnitude of p 
without changing the position of its maximum 
on the Hugoniot curve. For the purposes of this 
example, this constant has been taken as zero. 
With appropriate values of Ao, p would be 
smaller and its maximum, although less sharp, 
would identify the same detonation velocity. 
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The Ionization Potential of Acetone Vapor 


W. ALBERT Noyes, JRr., Department of Chemistry, Brown University 
(Received April 13, 1935) 


The ionization potential of acetone vapor has been measured with a two-electrode tube. The 
value found agrees within experimental error with that predicted from spectroscopic data. 


HE absorption spectrum of acetone vapor 

has been investigated recently both in the 
quartz region and with a vacuum grating to as 
far as 850A.' It consists of a continuous absorp- 
tion in the neighborhood of 3000A with a faint 
discrete absorption on the long wave end, a 
rather complex system of bands extending from 
about 1960A to about 1800A and finally several 
other band systems obviously involving several 
different upper electron states at still shorter 
wavelengths. Among the latter bands Duncan!» 
found one series which could be represented by 
the Rydberg formula 


> = 82,767 —R/(n—0.495)?(n=3, 4, 5-+-). (1) 


No other Rydberg series which were convincing 
from the standpoint of intensity relationships 
could be found, although it was thought at one 
time that two others existed. Eq. (1) predicts an 
ionization potential of 10.20 volts and the pres- 
ent work was undertaken with the object of 
verifying this value by the electron collision 
method. 


(A) EXPERIMENTAL PROCEDURE 


The purification of the acetone has already 
been described.'* It was stored over calcium 
chloride. 

1(a) W. A. Noyes, Jr., A. B. F. Duncan and W. M. 


Manning, J. Chem. Phys. 2, 717 (1934); (b) A. B. F. 
Duncan, ibid. 3, 131 (1935). 


Fig. 1 shows a diagram of the apparatus. The 
acetone vapor was admitted through the capil- 
lary tube H from a supply bulb kept at one of the 
following temperatures: — 77° (COve-ether) ; — 33° 
(liquid ammonia); 0° (ice water). At these 
temperatures the vapor pressures of acetone are 
approximately 0.15, 7 and 67 mm, respectively. 
In general the best results were obtained at the 
lowest vapor pressure, although the rate of 
evacuation was sufficiently high to permit 
measurements to be made at the higher tempera- 
tures. At B are two diaphragms made of glass 
perforated by holes about 0.5 mm in diameter. 
The traps A and C were immersed in liquid air 
during a run and E was attached to a high speed 
diffusion pump. Since the vapor pressure of 
acetone is immeasurable on a McLeod gauge 
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at the temperature of liquid air, the main re- 
moval of the acetone was accomplished by these 
traps rather than by the diffusion pump, the 
latter serving mainly to remove uncondensable 
gases either given out from various parts of the 
apparatus or by the decomposition of acetone 
by the hot filament. That such gases were effec- 
tively removed was indicated by the fact that 
during a run, even with the acetone at 0°C, the 
McLeod gauge attached to F always indicated 
a good ‘‘sticking’’ vacuum. With this experi- 
mental procedure a unidirectional stream of 
acetone vapor passed between the filament and 
the plate. 

The electrical circuit needs little discussion. 
The resistance D permitted the potential between 
various portions of the filament and the plate to 
be utilized. It was found that the negative end 
of the filament gave the results requiring the 
smallest correction. The plate current was read 
by the galvanometer G. The resistance R 
served merely to reduce somewhat the personal 
equation, the plate voltage being calculated 
from the value of this resistance, the resistance 
of the voltmeter and the reading of the volt- 
meter. The filament was of the spiral type made 
of thoriated tungsten. The entire tube was sur- 
rounded by an electrostatic shield grounded to 
the filament. 

As calibrating gases mercury vapor and hydro- 
gen were used. For the former liquid air was 
merely removed from the traps and the source 
of acetone vapor cut off, the diffusion pumps 
being kept in operation. Electrolytic hydrogen 
was admitted to the apparatus to a pressure of a 
few hundredths of a millimeter, liquid air being 
kept on the traps. 


(B) RESULTS 


Large scale graphs were made of the 3/2 
power of the accelerating voltage against the 
galvanometer deflections. Since both from theo- 
retical reasons and because of the nature of the 
experimental procedure, the exact position of the 
change in slope cannot be ascertained with high 
precision, readings were made at 0.5 (in some 
Cases smaller) intervals and enough runs were 
made so that a statistical average would have 
some significance. 
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35 45 55 
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Fic. 2. Top curve, R=1500; middle curve, R=2000; 
bottom curve, R=0. In the lower right-hand corner a plot 
is given of successive differences in galvanometer deflection 
against the voltage for the data used in the middle curve. 


Fig. 2 shows three typical curves. The fact 
that the current increased less rapidly with 
voltage after the “‘break’”’ than before can be 
explained on the basis of the geometry of the 
tube and the electrostatic shield surrounding 
it.2 At potentials just exceeding the ionization 
potential the positive ions would probably 
contribute little toward reducing the space 
charge in the neighborhood of the filament. 

Table I presents°a summary of the uncor- 
rected results. 


TABLE I. Average uncorrected values for ionization potentials 








Prob. error 
of arith. 
mean 


F 0.03, Vv 
4 . 0.07; v 
’ 0.07; v 


No. 
runs 


Ave. 
PD. 


10.3. Vv 0.3; 
10.4; v 0.3) \ 
16.25 v 0.45 v 


Ave. dev. 
from mean 


Gas 


Acetone 67 
Mercury 15 
Hydrogen 25 





V 








The ionization potential of mercury may be 
determined directly with high precision.* In 
the case of hydrogen the value obtained from 
Rydberg series.is 15.36 volts and Bleakney has 
obtained experimentally a value in excellent 
agreement with this figure.‘ However, applica- 
tion of the Franck-Condon principle leads to the 
prediction that the most probable transition 
produced by electron bombardment would form 

2Cf. A. R. Olson and T. F. Young, Phys. Rev. 25, 58 
(1925). 


3 Cf. E. O. Lawrence, Phys. Rev. 28, 947 (1926). 
*W. Bleakney, Phys. Rev. 40, 496 (1932). 
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a hydrogen molecule ion possessing vibrational 
energy, and that usually the observed ionization 
potential should be several tenths of a volt in 
excess of the above figure, at 15.9 volts.6 This 
value will be used. The two corrections obtained 
for calibration are, therefore, 0.03 and 0.36, 
average 0).2o. 

The ionization potential of acetone corrected 
by the above figure is, therefore, 10.1 with an 
uncertainty of several tenths of a volt. 


(C) Discussion OF RESULTS 


If the 0,0 bands have been chosen properly, 
the ionization potential predicted by the Ryd- 
berg series (Eq. (1)) would be that required to 
produce a positive acetone ion possessing only 
zero-point energy. While the value obtained in 
these experiments is not very precise, the agree- 
ment with the predicted value indicates that 
such an acetone ion is produced directly by elec- 
tron bombardment. 

The electrons in the carbonyl group may be 
described to a first approximation in the same 
manner as those for a diatomic molecule. Thus 
they would be given the designation a,’7,’ 
(a7?) ._,.(a*b?)-_- where the subscripts indicate 
the atoms or the bonds with which the elec- 
trons are concerned. a and 0 are used for want 
of a better designation for the carbon-carbon 
electrons. Of the electrons associated with the 
oxygen atom two (probably o,) would be classi- 
fied as nonbonding and would have ionization 
potentials not differing much from that of the 
oxygen atom,® two (probably 7z,) would be 
classified as antibonding and would have ioniza- 
tion potentials lower than that of the oxygen 
atom. Of the four bonding electrons two should 
probably be weakly and two strongly bonding, 

5 See O. W. Richardson, Molecular Hydrogen and Its 


Spectrum, Yale University Press, 1934, pp. 1, 154ff. 
6 Cf. R. S. Mulliken, Phys. Rev. 46, 549 (1934). 
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both sets probably having ionization potentials 
higher than either the oxygen or the carbon 
atom. One would predict, therefore, three and 
possibly four different ionization potentials for 
acetone, the lowest probably corresponding to 
the removal of a z, electron which will, in reality, 
not be associated entirely with the oxygen atom 
but belong to the group as a whole. The pre- 
dicted ionization potential would be lower than 
that of oxygen (13.55 volts) and probably lower 
than that of carbon (11.22 volts).? The values 
found both spectroscopically and in the present 
paper are in accord with this prediction. Other 
possible Rydberg series might be found. 

The removal of a zo antibonding electron 
probably would not lead to the production of an 
acetone ion in a highly excited vibration state. 
The two electronic band systems at 3000 and 
1960A which do not obviously fit into Rydberg 
series are probably to be ascribed to transitions 
involving either the bonding electrons or the 
nonbonding electrons. Since, if the analysis of the 
band system at 3000A is correct, the lowest de- 
formation frequency in the upper state differs 
considerably from that in the ground state, it is 
possible that transitions involving particularly 
the bonding electrons are involved. At 1960 since 
the fundamental frequencies do not seem to be so 
profoundly modified, it is possible that the 
nonbonding electrons are concerned. Higher 
states involving these two types of electrons 
do not seem to be stable and may be responsible 
for the weak continuous absorption observed 
throughout the Schumann region as well as be- 
tween 2500 and 3000A. These suggested interpre- 
tations of the spectrum of acetone will be ex- 
amined in the light of other evidence in the 
near future. 


7Cf. R. F. Bacher and S. Goudsmit, Atomic Energy 
States, McGraw-Hill Book Co., New York, 1932, pp. 333, 
101. 
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This section will accept reports of new work, provided these 
are terse and contain few figures, and especially few halftone 
cuts. The Editorial Board will not hold itself responsible for 
opinions expressed by the correspondents. Contributions to 
this section must reach the office of the Managing Editor 


Homogeneous Dissociation of Hydrogen Molecules 
by Collision with Positive Ions 


In a series of recent articles! we have shown that positive 
ions, emitted from heated sources, and accelerated by an 
electric field produce, after reaching certain critical poten- 
tials, a sharp reinforcement of the hydrogen as well as 
nitrogen clean-up. Calculations of the energy which under 
optimal conditions can be absorbed by a given molecule 
on collisions with the ions showed, that at different critical 
potentials, corresponding to ions of different masses, the 
molecule is always able to absorb the same definite quan- 
tity of energy corresponding to some excited state of the 
molecule (12.4 v for hydrogen). 

These data were further developed in a paper from our 
laboratory published in the Reports of the Academy of 
Science (USSR), (1934),2 concerning the synthesis of 
ammonia by bombardment of a hydrogen-nitrogen mixture 
with positive ions from synthetic spodumen-sources. For 
every sort of ions used there was found a definite critical 
potential, required for the synthesis of ammonia. The 
presence of the latter was shown both by means of pressure- 
measurements, and chemically by means of Nessler’s 
reagent. The calculation of energy that could be absorbed 
inthe optimum case by the molecules of the mixture shows, 
that for all sorts of ions used (Li+, Nat, K*+, Rb*t, Cs*) 
the nitrogen molecule was able to take up about 22 volts. 
Fig. 1 illustrates the sharpness of the effect. All these 
investigations without exception indicate the possibility of 
sharp chemical effects resulting, at critical potentials, from 
collisions of molecules with ions. The observation of the 
effects described requires very careful control of experi- 
mental conditions. These were described rather briefly 
in the foregoing articles, a circumstance which was perhaps 
the principle cause of some misunderstandings arising from 
attempts of Mitchell*.and of Kunsman and Nelson‘ to 
repeat our experiments. These authors, on the basis of 
their negative results, deny the correctness of our explana- 
tion of the data of Ho. 

In order to indicate possible causes of this disagreement, 
we shall describe some of the chief experimental conditions 
without which the observation of the observed effects is 
very difficult. 

It is very important to keep the level of liquid air in the 
Dewar flask constant, because the smallest oscillations of 
the latter sharply change the rate of pressure fall (because 
of the change in area of the cooled surface). We were forced 
to build an apparatus for automatic control of the con- 








not later than the 15th of the month preceding that of the 
issue in which the letter is to appear. No proof will be sent 
to the authors. The usual publication charge ($3.00 per 
page) will not be made and no reprints will be furnished free. 


stancy of the liquid air level, and for the enlargement of 
the observed effect, to keep the whole reaction vessel in a 
Dewar flask. It seems that these two conditions were not 
fulfilled in the work of Kunsman and Nelson. The form 
of the collector in the apparatus of Kunsman and Nelson 
(Fig. 1 of their article), indicates that hydrogen atoms had 
greater chance to recombine on the surface of the collector, 
than to arrive into the tube end cooled by liquid air, a 
circumstance which also should diminish the effect ob- 
served. 

It is further necessary to take into account that at 
pressures of the order of some tenths of a mm Hg the mean 
free path is 2-10 mm. Therefore, to secure a suitable 
probability of collisions of an ion with an H»2 molecule the 
distance grid-cathode must be taken of the order of 4-5 
cm. A smaller distance does not prevent thermal dissoci- 
ation of He, but it may diminish strongly the efficiency of 
the ions. From Fig. 1 of the article of Kunsman and Nelson 
one can assume, that in their arrangement this condition 
was not fulfilled. 

Further, before every experiment we had to outgas the 
ion source (by heating it to 750-800°C) and the vessel (by 
heating it to 350°C) by constant pumping for several 
hours. Only after our Pirani gauge no longer registered any 
pressure increase on heating the ion source could we begin 
our experiment. In the article of Kunsman and Nelson 
nothing is said about outgasing before experiments. 
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Fig. 1. The data on the curves indicate the voltages used. 


433 





a nme is mma TCS MES 


LETTERS TO 


\ 


\ 


40 60 80 100 6120 
Time (min.) 


Fig. 2. 


Finally to discover the described effect we had to choose 
the conditions (gas pressure, temperature of the source) 
so as to make the thermal dissociation on the source and 
the corresponding clean-up a minimum in comparison with 
the additional homogeneous effect. One must also re- 
member, that the clean-up process of H atoms, generated 
only by thermal dissociation does not follow a linear, but 
rather a logarithmic law. Therefore, to find the additional 
clean-up due to the action of ions, and to calculate the 
efficiency of an ion, we must plot log against ¢, and not 
simply p against t. In doing so with a part of curve (I) in 
Fig. 3 of the article of Kunsman and Nelson we get the 
curve, shown on Fig. 2 of the present article. Here we can 
see a rather marked increase of the clean-up rate at the 
time the field was turned on, and a slowing up of the clean- 
up at the time the field was turned off. Probably in this 
experiment of Kunsman and Nelson there happened to 
exist a sufficient coincidence of the required experimental 
conditions. 

It appears therefore, that the experiments of Kunsman 
and Nelson do not entirely contradict our explanation of 
the effects we have observed, but on the contrary, seem, 
in part, to prove the correctness of our experiments. 

A. SCHECHTER 

Institute for Chemical Physics, 

Leningrad, 
May 17, 1935. 

1 Leipunsky and Schechter, Zeits f. Physik 59, 857 (1930); Schechter, 
Zeits f. Physik 75, 671 (1932); Semenoff and Schechter, Nature 126, 
436 (1930). 

2 F. Fedoroff, Mocan, Roginsky and Schechter, Bildung von NH:; 
durch Stoss positiver lonen, Comptes rendus de l'Académie de Sciences 
de l'URSS, 1934, p. 367. 


3 Mitchel, J. Frank. Inst. 210, 269 (1935). 
4 Kunsman and Nelson, J. Chem. Phys. 2, 752 (1934). 
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Index of Refraction of HCl from 1 to 10u 


In view of the recent discussion in the literature con- 
cerning the effective charge of HCI,':? we believe it advis- 
able to communicate immediately the results we have 
obtained from dispersion measurements between 1 and 
10u. The dispersion curve is matched best by using a value 
of 1.00-10-* e.s.u. for the effective charge of the rotator- 
vibrator, and a value of 1.18-10-'8 e.s.u. for the electric 
moment of the rotator. The values for the effective charge 
and the electric moment are not mutually dependent to 
any great extent since the contribution of the pure rotation 
in the region near the rotation-vibration band is small, and 
vice versa. The uncertainty in the value of the effective 
charge of the rotator-vibrator is less than five percent. 
The value of the effective charge for the rotator-vibrator 
is in good agreement with Bourgin’s’ result obtained from 
measurements of the absorption intensities in the band at 
3.46u. 
Beyond about 4.5u, the contribution of the rotator to 
the index of refraction exceeds that of the rotator-vibrator, 
and at 10u the effect of pure rotation is nearly as great as 
that of the rotator-vibrator as near the 3.46y absorption 
band as we could accurately measure it. The measure- 
ments between 5 and 10y thus afford an excellent check of 
the negative terms of the Kramers dispersion formula, 
since the majority of the molecules are in excited states 
and the number in each state can be accurately computed. 
If the negative terms are omitted and a reasonable value 
of the electric moment is used, one obtains by computation 
a contribution due to pure rotation which is more than 
twice as large as that obtained by experiment. 
The value of the electric moment required to match the 
dispersion measurements is slightly larger than the value 
obtained by Zahn‘ (1.034-10-'8 e.s.u.) from dielectric 
constant measurements. This indicates that the low inten- 
sities obtained by Badger® and by Czerny® for the pure 
rotation spectrum are to be explained by the extreme dif- 
ficulty of making absolute intensity measurements, 
especially in the far infrared. A complete report of the 
work will be published shortly. 
R. ROLLEFSON 
A. H. ROLLEFSON 

Physical Laboratory, 

The University of Wisconsin, 
Madison, Wisconsin, 
May 18, 1935. 

E. C. Kemble, J. Chem. Phys. 3, 317 (1935). 

R. S. Mulliken, J. Chem. Phys. 2, 400, 712 (1934). 

D. G. Bourgin, Phys. Rev. 32, 237 (1928). 

= T. Zahn, Phys. Rev. 24, 400 (1924). 


M. Badger, Proc. Nat. Acad. Sci. 13, 408 (1926). 
M. Czerny, Zeits. f. Physik 34, 227 (1926). 
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An Attempted Concentration of the Heavy 
Nitrogen Isotope 


Urey and Greiff! have theoretically demonstrated that 
isotopic exchange reactions might be used for the separation 
or concentration of somé of the isotopes of the lighter 
elements. The isotopic reaction, 


N“H;+N"“H,OH=N"“H;+N®H,OH 
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TABLE I.* 








Ratio 


Sample NUN : NUNS 








* Between analyses, conditions in the mass spectrograph changed so 
that the background correction for the first three samples was not the 
same as for the last five. Thus the first three samples are comparable, 
as are the last five. 


was studied with the hope of concentrating the heavy 
nitrogen isotope. A twenty-five foot packed column, which 
was erected by Dr. W. G. Brown for the study of other 
exchange reactions, was used for this purpose. Gaseous 
ammonia was passed into the bottom of the column at a 
constant rate and water in at the top, at sucha rate that the 
through-put of ammonia was very small. After six hours of 
continuous operation with as large a reflux ratio as possible 
to maintain, a small sample of forward flow ammonia was 
collected for analysis. This is designated as sample No. 3 
and is to be compared with sample No. 1, a sample of tank 
ammonia. 

The ammonia molecule was not satisfactory for analysis 
in the mass spectrograph, so all the samples were converted 
to nitrogen by the reaction of ammonia with sodium hypo- 
bromite. The results of the analyses are recorded in Table I, 
together with the analyses of other samples described below. 

It isevident that the distribution of N" is slightly favored 
in the NH4,OH molecule; however, either the fractionation 
factor for this reaction is almost unity, making this reaction 
most impractical for the concentration of heavy nitrogen, 
or the column did not contain nearly as many theoretical 
plates as we had anticipated. 

Since a large Dewar flask filled with anhydrous ammonia 
served as a source of ammonia for the above experiment, it 
was possible at the same time to determine the ratio of 
vapor pressures of the molecules NH; and NH;3. The 
residue from this distillation is sample No. 4. The factor by 
which the sample was taken down in this experiment is 
somewhat uncertain. It may be as small as 76 and certainly 
not any larger than 800. If we use the first figure, a value of 
1.012 is obtained for the ratio of vapor pressures, and using 
the second the ratio becomes 1.0082. 

The experiment was repeated taking care to measure the 
volume of the residue more accurately. The factor by which 
the sample was reduced in this case was approximately 
32,800. The residue is designated as sample No. 8 and is to 
be compared with No. 7, a sample of tank ammonia. The 
ratio calculated from data obtained in this experiment is 

1.0052. This value is not in good agreement with those 
obtained above. In this experiment, a four-liter Dewar of 
liquid ammonia was evaporated under atmospheric pressure 
to4 liters of ammonia gas. This gas was then condensed ina 
small trap and again evaporated till only the trap filled 
with gas remained. It seems probable that when the liquid 
residue becomes small, spattering and lack of stirring defeat 
the concentration that should result from continued dis- 


tillation. In any case, the value 1.0052 is the minimum to 
be expected for the ratio of vapor pressures, and it is possible 
that the values obtained in the first experiment are more 
nearly correct. 

The mean deviations recorded in Table I are somewhat 
larger than those actually calculated. We have been some- 
what liberal in these estimates, since they are calculated 
from only five or six determinations in each case. The ratio 
of vapor pressures given as 1.0052 may be in error by about 
0.0013 because of the errors in analysis. Incidentally it may 
be remarked that these values for the nitrogen isotope ratio 
are in excellent agreement with that of Vaughan, Williams 
and Tate? who found N* : N®=265+8. 

A ten-foot vacuum-jacketed packed column one inch in 
diameter, was constructed for the distillation of liquid 
ammonia. A dry ice condenser was installed at the top of 
the column and the ammonia was boiled at the rate of two 
liters per hour. After several hours operation at total reflux 
a sample (No. 9) was taken by passing a slow stream of 
inert gas through the condenser. The analysis shows the 
column to have approximately six theoretical plates. This 
is less than the number expected, but is in accord with the 
findings of Cryder, Fenske and co-workers* on the distilla- 
tion of water with packed columns of this type. 

Samples of ammonia were obtained from the commercial 
stills of the Mathieson Alkali Works and analyzed to 
ascertain whether there might be a commercial source of 
ammonia containing an increased concentration of N™H3. 
Sample No. 5 is the starting material and No. 6 is a sample 
of the 3000-lb. residue obtained from 70,000 Ibs. of am- 
monia which was distilled in four batches, each batch 
being added to the preceding residue. It is apparent that 
the distillation is carried out so rapidly at such a small 
reflux ratio that the concentration is almost negligible. 

The vapor pressure ratio N“H;-N"H;, although not 
particularly encouraging, indicates that the distillation of 
liquid ammonia may be a means of obtaining increased con- 
centrations of N“H;. The difference between these vapor 
pressures is similar to that in the case of the H,O"* and 
H,0'8 waters.* 

We are indebted to Professor H. C. Urey for suggesting 
this problem and for suggestions offered during the course 
of the work, and to Professor Walker Bleakney for much 
helpful advice and the use of his apparatus in making the 
analyses. We also wish to thank the Mathieson Alkali 
Works, Inc. for sending us samples of ammonia from their 
large stills. 

M. H. Wank 
J. F. HurrmMan* 
Columbia University, 
New York, N. Y.; 
J. A. Hrppece, Jr. 

Princeton University, 

Princeton, New Jersey, 
June 13, 1935. 


* Indebted to E. K. Adams Fellowship Fund for financial assistance 
received during the course of this work. 

1H. C. Urey and L. J. Greiff, J. Am. Chem. Soc. 57, 321 (1935). 

ase L. Vaughan, J. H. Williams and J. T. Tate, Phys. Rev. 46, 327 
(1934). 

3D. S. Cryder, M. R. Fenske, D. Quiggle, C. O. Tongberg and A. R. 
Lux. Presented at New York meeting of Am. Chem. Soc. (1935). 

4M. H. Wahl and H. C. Urey, J. Chem. Phys. 3, 411 (1935). 
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Concerning Hydrocarbon Bond Additivity 


The data of Kistiakowsky, Ruhoff, Smith and Vaughan! 
on the heats of hydrogenation of some olefins furnish a test 
for the corrections proposed by the author? to obtain bond 
additivity. The two corrections proposed were (1) inclusion 
of the zero-point vibrational energy of each molecule in its 
total potential energy; (2) an additional interaction be- 
tween those hydrogen atoms attached to the same carbon 
atom. The latter is admittedly, at present, an empirical 
correction factor. 

Consider the reaction involved in hydrogenation; 
namely, C,H2,+H2~>C,H2n+2. The sum total effect is to 
dissociate 1 C=C bond and one molecule of He, and to form 
1 C—C and 2 C—H bonds. But, in addition, two further 
changes occur; first, the sum of the zero-point energies of 
the reactants differs from that of the product. Second, the 
H---H interaction is greater in the saturated hydrocarbon 
and, in addition, depends specifically on the position of the 
double bond in the olefin. For example, ethylene has 2 
H- --H interactions, ethane has 6, the change amounting 
to 12 Kcal.; propylene has 4 H---H, propane has 7, the 
change amounting to 9 Kcal.; butene (2) has 6 H---H, 
normal butane has 8, the change being 6 Kcal. 

In general, then, the heat of hydrogenation of C,Hen 
may be written as follows: 


AH = {Ho-c+Aoc-4 + Ayo (n...8) + Cpt 42} 
— {Ho-c+109.3 + cnte,} +** (1) 


where Ho-c, Ho-y, Hoec are the bond strengths of the 
bonds noted in subscript; Hsy...q is the change in the 
H---H interaction which is realized on hydrogenation; 
e denotes the zero point energy; 109.3 Kcal. is the heat of 
dissociation of Hz plus its vibrational zero-point energy. 

The zero-point energies of the saturated hydrocarbons 
are obtained from the frequency date of Kohlrausch and 
Koppel.’ The frequencies for ethylene were taken from the 
analysis of Mecke.t No Raman nor infrared data are avail- 
able at present from which one may estimate the frequen- 
cies for propene and the butenes. However, data do exist for 
the aldehydes and the ketones,® which furnish the means for 
a comparison. The skeleton vibrations in propene are ap- 
proximated by acetaldehyde; those in butene (1) by 
propylaldehyde; those in isobutene by acetone. From 
analogy with the cis and trans modifications of dichlor- 
ethylene, butene (2, cis) should have a frequency distribu- 
tion more like isobutene, while butene (2, trams) is more like 
the chain found in butene (1). Since we are expressing 
energies in kilocalories, this estimation is not very sensitive 
to the value of the frequency chosen, provided the total 
number of vibrational degrees of freedom are satisfied. 
For example, the stretching vibrations of the hydrogens 
alone contribute a little more than half the zero-point 
energy. Since the number of degrees of freedom for the 
stretching frequencies of C—H bonds equals the number of 
hydrogen atoms, and since these frequencies are all in the 
region 2800 to 3000 wave numbers, only a very small 
uncertainty is entailed if the assignment is in error. 

Using the values for the C—H bond and the H---H 
interaction from reference 2, and the observed heat of 
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hydrogenation for ethylene the C=C bond was obtained of 
138.5 Kceal. The remaining heats of hydrogenation were 
then calculated as given in Table I. 


TABLE I. 








AH AH (0°K) 


AgH-++H €CpHon €CnHen+s 
ie ee (Calc.) (Obs.) 





Ethylene 33.4 
Propylene 50.5 
Butene (1) 67.0 
Isobutene 66.0 
Butene (2, cis) 66.0 
Butene (2, trans) 67.0 


28.4 
28.6 
26.8 
26.6 
25.6 


NO 


00 00 00 90 me ohm 
Se Ser) 








The value of 138.6 for the C=C bond is higher than that 
calculated from the heat of combustion data, but is more 
probable because of the uncertainties in the combustion 
data for ethylene. The use of this value for the C=C bond 
does not change the conclusions arrived at in the first 
communication of the author.” 

V. Derrz* 

Urbana, Illinois, 

June 4, 1935. 

* National Research Fellow. 

1 Kistiakowsky, Romeyn, Smith and Vaughan, J. Am. Chem. Soc. 
57, 65 (1935); Kistiakowsky, Romeyn, Smith and Vaughan, J. Am. 
Chem. Soc. 57, 876 (1935). 

2V. Deitz, J. Chem. Phys. 3, 58 (1935). 

3 Kohlrausch and Koppel, Zeits. f. physik. Chemie B26, 209 (1934). 


4 Mecke, Zeits. f. Physik 64, 173 (1930). 
5 Kohlrausch and Képpel, Zeits. f. physik. Chemie B24, 370 (1933). 


Photodissociation of Polyatomic Molecules in the 
Schumann Ultraviolet 


Very valuable information concerning the way a molecule 
is disrupted on absorption of light may be obtained from the 


absorption spectrum; the interpretation is, however, 
seldom definite, especially in the case of broad continuous 
regions of absorption exhibited, as a rule, by polyatomic 
molecules. 

A direct observation of the photodissociation process 
becomes possible when the magnitude of the absorbed 
quantum is large enough to disrupt a molecule into excited 
radicals which can be then detected by their emission. From 
the study of this emission various important details regard- 
ing the mechanism of the disruption can be obtained, as was 
shown in a series of papers from this laboratory.! 

For the extension of this method to polyatomic organic 
molecules, radiation in the Schumann ultraviolet must be 
used, since large energies are required to disrupt bonds in 
these molecules with simultaneous excitation of the 
radicals set free. 

During some researches on gaseous photochemical re- 
actions in the Schumann ultraviolet,? we observed that in 
the vapors of H2O, CH;OH and CH;CN, the bands belong- 
ing to OH and CN radicals are emitted under the action of 
this short radiation.* 

The work has been lately extended to organic acids, 
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ammonia and methylamine. In all these experiments the 
vapor of the substance studied was subjected to intense 
illumination through a thin fluorite window from a power- 
ful hydrogen discharge tube, or a spark in an argon at- 
mosphere. The shortest acting wavelength was 1300A. By 
interposing a thin film of fused quartz, wavelengths shorter 
than 1500A could be easily cut out and the range of active 
wavelengths thus approximately determined. In nearly all 
cases the observed emission disappeared when the quartz 
film was inserted in the path of the exciting light. 

The results so far obtained are summarized in the follow- 
ing table, where under D+E£ are given the threshold 
energies needed to produce the dissociation process repre- 
sented in the preceding column; the wavelengths corre- 
sponding to these energies are given in Table I under Xo. 


TABLE I. 








Process: 


H +OH* 
CH;s+OH* 
C2Hs +OH* 
HCO+OH* 

CHsCO+OH* 
CH3sCN CH3s+CN* 
NHs H +NH2* 

Ie I+I* 


D+E (kcal.) 


115 +92 =207 
90 +92 = 182 
90 +92 = 182 
90 +92 = 182 
90 +92 = 182 
105 +73 =178 
117 +55 =172 
35 +159 = 194 


Substance: 





H20 
CH;0H 
C:H;sOH 
HCOOH 
CHsCOOH 








The hydroxyl radicals were detected by the emission of 
the known band 3062A (92 kcal.), the cyanogen radical by 
the characteristic bands begimming with the 3883A (73 
kcal.) one. In ammonia and methyl amine vapors a visible 
emission appears under the action of Schumann radiation, 
consisting of a broad band lying in the green-yellow part of 
the spectrum (probably a-bands of ammonia). This emis- 
sion which is now being studied with larger dispersion, 
presumably belongs to an excited NHg2 radical liberated 
from the ammonia molecule in the dissociation process 
represented above. 

Under the action of these short wavelengths iodine 
vapor emits the atomic line 2062A of iodine. This points to 
a photodissociation of the iodine molecule with the excita- 


tion of one of the iodine atoms to the 6s ?P3,. state, from 


which the emission of this line occurs. 

The direct observation of free radicals with the help of 
their emission opens the way to a spectral study of the 
kinetics of various reactions they can undergo in collisions 
with other molecules. We are therefore studying now the 
quenching of this emission, as well as the efficiency of vibra- 
tion energy transfers from these excited radicals to various 
other molecules. A full account of the whole work will be 
given at another place. 

A. TERENIN 
H. NEUJMIN 
Photochemical Laboratory of the Optical Institute, 
Leningrad 53, USSR, 
May 25, 1935. 

1A. Terenin, Zeits f. Physik 44, 713 (1927); cf. also; Physik. Zeits. d. 
Sowjetunion 2, Heft 4-5 (1932); A. Terenin, Nature 135, 543 (1935); 
Comptes rendus Acad. Leningrad, 1935, 485. 

*H. Neujmin and B. Popov, Zeits. f. physik. Chemie B27, 15 (1934); 
see also forthcoming paper of B. Popov in Acta phys. chim. URSS on 


the reaction O’ ++CO. 
+A. Terenin and H. Neujmin, Nature 134, 255 (1934). 


THE EDITOR 


Isotopic Exchange Reactions with Iodine 


The reaction RX+I-~=RI+X~- (R=organic radical, 
X =Cl, Br, I) has been the subject of several well-known 
investigations.. When X =I the only possibility of ob- 
serving the reaction has been to make R an optically 
active radical, and to follow the racemization occurring 
in the reaction. This prevents the use of the method in the 
case of the simple alkyl radicals (besides involving an 
assumption concerning the connection between the race- 
mization and the reaction mechanism). The use of the 
radioactive iodine isotope provides a method free from any 
theoretical assumption and of general applicability. Using 
a Geiger-Miiller counter, we have observed some of these 
isotopic exchange reactions. 

Sodium iodide and iodine were activated in sealed tubes 
sunk in a paraffin block by bombardment with neutrons 
from a source of 100 millicuries of radon sealed up with 
beryllium powder. 

We will not consider here the effect of momentum of the 
neutron absorbed by the iodine nucleus and of the gamma- 
ray emitted, upon the form (ion, atom or molecule) in 
which the radio-iodine is obtained when Na*I~ or I: are 
activated, because the question is complex. But on the 
strength of the two following results we conclude that when 
activated (solid) Na*I~ is dissolved in water, the solution 
contains the radio-iodine in the form of RI~ ions. Activated 
NatI~ was dissolved in water, and inactive iodine added. 
On extracting the iodine by CS» the activity was found to 
be divided between the iodide and the iodine. Similarly, 
starting with activated iodine and inactive NatI~ the 
activity was shared. 

We have examined two series of aliphatic iodine com- 
pounds: 


CH; 


CH.I, CH;.CH2.CH:.I, CHz=CH.CH2I, 

CH; 

CH;.CH2.I, CHsl CH;I, CHelI2, CHI;. 

The procedure for the exchange experiments was dictated 
by the short half-life period of radio-iodine (ca. 24 minutes). 
Activated sodium iodide was dissolved (ca. }.N) in alcohol 
with a few percent of water. To this solution the alkyl 
iodide (ca. 2N) was added. After about 1 minute the 
organic iodide was precipitated by excess water, separated 
from the water layer and repeatedly washed. For the 
counter measurements AgI was prepared from both the 
organic iodide and the Nal solution by adding respectively 
alcoholic or aqueous silver nitrate. In the case of methylene 
iodide the iodine was first liberated by boiling with con- 
centrated sulphuric acid. The iodoform was dissolved in 
acetone (containing the activated Nal) and after pre- 
cipitation tested directly. 

With CH;I and CH2=CH.CH.I, in one minute at room 
temperature and under the conditions of solvent and con- 
centration stated, the exchange took place to a consider- 
able extent (>50 percent). With the other substances the 
exchange was negligible or small (certainly less than 10 
percent). 


and 
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Using the ordinary expression for the absolute rate of a 
bimolecular reaction, we find for the interchange reaction 
CH;I+I~ in alcohol an activation energy of 16.5 Kcal. 
as an upper limit. This is of interest in connection with the 
recent theoretical discussion! of such reactions by Ogg and 
Polanyi. 

F. JULIUSBURGER 
B. TOPLEY 
J. WEIss 
The Sir William Ramsay Laboratories 
of Inorganic and Physical Chemistry, 
University College, 
Gower Street, 
London, W.C. 1, 
May 17, 1935. 


1 See Ogg and Polanyi, Trans. Faraday Soc. 31, 604 (1935). 


The Difference in the Heats of Formation of the Two 
Isomers of Butane. An Estimate of the Heat of Forma- 
tion of Tetramethylmethane (Neopentane)! 


Data recently obtained in this laboratory on the heat 
of combustion of isobutane can be combined with similar 
data on normal butane? to yield a reliable value for the 
difference in the energies of formation of the two isomers. 

The calorimetric experiments on both normal butane and 
isobutane were made with the same apparatus and pro- 
cedure, and both hydrocarbons used in the measurements 
were of established purity. For these reasons, the reliability 
of the relatively small difference in the two values of the 
heats of combustion will be determined by the ‘‘precision”’ 
uncertainty of the calorimetric combustion experiments, 
rather than by the ‘‘accuracy”’ uncertainty. In the two 
series of calorimetric combustion experiments, the ‘‘pre- 
cision”’ uncertainties were +0.12 and +0.09 kilocalorie per 
mole for normal butane and isobutane, respectively. Taking 
the resulting uncertainty in the difference of the two heats 
of combustion as the square root of the sum of the squares 
of the foregoing quantities, the data yield, at 25°C and 1 
atmosphere, for the reaction, 


n— C4Ho(gas) = iso — C4Ho(gas), 
AHo9s.1 = —1.63+0.15 kilocalories per mole. 


Since for these isomers, at room temperature and a pressure 
of 1 atmosphere, the difference in the values of PV, and of 
E(p=1)—E(p=0), is not significant, one may write for the 
same reaction 
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AH 295.1 = E%293.1 = —1.63+0.15 kilocalories per mole. 


AH and AE are, respectively, the increments in the heat 
content and in the internal or intrinsic energy. 

The value of the heat of combustion of isobutane (tri- 
methylmethane) can be used in conjunction with the values 
already reported*: * for methane, ethane (monomethyl- 
methane), and propane (dimethylmethane) to deduce a 
fairly reliable value for the heat of combustion of tetra- 
methylemethane (neopentane). By plotting (Qc—A)/n 
(where Qc is the heat of combustion, A is a constant whose 
value was taken as 52.00, and is the number of carbon 
atoms in the molecule) against 1, one obtains a smooth 
curve which can be extrapolated to »=5 with some 
confidence. Assigning a reasonable limit of error, one ob- 
tains for the heat of combustion of gaseous tetramethyl- 
methane, C(CH3;)4, in oxygen to form gaseous carbon 
dioxide and liquid water, at 25°C and a constant pressure of 
1 atmosphere, the value 840.4+1.0 kilocalories per mole. 
This value is about 6 kilocalories per mole less than the one 
experimental value extant in the literature, that of 
Thomsen‘ reported in 1886. Thomsen’s tetramethylme- 
thane was impure, his samples having the empirical 
formulas of C;Hi.64 and C;Hu.7s, instead of C;Hi2. This 
thermochemical constant will be redetermined experi- 
mentally as soon as material of the necessary purity is 
obtained. 

Combination of the foregoing value for tetramethyl- 
methane with the values recently reported for normal 
pentane,” * yields the following: 


n— C;Hi2(gas) = C(CHs)4(gas), 

AH 295.1 = AE 293.1 = —4.9+1.0 kilocalories per mole; 
5C(c, B graphite) +6H2(gas) = C(CHs)4(gas), 

AH 293.1 = —40.7 +1.2 kilocalories per mole; 
5C(c, diamond) +6H2(gas) = C(CHs)4(gas), 

AH 93.) = —41.8+1.2 kilocalories per mole. 


FREDERICK D, RossiNI 
National Bureau of Standards, 
Washington, D. C., 
June 10, 1935. 


1 Publication approved by the Director of the National Bureau of 
Standards, U. S. Department of Commerce. 

2 Rossini, Bur. Standards J. Research 12, 735 (1934). 

3 Rossini, J. Research Nat. Bur. Standards 13, 21 (1934). 

4Thomsen, Thermochemische Untersuchungen, Vol. IV 
Leipzig, 1886). 
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